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STUDIES  OF  GAS-PHASE  ION/MOLECULE  REACTIONS  IN  RELATION  TO 
A  PROPOSED  IONIC  MECHANISM  OF  SOOT  FORMATION 

By 

FEZA  OZTORK 

August.  1988 


Chairman:  John  R.  Eyler 
Major  Department:  Chemistry 


The  reactions  of  small  hydrocarbon  ions  such  as  C3H3"'", 
C5H3*,  C5H5*,  and  CrjHj'*'    with  acetylene  and  dlacetylene  have 
been  Investigated  using  Fourier  Transform  Ion  Cyclotron 
Resonance  (FTICR)  mass  spectrometry  to  provide  qualitative 
and  quantitative  Information  about  a  proposed  Ionic 
mechanism  of  soot  formation.   Ion/molecule  reaction  pathways 
and  rate  coefficients  have  been  determined  for  several 
Isomers  of  each  of  the  Ions  listed  above,  formed  from  a 
variety  of  precursors,  and  reacting  with  precursors, 
acetylene  and  dlacetylene.   Further  understanding  of 
reaction  mechanisms  of  some  of  the  reactions  proposed 
experimentally  was  obtained  from  kinetic  modeling  studies. 

All  of  the  small  hydrocarbon  ions  studied  reacted  with 
dlacetylene  extensively  while  almost  no  production  of  larger 
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lons  was  observed  with  acetylene.   Linear  C3H3'*"  was  formed 
by  charge  transfer  Ionization  of  propargyl  Iodide  using  Xe''', 
and  was  found  to  Isomerlze  to  the  cyclic  form  of  C3H3'*^  in 
reactions  with  both  acetylene  and  dlacetylene.   The 
isomerlzatlon  was  shown  to  take  place  via  a  long-lived 
CgHg"^   complex  by  Isotope  exchange  reactions  between  linear 
C3H3'*^  and  deuterated  acetylene.   The  reaction  rate 
coefficients  for  the  reaction  of  C3H3'*"  with  deuterated 
acetylene  and  dlacetylene  were  (4.5  +  1.9)  x  10"^^  cm^/s  and 
(1.4  +  0.7)  X  10"^  cm^/s,  respectively. 

While  different  structures  could  be  attributed  to 
^3^3^.  ^5113'*^,  and  C^H^^  Ions  produced  from  different 
precursors  on  the  basis  of  reactivity,  this  could  not  be 
done  for  C5H5'*^  Ions,  whose  reactivity  with  acetylene  and 
dlacetylene  was  similar  within  experimental  error  regardless 
of  precursor.   Presence  of  two  structural  Isomers  of  C5H3''" 
was  determined  using  different  precursors  and  several 
structures  for  these  Isomers  were  proposed.   The  reactive 
Isomer  was  observed  to  react  with  dlacetylene  with  a  rate 
constant  of  (5.6  +  1.7)  x  10"-'^  cm^/s.   Rate  coefficients 
for  the  Ion/molecule  reactions  of  CrjE,j'*'    were  found  to  be 
different  when  It  was  formed  from  various  precursors,  which 
Implied  the  Involvement  of  different  C^Hy*  structures.  All 
experimental  results  are  discussed  in  relation  to 
theoretical  work  involving  structures  of  the  Ions  studied 
and  a  proposed  ionic  route  to  soot  formation. 


vill 


CHAPTER  1 
INTRODUCTION 


Soot  Is  a  coBbustion  product  which  has  undesirable 
effects  both  on  human  health  and  on  the  efficiency  of  fuel- 
powered  engines  due  to  Increased  heat  transfer  to  critical 
engine  components.   Considering  the  fact  that  the  use  of 
alternative  fuels  such  as  shale  or  coal-derived  gas  in  the  - 
near  future  would  lead  to  Increased  soot  emissions,  search 
for  an  effective  way  to  reduce  soot  has  been  going  on  for  a 
long  time.   The  design  and  development  of  advanced  engines 

/•'       to  achieve  soot  suppression  requires  a  sufficient 

•  ■  > 

fundamental  understanding  of  the  mechanisms  governing  soot 
,•  '      formation.   Soot  formation  in  hydrocarbon  flames  involves  a 
succession  of  fast  processes  that  occur  within  a  few 
milliseconds  during  the  combustion  of  the  hydrocarbon. 
i'       Although  many  investigations  have  been  undertaken  in 

relation  to  the  process  of  soot  formation  in  fuel-rich 
:^:        hydrocarbon  flames,  understanding  of  this  process  is  still 
limited.   Experimental  techniques  such  as  mass  spectrometry 
(Calcote,  1963;  Mlchaud  et  al.,  1981),  laser-light 
scattering  (Kent  et  al . ,  1981;  D'Alessio  et  al . .  1977), 
laser  Induced  ionization  (Smith  and  Mallard,  1981),  and 

l^/i,'.  y  laser  Induced  fluorescence  (DiLorenzo  et  al  .  ,  1981)  have 

1/  ■  ■    ' 
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been  used  to  study  the  molecular  species  and  particles  In 
flames  . 

Soot  collected  from  flames  consists  of  chainlike 
aggregates  of  spherical  units  having  diameters  of  10-50  nm 
with  a  carbon/hydrogen  atom  ratio  In  the  range  8:1-12:1 
(Calcote,  1981).   Three  distinct  steps  of  soot  formation 
have  been  recognized  over  many  years  of  research  (Calcote, 
1981) : 

1.  Nucleatlon-  a  phase  change  from  molecular  species  where 
chemical  reactions  dominate  to  Incipient  soot  particles 
where  physical  processes  dominate; 

2.  Growth  to  spherical  particles  of  10-50  nm  in  diameter; 
and 

3.  Aggregation  or  agglomeration  of  the  spherical  units  to 
form  chains. 

Among  these,  the  first  step,  nucleatlon,  is  the  least 
understood  in  the  process.   Although  many  diverse  theories 
have  been  put  forward  to  explain  the  nucleatlon  process, 
only  two  of  these  have  received  quantitative  support: 

1.  Mechanisms  involving  neutral  and  free  radical  species; 
and 

2.  Ionic  mechanisms. 


Neutral  and  Free  Radical  Mechanisms.   A  large  number  of 
neutral  and  free  radical  mechanisms  have  been  proposed  to 
describe  the  transition  from  molecular  species  to  soot 


(Palmer  and  Cullls,  1965;  Lahaye  and  Prado ,  1978;  Gaydon  and 
Wolfhard,  1979). 

One  of  the  proposed  paths  Involves  polyacetylenes 
where  acetylene  reacts  to  form  higher  species  (Bonne  et  al . , 
1965)  . 


^2^2 


^2*^2 


C4«2 


l-H 

I 
•>  C4H3  (+H) 


^6'^2 


^2^*2 


>  CgHg  (+H2) 


^2^2 


^8^2 

l-H 

I 

>   CgHg   (+H2) 


It  was  suggested  that  the  continuation  of  this  series 
leads  to  larger  radicals  which  further  react  with  each  other 
and  with  higher  polyacetylenes  forming  even  larger 
molecules.   However,  because  this  mechanism  leads  to 
formation  of  a  giant  chain  molecule  Instead  of  a  polycycllc 
carbon  structure,  It  was  not  accepted  as  being  responsible 
for  soot  formation  (Homann,  1967).   To  account  for  the 
problem  of  ring  formation,  Homann  and  Wagner  (1967,  1968) 
suggested  that  radicals  such  as  CgH  attack  the 
polyacetylenlc  species  forming  side  chains  which  lead  to 
branched  polymerization,  and  presumably  ring  closures.   The 
experimental  support  of  this  scheme  comes  from  the 
observation  that  In  the  oxidation  region  of  the  flame,  the 
concentration  of  polyacetylenes  decreases  while  that  of  the 
polycycllc  aromatlcs  Increases.   But  the  growth  of 
polyacetylenes  was  considered  too  slow  by  Calcote  (1981),  to 


account  for  the  rapid  formation  of  soot  particles. 
Furthermore,  as  pointed  out  by  Bonne  et  al .  (1965),  Thomas 
(1965)  and  Cullls  (1976),  the  rearrangement  of  a 
polyacetylene  to  an  aromatic  graphlte-1 Ike  structure  would 
also  be  a  slow  process. 

Another  hypothesis  proposed  by  Classman  (1980) 
emphasized  the  Importance  of  strongly  conjugated  molecules 
such  as  butadiene  In  the  formation  of  ring  structures  with 
side  chains  to  which  butadiene-type  molecules  keep  adding. 
This  mechanism  did  not  receive  much  support  on  the  basis  of 
concentration  profiles  In  which  the  butadiene  concentration 
is  observed  to  fall  very  rapidly  even  before  the  Initiation 
of  soot  formation  (Calcote , 1981 ) . 

Some  models  considered  Cg.  CgH  and  C3  radicals  as 
initial  nuclei  for  soot  formation  which  undergo  condensation 
reactions  with  each  other  or  with  acetylene  (Jensen,  1974; 
Smith,  1940;  Carbannes,  1956).   The  main  objection  to  these 
mechanisms  based  on  Cg  and  C3  condensation  reactions  was  the 
presence  of  a  very  strong  C2  radiation  in  rich  cyanogen- 
oxygen  flames  which  do  not  form  soot. 

Finally,  direct  condensation  of  aromatic  rings  was 
proposed  as  being  the  most  favored  thermodynamic  path  to 
soot  formation  (Bonne  et  al . ,  1965;  Graham  et  al . ,  1975; 
Stein,  1978).   However,  by  a  detailed  analysis  of  the 
concentration  and  the  flux  profiles  of  the  species,  Bittner 
and  Howard  (1981a,  1981b)  were  able  to  show  that  soot  is 


produced  via  nonaromatlc  hydrocarbon  intermediates  rather 
than  via  Intact  aromatic  rings.   But  after  a  consideration 
of  many  species,  only  a  few  intermediates  were  found  to 
favor  the  rearrangement  to  a  ring  structure  over  other 
decomposition  channels. 

To  conclude  the  discussion  above,  it  is  clear  that  all 
the  proposed  neutral  species  mechanisms  for  soot  formation 
have  serious  problems  with  the  nucleation  step,  and  thus 
none  can  be  considered  completely  satisfactory  in 
understanding  the  soot  formation  process.   Particularly  it 
is  difficult  in  any  of  the  proposed  neutral  schemes  to 
produce  small  cyclic  and  polycyclic  species  with  a 
sufficiently  high  rate  to  model  the  quite  rapid  soot 
formation  seen  in  flames.   Recently,  however,  a 
neutral/radical  soot  formation  model  has  been  proposed  which 
contains  many  steps  and  incorporates  quite  recently  obtained 
shock  tube  kinetic  data  (Frenklach  et  al . ,  1985).   This 
model  is  capable  of  yielding  small  ring  compounds  at  a 
higher  rate  than  previously  suggested  radical  mechanisms, 
although  there  is  some  question  as  to  the  correct  rate 
constant  for  one  or  two  of  the  crucial  steps  in  formation  of 
the  first  ring. 


Ionic  Mechanisms.   A  number  of  experimental 
observations  in  flames  have  led  to  a  consideration  of  the 
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iaportance  of  ions  In  soot  formation.   Such  observations  in 
flames  include  (Olson  and  Calcote,  1981a;  Calcote,  1981): 

1.  Presence  of  a  unit  charge  on  soot  particles; 

2.  A  positive  correlation  between  the  growth  of  large  Ions 
and  the  appearance  of  soot; 

3.  The  increase  both  in  the  concentration  of  large  ions  and 
the  onset  of  soot  formation  by  the  addition  of  electrophilic 
molecules ; 

4.  Identical  location  for  the  peak  concentrations  of 
molecular  ions  and  soot  precursors; 

5.  Production  of  larger  particles  by  an  increase  in  the 
residence  time  of  the  positive  ions;  and 

6.  High  ion  concentrations  in  hydrocarbon  pyrolysis. 

As  a  consequence  of  the  above  observations,  several 
ionic  mechanisms  have  been  suggested  (for  references,  see 
Calcote.  1981).   The  postulated  general  growth  of  soot 
particles  from  primary  ions  is  shown  in  Figure  1,1.   Calcote 
(1962)  proposed  the  chemlionlzatlon  reactions  (1.1)  and 
(1.2)  for  the  formation  of  primary  ions: 


CH  +  0 >  CHO"^  +  e" 


CH   +  CgHg 


->  CgHg"^  +  e 


(1.1) 
(1.2) 


The  C3H3   ion  is  the  dominant  Ion  observed  in  fuel-rich 
flames  (Knewstubb  and  Sugden,  1959;  Goodlngs  et  al . ,  1979; 
Olson  and  Calcote,  1981b:  Michaud  et  al . ,  1981)  and  its   1 
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Figure  1-1.   Growth  of  Soot  Particles  from  Primary  Molecular 
Species . 
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concentration  falls  very  rapidly  at  the  critical  equivalence 
ratio  of  soot  formation  where  the  concentration  of  the 
larger  Ions  starts  to  Increase  (see  Figure  1.2) 
(Calcote,  1981).   Unfortunately,  the  mechanism  of  CgHg"^ 
formation  Is  still  not  clear.   An  alternative  mechanism 
(Calcote,  1972)  for  CgHg^  formation  is  shown  by  reactions 
(1.3)  and  (1.4): 

(1.3) 


CHO"^  +  CgHgO >  CgHgO"^  +  CO 

CoH-jO"^  +  CoHo >  CoHo"^  +  CHoO 


'2"3' 


'2"2 


3"3 


(1.4) 


Reaction  (1.5)  is  postulated  (Calcote,  1972)  to  account  for 
the  dominant  HgO"*"  observed  in  near  stoichiometric  and  lean 
flames  : 


CHO"^  +  HgO >  HgO"^  +  CO 


(1.5) 


The  CoHg"^  ion  was  also  considered  as  a  possible  nucleating 
center  for  soot  formation  in  flames  and  its  formation  is 
shown  by  reaction  (1.6)  as  proposed  by  Vinckler  et  al, 
(1977)  . 


CHO"^  +  CgHg >  CgHg"^  +  CO 


(1.6) 


The  primary  ions  mentioned  above  are  proposed  to  react 
with  CjjHy  molecules  by  fast  ion/molecule  reactions  (Calcote, 
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Figure  1.2.   Effect  of  Equivalence  Ratio  on  Peak  Ion 
Currents  for  2.0  kPa  Acetylene-Oxygen  Flames.   The  shaded 
area  Indicates  the  minimum  equivalence  ratio  for  soot 
formation . 


Reprinted  by  permission  of  Elsevier  Science  Publishing  Co., 
Inc.  from  the  article  by  H.  Calcote  (1981).  Copyright  1981 
by  the  Combustion  Institute. 
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1981;  Olson  and  Calcote,  1981a).   The  typical  set  of 
reactions  for  C3H3^  are  shown  below: 


C3H3  +  C2H2 

C3H3  +  CgHg 

C3H3  +  C4H2  ■ 

C3H3  +  C4H2  • 


--->  C5H3- 
— >  C5H5' 


■-->  C5H3   *  H 

>  C5H3   +  C2H2 


(1.7) 

(1.8) 

(1.9) 

(1.10) 

(1.11) 


The  next  step  is  the  further  addition  of  neutral 
building  blocks,  acetylene  and  diacetylene  to  the  product 
ions  producing  larger  ions  with  increased  C:H  ratio. 
Following  are  a  few  reactions  representative  of  this  step 


^5^3   "^  ^2^2 


^5^3   ^  ^4^2 


■->  C^Hg- 


>  CgHg- 


CyHg     +   C4H2  >   CgHg     +   C2H2 


(1.12) 
(1.13) 

(1.14) 


This  series  of  reactions  is  suggested  to  lead  to  the 
formation  of  polycyclic  structures  such  as  C,oHg^  and 


^19"ll 


L 


dop 
oia 


^13^9 


C  1  oH 


19"  11 


11 

The  problem  of  rearrangement  of  linear  structures  to 
polycyclic  structures  encountered  in  the  neutral  mechanism 
is  overcome  in  the  ionic  mechanism  because  of  the  general 
observation  that  gaseous  ions  usually  rearrange  upon 
formation  to  their  most  stable  structure  (Calcote,  1981). 

Mlchaud  et  al .  (1981)  have  made  an  alternate  suggestion 
that  direct  reactions  of  CgHg"^  with  aromatic  neutrals  such 
as  benzene,  toluene,  naphthalene,  methylnaphthalenes ,  and 
indene  may  be  more  important  in  forming  polycyclic  ions  than 
sequential  reactions  involving  acetylene  and  diacetylene. 
In  fact,  recent  ion  cyclotron  resonance  (ICR)  studies  of 
reactions  of  CgHg"*"  with  aromatic  neutrals  showed  (Baykut  et 
al..  1986)  that  these  reactions  are  fast  enough  to  be 
considered  as  possible  bypass  channels  in  ionic  soot 
formation  pathways. 

The  main  objection  to  the  ionic  mechanism  is  the  lower 
concentration  of  ions  in  flames  compared  to  that  of  neutral 
species.   Figure  1.3  (Calcote,  1981)  shows  the  concentration 
profiles  of  primary  neutrals  and  ions  observed  in  flames  as 
a  function  of  distance  above  the  burner.   One  of  the  most 
important   features  of  the  figure  is  that  the  soot 
concentration  is  lower  than  that  of  the  large  positive  ions 
which  implies  that  there  are  probably  enough  of  these  large 
Ions  to  produce  the  observed  levels  of  soot.   Several  orders 
of  magnitude  higher  concentrations  of  neutral  species 
compared  to  those  of  ions  in  flames  can  be  explained  as  a 
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requirement  for  the  ion/molecule  reactions  to  proceed  fast 
enough  to  account  for  the  rapid  formation  of  soot. 

In  this  study,  the  reactions  of  primary  ions  with  flame 
neutrals  have  been  investigated  to  identify  different 
reactant  ion  structures,  reaction  products  and  mechanisms. 
One  of  the  goals  of  this  study  was  to  determine  how  fast 
these  ion/molecule  reactions  proceed  under  laboratory 
conditions,  which  hopefully  leads  to  a  better  understanding 
of  the  complex  phenomena  occurring  under  flame  conditions. 
Fourier  Transform  Ion  Cyclotron  Resonance  (FTICR)  mass 
spectrometry,  used  in  this  study,  is  well  suited  to 
Investigate  ion/molecule  reactions  due  to  its  inherent 
qualities  which  differ  from  other  mass  spectrometers.   These 
qualities  Include  very  clean  isolation  of  the  mass  of 
Interest  from  a  complicated  reaction  mixture  and  accurate 
rate  constant  measurements  along  with  reaction  mechanism  and 
branching  ratio  studies.   The  basics  of  this  technique  are 
discussed  in  Chapter  2  along  with  an  explanation  of  the 
determination  of  some  experimental  pressure  measurement 
parameters  such  as  Baratron  factors  and  the  system  factor 
inherent  to  the  instrument.   The  use  of  FTICR  mass 
spectrometry  for  studying  ion/molecule  reactions  is 
discussed  in  Chapter  3  along  with  some  theoretical  models 
developed  for  ion-molecule  reactions.   The  main  body  of  the 
experimental  work,  reported  in  Chapters  4-7,  deals  with 
proposed  nucleatlon  steps  of  the  ionic  mechanism  of  soot 
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formation.   An  overall  conclusion  of  the  experimental 
results  is  Included  in  Chapter  8. 


CHAPTER  2 
INSTRUMENTATION 


Fourier  Transform  Ion  Cyclotron  Resonance  (FTICR)  Mass 
Spectrometry  was  introduced  in  the  mid-1970s  by  Comisarow 
and  Marshall  (1974a.  1974b.  1974c,  1975).   It  involves  the 
application  of  Fourier  multiplex  concepts   (Griffiths,  1978; 
Marshall,  1982)  to  ICR  mass  spectrometry.   High  speed,  high 
sensitivity,  computerization,  ultra-high  mass  resolution, 
and  wide  mass  range  are  some  of  the  advantages  of  Fourier 
Transform  over  continuous-wave  spectrometry  (Marshall, 
1985).   Basic  principles  and  the  vast  applications  of  this 
relative  newcomer  to  mass  spectrometry  are  summarized  in 
several  recent  review  articles  (Johlman  et  al . ,  1983;  Gross 
and  Rempel,  1984;  Comisarow,  1985;  Marshall,  1985;  Baykut 
and  Eyler,  1986;  Laude  et  al.,1986).   In  this  chapter,  the 
basic  theory  of  operation  of  ICR  which  led  to  the  concepts 
and  development  of  the  FTICR  technique  will  first  be 
discussed.   Following  that  is  a  description  of  the 
instrument  along  with  typical  experimental  parameters  and 
the  general  sequence  of  operation. 

Theory .   An  ion  of  charge,  q,  moving  at  velocity,  v 
when  put  in  a  uniform  magnetic  field,  B,  is  subject  to  the 
Lorentz  force, 
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ma  =  q  ( V  X  B ) 


(2.1) 


which  acts  perpendicular  to  the  direction  of  Ion  motion. 
Under  the  Influence  of  this  force,  the  Ion  follows  a  helical 
path  (Lawrence  and  Livingston,  1932)  which,  when  projected 
into  a  plane  perpendicular  to  the  magnetic  field  gives  a 
circle  with  a  radius,  r  proportional  to  the  velocity  of  the 
Ion,  as  shown  in  Figure  2.1.   Substitution   of  the 
centripetal  acceleration  in  terms  of  angular  frequency,  ©, 
into  equation  (2.1)  gives 


and 


IFI  =  mlal  =  mv'/r  =  mrco*^  =  qorB 
CO  =  qB/m  (mks  units) 


(2.2) 
(2.3) 


which  relates  the  ion's  characteristic  cyclotron  frequency 
to  its  q/m  (charge/mass).   Equation  (2.3)  is  called  the 
cyclotron  equation  and  provides  the  basic  principle  of  mass 
measurement  in  terms  of  cyclotron  frequency  in  the  ICR  mass 
spectrometer.  For  a  magnetic  field  of  3.0  Tesla,  singly 
charged  ions  with  masses  In  the  range  of  18-5000  a.m.u.  have 
cyclotron  frequencies  in  the  radiof requency  range  (10  kHz- 
2.6  MHz),  within  which  frequency  can  be  measured  with  high 
precision . 

Operation .   A  schematic  representation  of  the 
commercial  Nicolet  FTMS-1000  mass  spectrometer  with  a 
superconducting  solenoid  magnet  (3  Tesla)  Is  shown  in  Figure 
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Figure  2.1.   Motion  of  an  Ion  with  Initial  Velocity  v  In  a 
Magnetic  Field  B. 
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2.2.   All  data  acquisition  and  processing  and  the  selection 
of  various  experimental  parameters  (except  emission  current) 
are  under  computer  control.   Figure  2.3  shows  a   simplified 
block  diagram  of  an  FTICR  mass  spectrometer.  The  simplicity 
of  operation  results  from  the  ability  to  form,  excite,  react 
and  detect  ions  in  the  same  1-inch  cubic  stainless  steel 
cell  shown  in  Figure  2.4.   This  analyzer  cell  is  located  in 
the  center  of  a  vacuum  chamber  which  passes  through  the  bore 
of  the  magnet.   The  background  pressure  of  10"^  torr  is 
achieved  by  the  use  of  a  four-Inch  diameter  oil  diffusion 
pump  and  by  baking  out  the  system  at  250*C  for  several  hours 
each  night. 

An  inlet  system  which  is  evacuated  with  a  three-inch 
diameter  diffusion  pump  permits  the  introduction  of  gas  and 
liquid  samples  into  the  high  vacuum  through  precision  leak 
valves.  The  liquid  samples  are  used  after  multiple  freeze- 
pump-thaw  cycles  to  remove  non-condenslble  gases.   A  solids 
insertion  probe  is  used  for  the  introduction  of  solid 
samples  into  the  vacuum  chamber  and  can  be  heated  up  to 
250°C  to  promote  sublimation  of  the  solid. 

Ion  formation  is  typically  achieved  in  5  ms  by 
collisions  of  electrons  accelerated  to  10-70  eV  with 
neutrals  leaked  into  the  vacuum  chamber.   The  number  of  ions 
is  controlled  by  monitoring  the  emission  current,  measured 
by  means  of  an  electron  collector  located  behind  the  trap 
plate  opposite  to  the  filament.   When  ions  are  formed,  they 
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Figure  2.3.   Simplified  Block  Diagram  of  an  FTICR  Mass 
Spectrometer . 
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Figure  2.4.   Cubic  FTICR  Cell. 
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start  cyclotron  motion  at  their  characteristic  frequency 
with  a  random  phase.   The  trapping  plates  shown  In  Figure 
2.4  constrain  the  ion  motion  in  the  direction  parallel  to 
magnetic  field  lines.  Positive  and  negative  ions  are  trapped 
by  the  application  of,  typically,  positive  and  negative  1 
volt  potentials  to  the  trapping  plates,  respectively. 

To  detect  the  ions,  a  fast  sweep  of  oscillating  voltage 
in  the  radiof requency  range  is  applied  to  the  transmitter 
plates  shown  in  Figure  2.4.   Ions  absorb  energy  at  their 
cyclotron  frequency  and  are  driven  into  coherent  motion  with 
a  larger  radius,  thus  inducing  an  image  current  oscillating 
at  the   frequency  of  their  cyclotron  motion  on  the  receiver 
plates  (see  Figure  2.5).   The  image  current  is  amplified, 
digitized  and  stored  in  the  computer's  memory  as  a  time- 
domain  signal  which  contains  the  superimposed  frequency 
Information  of  all  the  ions  with  different  masses  in  the 
analyzer  cell.   A  fast  Fourier  Transform  (Cooley  and  Tukey, 
1965)  is  applied  to  the  time-domain  data  to  produce  the 
frequency-domain  spectrum  which  can  be  plotted  In  terms  of 
the  ion  masses.   Figure  2.6a  shows  a   time-domain  spectrum 
of  two  superimposed  sine  waves  of  frequencies  <Uj  and  ©2  . 
The  decay  of  the  signal  results  from  the  dephasing  of  ions 
due  to  collisions  with  the  neutral  molecules.   Fourier 
transformation  of  this  spectrum  gives  the  mass-domain 
spectrum  shown  in  Figure  2.6b.   In  this  procedure,  the 
lowest  mass  detected  is  limited  by  the  sampling  frequency  of 
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(a) 


(b) 


(c) 


Figure  2.5.   Excitation  and  Detection  of  Ion  Motion. 

(a)  Ion  motion  In  magnetic  field  only,  (b)  motion  during 

radio-frequency  pulse  applied  to  the  transmitter  plates, 

(c)  larger  radius  of  Ion  motion  after  excitation,  and 

(d)  generation  of  an  Image  current  In  the  receiver  plates  by 
a  rotating  Ion  "clump". 
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Figure  2.6.   Ion  Detection  by  Fourier  Transform.   (a)  Time- 
domain  spectrum  (Abscissa  is  time;  ordinate  is  ion  signal), 
(b)  Frequency-Domain  Spectrum  after  Fast  Fourier  Transform 
of  the  time-domain  data  (Abscissa  is  frequency;  ordinate  Is 
1  on  signal ) . 
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the  analog-to-dlgltal  converter,  which  should  be  twice  the 
signal  frequency  according  to  the  Nyquist  criterion.   Since 
the  highest  cyclotron  frequency  detected  is  one-half  the  ADC 
frequency,  for  a  magnetic  field  of  3  Tesla  used  in  this 
work,  the  singly  charged  ion  mass  corresponding  to  this 
frequency  is  18  a.ra.u.,  which  sets  the  low  mass  limit. 

A  typical  experimental  sequence  is  shown  in  Figure  2.7. 
Ions  can  be  manipulated  in  various  ways  between  the 
formation  and  the  excitation  pulses.   Ion  ejection  is  used 
to  clean  the  cell  of  all  ions  except  one  of  interest.   Ions 
are  ejected  by  the  same  principle  as  excitation.   The  only 
difference  is  that  a  higher  amplitute  of  excitation  voltage 
is  used.   In  this  way  ions  achieve  large  enough  radii  to 
strike  the  cell  walls  and  be  neutralized.   Application  of 
either  a  single  ejection  pulse  at  a  single  cyclotron 
frequency  or  a  swept  pulse  covering  a  range  of  frequencies 
is  determined  according  to  the  user's  purpose. 

After  selection  of  the  ion  of  Interest,  different  kinds 
of  ion  chemistry  can  be  performed.  Among  these,  ion/molecule 
reactions,  col 1 isional ly  activated  dissociation,  and 
photodissociatlon  are  the  most  commonly  used  processes  in 
characterization  of  ion  structures.   Since  the  number  of 
stages  of  ion  ejection  and  manipulation  is  not  limited  by 
hardware  or  software,  (MS)"  (multiple  stages  of  mass 
selection  and  manipulation)  experiments  can  be  performed  by 
repeating  the  same  processes  for  the  various  generations  of 
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offspring  ions.   Laser  desorptlon  of  solid  samples  and  laser 
ionization  are  other  techniques  used  in  conjuction  with 
FTICR. 


p.'-' 


M: 


Experimental  Parameters.   In  this  study,  the  FTICR 
technique  was  used  mainly  to  determine  ion/molecule  reaction 
rates  and  mechanisms.   Ion  formation  was  accomplished  by 
dissociative  ionization  using  various  charge  transfer  agents 
instead  of  electron  impact  ionization,  in  order  to  form  the 
ions  with  well-defined  internal  energies.   Reaction  pathways 
were  delineated  by  using  the  ejection  capabilities  of  FTICR 
which  make  it  possible  to  eject  one  ion  from  a  complicated 
reaction  mixture  to  determine  its  contribution  to  the  mass 
spectrum  of  all  the  other  ions. 

Reaction  rate  coefficients  were  determined  by 
monitoring  the  intensity  of  the  reactant  ion  as  a  function 
of  time  (typically  for  at  least  2s)  after  ejection  of  all 
other  ions  from  the  analyzer  cell.   Neutral  gas  pressures 
were  measured  with  an  ionizing  gauge.   Ionization  gauge 
readings  were  then  corrected  by  constructing  calibration 
curves  of  ionization  gauge  vs.  capacitance  manometer  (MKS- 
Baratron)  readings  in  the  1  x  10~®  to  1  x  10"^  Torr  range. 
In  order  to  correct  for  the  fact  that  the  ionization  gauge 
and  capacitance  manometer  were  located  at  different  points 
on  the  vacuum  system,  somewhat  removed  from  the  FTICR 
analyzer  cell,  a  correction  factor  was  required.   This  was 


28 
obtained  by  determining  the  rate  coefficient  of  a  well- 
studied  reaction  (CgH^"^  +  C2H4 >  products),  where  the 

ionization  gauge  pressure  readings  were  corrected  by  using 
the  capacitance  manometer.   This  experimentally  determined 
rate  coefficient  was  then  compared  with  the  average  of 
published  values,  kj^^  =  (1.0  +  0.3)  x  10"^  cm^/s  (Herod  and 
Harrison,  1970;  Gross  and  Norbeck,  1971;  Sleek  and  Ausloos, 
1972;  Warneck,  1972;  Le  Breton  et  al . ,  1975)  and  the  ratio 
of  the  measured  value  to  the  published,  which  was  0.3  +  0.1, 
was  used  as  a  correction  factor.   This  factor  was  used  in 
calculating  the  absolute  rate  coefficients  reported  in  this 
work.   The  large  uncertalnities  (95*  confidence  limits) 
reported  for  the  correction  factor  and  rate  coefficients 
calculated  using  it  are  primarily  due  to  the  wide  range  of 
reported  values  for  the  CgH^'*'  +  ^2^\    reaction  used  to 
determine  the  correction  factor.   The  even  more  widely 
studied  "standard"  calibration  reaction  of  CH^"^  with  CH^ 
could  not  be  employed  because  of  the  lowest  accessible  mass 
limit  in  the  FTICR  Instrument  mentioned  earlier.   All 
calculations  of  rate  coefficients  and  95%  confidence  limits 
were  performed  with  a  menu-driven  Fortran  computer  program 
(given  in  Appendix  I). 
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CHAPTER  3 
ION/MOLECULE  REACTIONS 


In  this  chapter  the  various  theoretical  models  which 
have  been  developed  for  Ion/molecule  reactions  will  first  be 
discussed.   Then  follows  a  description  of  practical  aspects 
of  studying  such  reactions  by  the  use  of  FTICR  mass 
spectrometry. 

Theory  of  Ion/Molecule  Reactions 
A  number  of  both  classical  and  statistical  Ion/molecule 
collision  theories  have  been  introduced  during  the  last  two 
decades  to  provide  an  adequate  model  explaining  the 
experimental  observations  of  ion/molecule  reactions. 
Classical  treatments  mainly  Include  the  pure  polarization 
(ion-induced  dlpole)  theory  first  developed  by  Langevln 
(1905)  and  the  various  lon-dipole  theories. 

A.  Pure  Polarization  Theory.   Detailed  discussions  of 
this  model  can  be  found  in  a  number  of  articles  (Gloumousis 
and  Stevenson,  1958;  McDaniel,  1964;  Futrell  and  Tiernan, 

'  ^       1968;  Henglein,  1970).   Langevin's  model  assumes  that  the 
neutral  has  no  permanent  dlpole  moment,  and  that  both  the 

^-  >  .  .   Ion  and  the  neutral  molecule  are  point  particles  with  no 
Internal  energy.   The  classical  charge-induced  dlpole 

^i'^r:  potential  at  an  ion-molecule  separation  r  is 

t'"'  ■'•  2  9 
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V(r)  =  -  aq2/2r* 


(3.1) 


where  q  Is  the  charge  on  the  Ion  and  a    is  the  polarlzabl 1 ity 
of  the  neutral.   Since  the  energy  of  relative  rotation, 
Ej.Q^(r)  of  the  particles  is  associated  with  an  outwardly 
directed  centrifugal  force,  the  effective  potential  of  the 
ion/molecule  system  can  be  given  by 


V,,.(r)  =  V(r)  +  E^„.(r) 


rot 


or 


Vgff(r)  =  -(q^a/2r^)    +  (L^/2fiT^) 


(3.2) 


(3.3) 


where  L  is  the  classical  orbiting  angular  momentun  of  the 
two  particles  and  fi    is  the  reduced  mass.   The  total  relative 
energy  of  the  system  is  a  sum  of  the  translational  energy 
and  the  effective  potential  energy. 


^r  "  ^trans 


(r)  +  Veff(r) 


(3.4) 


Figure  3.1  shows  a  plot  of  Vg^^Cr)  versus  r  at  constant  E^. 
for  three  different  values  of  the  Impact  parameter,  b.   When 
b=0,  since  there  is  no  centrifugal  potential  energy,  the 
effective  potential  is  attractive  at  all  ion/molecule 
separations  leading   to  a  collision  between  two  particles. 
For  b  >  0,  there  is  a  critical  value  of  the  impact 
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Figure  3.1.   Plot  of  Vg^^  vs  r  from  Equation  (3.3)  for  Ng 
Colliding  with  Ng . 
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paramater,  b^  at  which  the  particles  orbit  around  the 
scattering  center  with  a  constant  radius,  r^ .   At  this 
lon/Bolecule  separation,  there  Is  no  contribution  from  the 
attractive  potential  and  Vg^j.  Is  equal  to  centrifugal  energy 
which  creates  a  "centrifugal  barrier"  to  a  capture 
collision.   For  all  b  <  b^,  a  capture  collision  occurs, 
whereas  It  Is  precluded  by  the  centrifugal  barrier  for  all 
b  >  bjj. 

Capture  cross  section  Is  defined  as  the  area  of  a 
circle  with  radius  b^  perpendicular  to  the  line  of  collision 
at  infinite  ion/molecule  separation  and  a  capture  collision 
occurs  for  all  the  ions  that  pass  through  the  circle  when 
approaching  the  neutral  molecule.   Capture  cross  section  can 
be  derived  for  a  given  relative  velocity  in  terms  of  the 
charge  of  ion,  q,  polar Izabllity  of  the  neutral,  a,  and  the 
reduced  mass,  [i    as  shown  in  equation  (3.5). 


a_(v)  -  2iiq{a/fi)^^^/v 


(3.5) 


Thus  the  collision  rate  constant  is  given  by 


kg  =  va-  =  2nq{a/(i) 


1/2 


(3.6) 


The  rate  coefficients  obtained  from  this  expression  are 
generally  good  for  some  simple  low  energy  ion/molecule 
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reactions  but  underestimate  the  rate  constants  of  most 
ion/polar  molecule  collisions. 

B.  lon/Dlpole  Theory.   Several  different  models  were 
proposed  to  describe  the  theory  of  ion/dipole  collisions. 
Among  these  are  "The  Locked  Dipole  Approximation"  (Moran  and 
Hamill,  1963).  "The  Frozen  Rotor  Approximation"  (Dugan  and 
Magee,  1966),  "lon/Dipole  Trajectory  Calculations"  (Dugan 
and  Magee,  1967),  "The  Average  Dipole  Orientation  (ADO) 
Theory"  (Bowers  and  Laudenslager ,  1972),  "Barker-Ridge  (BR) 
Model"  (Barker  and  Ridge,  1976)  and  "The  Average  Dipole 
Orientation  Theory  with  Conservation  of  Angular  Momentum 
(AADO)"  (Su  et  al . ,  1978)  . 

Both  the  "Locked  Dipole"  and  "The  Frozen  Rotor" 
approximations  assume  that  the  polar  molecule  is  fixed  at 
angle  0    {0    =    0    for  the  "Locked  Dipole"  model)  with  respect 
to  the  line  of  centers  of  the  collision.   The  resulting  rate 
constants  provide  upper  limits  to  the  ion/dipole  capture 
collisions  and  overestimate  the  dipole  effect  (Su  and 
Bowers ,  1973 ) . 

lon/Dipole  Trajectory  Calculations  provide  numerical 
solutions  for  the  equations  of  motion  for  the  collision  of 
an  ion  with  a  rotating  polar  molecule.   As  a  result  of  these 
calculations,  the  numerical  capture  cross  section  was  found 
to  depend  on  relative  translat ional  energy.   When  compared 
with  experimental  results,  this  numerical  approach  seems  to 
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be  useful  In  the  Investigation  of  nonreactive  ion/dipole 
phenomena  such  as  energy  transfer,  but  does  not  give  a  good 
estimation  of  capture  collision  rates. 

ADO  theory  assumes  that  there  exists  an  overall  average 
orientation  of  the  dipole  with  respect  to  the  ion  and  the 
extent  of  locking  of  the  dipole  Increases  as  the  dipole 
moment  of  the  molecule  increases.   Su  and  Bowers  (1973)  have 
parametrized  the  ADO  theory  to  calculate  the  ADO  capture 
rate  constants.   The  results  indicate  that  the  ADO  theory 
gives  the  best  prediction  of  rate  constants  in  most  cases 
with  respect  to  the  other  theories. 

BR  theory  uses  a  simple  statistical  approach  to  study 
the  ion/dipole  collisions  and  is  based  on  the  Langevin 
model.   BR  results  overestimate  the  dipole  effect  by  almost 
a  factor  of  2  for  charge  transfer  and  proton  transfer  as 
well  as  momentum  transfer  reactions. 

Su  et  al.  (1978)  modified  the  ADO  theory  by  the 
application  of  the  conservation  of  angular  momentum  with  the 
assumption  that  there  is  no  net  angular  momentum  transfer 
between  the  rotating  molecule  and  the  ion/molecule  orbital 
motion.   This  modified  version  of  ADO  theory  is  termed  the 
AADO  theory.   The  capture  rate  constants  obtained  from  AADO 
model  are  larger  than  the  ADO  predictions,  providing  better 
agreement  between  theory  and  experiment. 
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C.  lon/Quadrupole  Theory.   This  theory  makes 
theoretical  predictions  of  lon/quadrupole  Interactions  for 
the  molecules  having  D„  symmetries.   Similar  to  ADO  theory, 
the  average  quadrupole  orientation  theory,  AQO,  was 
developed  by  Su  and  Bowers  (1975).   This  model  predicts 
significant  quadrupole  effects  for  molecules  with  high 
quadrupole  moments  and  Improves  the  agreement  between  theory 
and  experiment. 

Statistical  Treatments  of  Ion/molecule  reactions 
Include  the  statistical  theory  and  the  orbiting  transition 
state  theory.   Statistical  Theory  is  divided  Into  three 
categories  according  to  the  type  of  reaction  dynamics  as 
shown  In  Figure  3.2.  - 

A.  Transition  State  Theory.   This  theory  Is  developed 
for  direct  blmolecular  reactions  (Pelzer  and  Wlgner,  1932; 
Evans  and  Polanyi,  1935;  Eyring,  1935).   This  method  was 
later  refined  by  Pechukas  and  McLafferty  (1972)  and  by 
Miller  (1974).   The  mlcrocanonical  rate  constant  is  given  by 


k(E)  =  W*(E-Eo)/hp(E) 


(3.7) 


where  w'''(E-Eq)  is  the  total  sum  of  states  (vibrational  and 
rotational)  of  the  transition  state  complex  with  nonflxed 
energy  less  than  or  equal  to  E-Eq.  where  E  is  the  total 
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Figure  3.2.   Schematic  Potential  Surfaces  for  Various 
Ion/Molecule  Reactions,  (a)  Direct  blmolecular ,  (b) 
Unlmolecular ,  and  (c)  Complex  formation  reactions. 
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system  energy,  and  p(E)  Is  the  total  density  of  states  of 
the  reactants.   For  a  given  potential  surface,  transition 
state  theory  provides  a  rigorous  upper  limit  to  the  true 
rate  constant. 

B.  Unlmolecular  Theory.   This  theory  involves  two 
approaches,  RRKM  and  QET  ( quasi-equilibrlum  theory)  which 
arrive  at  similar  results  for  the  unlmolecular  reaction  of 
an  isolated  energized  molecule.   These  methods  were 
developed  Independently  by  Marcus  and  Rice  (1951)  using  the 
assumptions  of  transition  state  theory  and  by  Rosenstock  et 
al.  (1952)  using  the  quasi-equl llbrlum  hypothesis. 
Experimental  evidence  suggests  that  most  unlmolecular 
systems  appear  to  behave  according  to  the  quasi-equilibrlum 
hypothesis  which  assumes  a  random  distribution  of  excess 
energy  Independent  of  the  pattern  of  initial  energization 
prior  to  unlmolecular  decomposition. 

C.  Complex  Formation  Theory.   This  theory  Involves  the 
statistical  treatment  of  reactions  which  form  a  long-lived 
complex  assuming  quasi-equilibrlum  for  the  intermediate 
complex  (Keck,  1958;  Light,  1967).   Thus,  the  rate  constant 
for  decomposition  to  each  available  channel  is  assumed  to  be 
proportional  to  the  flux  through  the  transition  state  for 
that  channel . 
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The  Orbiting  Transition  State  Theory  Is  developed  to 
overcome  the  problem  of  determining  the  energy  distribution 
of  reactants  that  leads  to  the  formation  of  a  given 
transition  state  complex.   In  this  method  the  transition 
state  Is  located  at  the  maximum  of  the  long  range  effective 
potential  where  the  reactants  retain  their  original 
identity.   The  orbiting  rate  constants  for  both  blmolecular 
reactions  (Miller,  1976;  Light.  1967),  and  unlmolecular 
reactions  (Chesnavich  and  Bowers.  1977)  can  be  obtained  as  a 
function  of  translational .  rotational,  and  vibrational 
energy  of  the  reacting  system.   Comparison  with  experiment 
in  general  reveals  that  orbiting  unlmolecular  rate  constants 
are  unreallstically  high  while  the  orbiting  blmolecular  rate 
constants  are  similar  to  experimental  results  at  thermal 
energies.   This  model  also  predicts  the  probability  of  a 
given  translational  energy  release  (Chesnavich  and  Bowers, 
1976)  which  in  most  cases  is  in  good  agreement  with 
experimental  results  for  reactions  with  no  reverse 
activation  energy. 

Experimental  Study  of  Ion/Molecule  Reactions  Using  FTICR 


Rate  Coefficients.   ICR  techniques  cannot  be  used  for 
experimental  tests  of  statistical  theories  since  the 
application  of  statistical  models  requires  a  detailed 
knowledge  of  energy  states  of  the  system,  which  cannot  be 
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obtained  from  ICR.   On  the  other  hand,  theoretical 
blmolecular  rate  constants  using  Langevln's  model  are 
usually  found  to  agree  with  the  experimental  results 
obtained  for  low  energy  Ion/molecule  reactions  using  ICR 
technique.   Reactant  ions  with  energies  close  to  thermal  can 
be  prepared  by  charge  transfer  ionization  of  a  precursor 
molecule  followed  by  colllsional  relaxation  of  the  ions  by 
inert  gas  molecules.  Thus,  this  method  provides  the 
determination  of  thermal  blmolecular  rate  constants  which 
can  be  used  to  test  the  basic  assumptions  of  the  Pure 
Polarization  method. 

The  first  step  in  the  experimental  determination  of 
rate  coefficients  by  FTICR  involves  the  ejection  of  all  ions 
except  the  one  of  interest  from  the  analyzer  cell  following 
charge  transfer  ionization.   Next  the  decay  of  the  ion 
signal  is  monitored  as  a  function  of  time  long  enough  to 
permit  the  ion  of  interest  to  react  away  completely. 


Method  of  Data  Analysis.   The  ion/molecule  reactions 
studied  in  FTICR  most  often  involve  blmolecular  collisions 
and  follow  second  order  kinetics.   The  rate  expression  for  a 
simple  second  order  reaction  is  given  by 


-d[I]/dt  =  k2[I][R]  (3.8) 

where  I  is  the  ion  and  R  is  the  reactant  neutral.   Since  the 
ion  number  density  created  in  the  FTICR  instrument  is  at 
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least  three  orders  of  magnitude  lower  than  the  neutral 
nunber  density,  a   "pseudo-first  order"  rate  coefficient, 
kj   can  be  substituted  for  the  terra  k2[R]  in  equation  (3.8) 


-d[I]/dt  =  kj  [I] 


(3.9) 


After  integration,  the  rate  expression  giving  the  number 
density  of  the  reactant  at  any  given  time  is  obtained: 


[I](t)  =  [I](0)  exp(-kj  t) 


(3.10) 


where  [I](0)  is  the  initial  number  of  ions  present.   Thus, 

kj   can  be  determined  from  the  slope  of  the  simple  ln(ion 

signal)  vs  time  plot.   The  true  rate  coefficient,  k2  can 

then  be  obtained  from  kj   by  dividing  kj   by  the  number 

density  of  the  reactant  neutral.   The  number  density  of  the 

neutral  Is  determined  using  the  expression,  (n/V)  =  P/RT, 

where  V  and  T  are  the  volume  and  the  temperature  of  the 

analyzer  cell  respectively,  and  P  is  the  pressure  of  the 

reactant  neutral  measured  by  a  capacitance  manometer. 

In  an  experimental  situation,  there  is  often  a  second 

neutral  present  in  the  cell,  which  is  the  precursor  neutral 

of  the  ion  of  interest.   Thus,  the  decay  of  the  ion  signal 

Includes  Ion/molecule  reactions  with  both  the  precursor  and 

the  neutral  molecules  of  interest.   In  this  case  the  rate 

expression  becomes 
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-d[I]/dt  =  (kglR]  +  k3[P])[I] 


(3.11) 


where  kg  is  the  blmolecular  rate  coefficient  for  the 
reaction  of  the  ion  with  the  precursor  molecule,  P.   The 
Integrated  form  is 


[I](t)  =  [I](0)  exp  [-(kgER]  +  k3[P])t] 


(3.12) 


In  order  to  determine  the  rate  coefficient  of  reactions  of 
the  ion  with  the  neutral  reactant,  it  is  necessary  first  to 
monitor  the  reactions  of  the  ion  with  precursor  neutral  with 
no  reactant  neutral  present  in  the  analyzer  cell  and  then  to 
subtract  the  rate  coefficient  for  this  reaction  from  the 
total  rate  coefficient  observed  in  the  presence  of  both  the 
precursor  neutral  and  the  reactant  of  Interest.   If  the 
precursor  neutral  pressure  is  kept  constant  for  the  two 
cases,  only  the  pseudo-first  order  rate  constants  are  used 
in  the  subtraction  step. 

When  two  or  more  isomers  of  a  particular  ion  are  formed 
from  the  same  parent  neutral,  a  single  ion  signal  is 
observed  in  the  mass  spectrometer  whose  intensity  is  the  sum 
of  the  intensities  of  all  the  isomers.   For  a  simple  case  of 
only  two  ion  structures. 


[I](t)  =  [A](t)  +  [B](t) 


(3.13) 
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Thus,  the  integrated  form  of  the  rate  expression  Is  given  by 


[I](t)  =  [A](0)  exp[-(k2[R]  +  k3[P])t] 

+  [B](0)  exp[-(k2'[R]  +  k3' [P])t] 


(3.14) 


where  A  and  B  are  the  isomers  of  ion  I.   Consequently,  the 
decay  curve  of  the  ion  signal  is  composed  of  segments  with 
two  different  exponential  time  constants,  each  slope 
corresponding  to  the  reaction  rate  of  a  different  ion 
structure.   The  sum  of  the  rate  coefficients  (for  reactions 
with  precursor  and  reactant  neutrals)  of  the  slower  reacting 
structure  is  calculated  from  the  tail  of  the  double  decay 
curve  while  that  of  the  faster  reacting  structure  is 
determined  by  subtracting  the  amount  of  the  slow  reacting 
isomer  at  each  point  in  the  initial  part  of  the  curve  using 
Equation  (3.14).   Often,  one  of  the  structures  is  completely 
unreactive,  while  the  other  Is  reactive  as  shown  in  Figure 
3.3,  simplifying  the  data  analysis  for  structure 
identification. 


Reaction  Mechanisms.   Reaction  pathways  leading  to 
final  product  ions  from  a  ion/molecule  reactant  system  can 
be  studied  using  the  ejection  capabilities  of  FTICR.   The 
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method  Involves  the  ejection  of  each  first  generation 
product  ion  continuously  during  the  reaction  time  to  see  Its 
effect  on  the  Intensity  distribution  of  the  mass  spectrum  of 
the  final  product  ions.   This  method  also  can  be  used  to 
calculate  the  branching  ratios  for  the  product  Ions. 
Studying  reaction  pathways  and  branching  ratios  of  the  ions 
produced  in  different  ways  (either  from  different  precursors 
or  different  methods  of  Ionization)  provides  an  alternative 
means  to  rate  coefficient  determination-  for  the 
differentiation  of  structural  Isomers  of  the  ion  of 
Interest . 


CHAPTER  4 
REACTIONS  OF  C3H3'^  WITH  ACETYLENE  AND  DIACETYLENE 

IN  THE  GAS  PHASE 


Introduction 


The  CgHg"^  Ion  has  received  considerable  attention  In 
recent  years  as  a  possible  soot  precursor  because  It  is 
found  In  quite  high  abundance  In  fuel-rich  and  sooting 
flames  (Goodlngs  et  al ,  1979;  Olson  and  Calcote,  1981b; 
Mlchaud  et  al . ,  1981).   Although  substantial  uncertainty 
remains  as  to  CgHg"*"  formation  mechanisms  In  flames  (Calcote, 
1981),  the  Ion  Is  postulated  (Calcote,  1981;  Olson  and 
Calcote,  1981a)  to  react  with  neutrals  such  as  acetylene, 
dlacetylene,  and  CgH  In  rapid,  sequential  condensation  and 
condensation/elimination  reactions,  forming  successively 
larger  Ions,  which  can  rearrange  to  cyclic  species  during 
the  ion/molecule  reaction  chain. 

Two  Isomeric  structures  of  CgHg"*^  shown  in  Figure  4.1 
are  important  in  discussing  the  role  of  C3H3*  in 
ion/molecule  reaction  mechanisms.   The  first  is  the 
cyclopropenylium  ion,  recognized  as  the  most  stable  Isomer, 
which  has  a  theoretically  calculated  heat  of  formation  of 
253  kcal/mol  (Radom  et  al . ,  1976).  in  quite  good  agreement 
with  the   256  kcal/mol  determined  by  experiment  (Lossing, 
1972).   A  second  and  potentially  more  important  CgHg^ 
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Figure  4.1.   The  Most  Stable  C3H3   structures 
Cyclopropenyl ium  (4.1)  and  Propargylium  (4.2) 
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structure  Is  that  of  the  linear  propargyllum  Ion  with  a 
calculated  (Radom  et  al . ,  1976)  heat  of  formation  31  to  34 
kcal/aol  higher  than  that  of  the  cyclopropenyllum  ion,  in 
fair  agreement  with  the  25  kcal/mol  difference  found 
experimentally  by  Lossing  (1972).   Recent  calculations  have 
shown  several  other  stable  isomeric  structures  of  the  C3H3* 
Ion  with  higher  heats  of  formation  relative  to  the 
cyclopropenyllum  and  propargyllum  ions. 

The  importance  of  different  precursors  in  affecting  the 
reactivity  of  C3H3*  ions  was  reported  in  a  thermochemical 
study  by  Holmes  and  Lossing  (1979).   In  an  ICR  study  of 
C3H3'*^  reactions,  Ausloos  and  Lias  (1981)  showed  that 
significant  fractions  of  the  linear  isomer  can  be  produced 
by  charge  transfer  reactions  of  small  ions  (Ar  ,  Xe  ,  CO  , 
Ne"^,  etc.)  with  propargyl  chloride  and  bromide.   Later  it 
was  reported  (Baykut  et  al . ,  1986)  that  even  higher 
proportions  of  the  propargyllum  isomer     relative  to  the 
cyclopropenyllum  Isomer  can  be  obtained  with  propargyl 
iodide  either  by  electron  impact  or  charge  exchange  using 

xe^ 

A  study  by  Smyth  et  al.(1982)  demonstrated  that  the 
cyclopropenyllum  ion  was  relatively  unreactive  toward  simple 
hydrocarbon  fuels  whereas  the  linear,  propargyllum  ion  was 
quite  reactive.   In  particular,  reaction  of  propargyllum 
ions  with  acetylene  was  reported  to  produce  C5H3'*"  and  C^Hg* 
Ions  with  an  overall  C3H3'^  disappearance  rate  coefficient  of 
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1  X  10  '  cm''/s.   In  the  work  discussed  In  this  chapter, 
results  substantially  different  from  those  reported  earlier 
were  found  and  the  research  was  thus  expanded  to  Investigate 
possible  production  of  C5H5'^  Ions  from  ionic  sources  other 
than  C3H3'''  present  in  the  reaction  media  using  the  ejection 
capabilities  of  Fourier  transform  Ion  cyclotron  resonance 
(FTICR)  mass  spectrometry.   Reactions  of  propargyllum  ions 
with  C2D2  were  also  studied  in  order  to  delineate  further  a 
proposed  mechanism  for  the  C3H3'*'/acetylene  Interaction. 

Dlacetylene  is  another  Important  flame  neutral  which 
has  been  postulated  to  react  with  63113'''  ions  in  an 
ion/molecule  soot  formation  mechanism  (Olson  and  Calcote, 
1981a).   Reaction  pathways  and  the  rate  coefficient  for  the 
reaction  of  propargyllum  ions  with  dlacetylene  near  room 
temperature  were  thus  also  investigated.   In  this  chapter, 
results  for  the  reactions  of  C3H3'''  ions  with  acetylene, 
deuterated   acetylene,  and  dlacetylene  are  reported  and 
discussed  in  relation  to  previous  work  involving  C3H3''' 
reactions  (Ausloos  and  Lias,  1981;  Smyth  et  al . ,  1982; 
Baykut  et  al . ,1986) . 


Experimental 

Reactive  and  non-reactive  C3H3''"  ions  were  produced  by 
charge  transfer  reactions  of  various  precursors  with  Xe'*' , 
formed  with  an  electron  beam  pulse  of  5  ms  duration  at  an 
ionizing  electron  energy  of  ca.  15  eV.   Propargyl  iodide  was 
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used  as  a  precursor  for  C3H3*  In  studies  of  the  reaction  of 
this  Ion  with  acetylene,  deuterated  acetylene,  and 
dlacetylene,  while  a  number  of   different  precursors 
(propargyl  Iodide,  propargyl  bromide,  propargyl  chloride, 
propyne  and  allene)  were  used  to  Investigate  the  reaction 
mechanisms  leading  to  CgHg"^  ion  formation.   In  some  studies 
C3H3'^  Ions  were  produced  directly  from  the  above-mentioned 
precursors  by  electron  Impact  to  compare  the  effect  of 
Ionization  technique  on  the  ratio  of  reactive  to  unreactive 
Isomers . 

Propargyl  Iodide  was  prepared  from  propargyl  chloride 
by  a  halogen  exchange  reaction  (March,  1977).   The  details 
of  purification  are  given  elsewhere  (Baykut  et  al  .  ,  1986). 
Propargyl  chloride,  propargyl  bromide,  allene,  propyne  and 
acetylene  were  obtained  commercially  and  their  purity  was 
checked  by  obtaining  wide  mass  range  spectra.   In  the  case 
of  acetylene,  some  production  of  protonated  acetone  was 
observed  at  long  delay  times  indicating  the  presence  of  a 
small  amount  of  acetone  as  an  impurity.   Propargyl  bromide 
was  distilled  before  use  in  order  to  remove  toluene  which 
was  present  as  a  stabilizing  agent.   Deuterated  acetylene 
was  prepared  from  DgO  and  CaCg .   Dlacetylene  was  prepared  by 
the  method  of  dehydrochlorlnatlon  of  1 , 4-dichloro-3-butyne 
in  aqueous  potassium  hydroxide/dioxane  solution  (Snow, 
1985).   All  the  samples  were  used  after  multiple 
f reeze-pump-thaw  cycles. 
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,  Results 

C3H3   Reactions  with  Acetylene.   Despite  an  earlier 
report  (Smyth  et  al .  ,  1982)  that  C3H3'^  Is  quite  reactive 
with  acetylene,  only  very  low  intensities  of  C5H3'*'  and  CgHg* 
produced  via  this  reaction  could  be  found  In  this  work. 
Experimental  conditions  of  the  earlier  study  were  duplicated 
as  closely  as  possible,  and  then  varied  substantially  with 
respect  to  relative  pressures  of  neutrals  (from  1:1  to  8:1 
^2^2"^3^3^^  and  overall  system  pressure  (from  5  X  10    to  3 
X  10"^  torr).   The  C3H3"''  Ions  were  formed  from  propargyl 
chloride,  bromide,  and  iodide  by  both  charge  transfer  using 
Xe"*"  and  electron  impact.   In  order  to  determine  other 
possible  sources  of  CgHg"*"  observed  under  the  earlier 
reaction  conditions,  binary  mixtures  of  acetylene  and  one  of 
the  €3113'^  precursors  reported  (Ausloos  and  Lias,  1981;  Smyth 
et  al  .  ,  1982;  Baykut  et  al.,  1986)  earlier  were  used. 
Intensities  of  CgHg"^  and  C5H3'*"  were  first  measured  after  a 
125  ms  reaction  time.   Then  the  parent  ion,  C2H2^,  and  C3H3"*" 
were  each  ejected  separately  during  the  125  ms  reaction 
period  to  assess  their  contribution  to  CgHg"^  and  C5H3'*' 
formation. 

For  each  different  neutral  precursor,  the  sources  and 
amounts  of  C5H5*  ions  produced  were  found  to  be  different. 
Propyne  and  allene  were  similar  in  producing  large  amounts 
of  CgHg^  and  no  CgHg^  ion.   However,  the  C3H3*  +  CgHg 
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reaction  was  not  responsible  for  C5H5*  formation.   The  main 
reactions  leading   to  C5H5'^  were 


C3H4+  +  CoHp  >  C.H^*  +  H 


'2"2 


5"5 


(4.1) 
(4.2) 


In  both  cases.   On  the  other  hand,  when  propargyl  Iodide, 
propargyl  bromide  and  propargyl  chloride  were  used  as 
precursors,  relatively  smaller  amounts  of  CgHg^  Ion 
formation  were  observed  along  with  C5H3'^  Ion  production.   In 
the  propargyl  chloride  case 


C3H3CI   +  CgHg 


>  CgHg^  +  CI 


C^Hp"^  +  C^H-,C1 >  CcHc"^  +  CI 


'3"3 


5"5 


(4.3) 
(4.4) 


were  the  major  reactions  leading  to  CgHe"^  formation.   For 
propargyl  bromide  the 


CaHgBr*  +  CgHg 


>  CgHg^  +  Br 


(4.5) 


reaction  was  the  only  source  of  CgHg"^  Ions  observed.   Any 
contribution  to  CgHg"^  formation  from  linear  C3H3'^  was  less 
than  the  experimental  uncertainty.   Finally,  very  little 
(almost  negligible)  amounts  of  CgHg"^  Ions  were  observed  when 
propargyl  iodide  was  used  as  a  precursor  and  the  reactions 


?ifefe 
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C^H^"^    +    C<,H<>I >    CrHc"^    +    I 


'2"2 


'3"3 


5"5 


C3H3       +    CgHg  >    C5H5 


(4.6) 
(4.7) 


were  the  major  contributors  In  this  case.   An  upper  limit 
for  the  rate  constant,  k,  for  reaction  (4.7)  was  estimated 
as  5  X  10"^^  cm^/s  by  assuming  that   the  very  small  CgHg"*" 
signal  observed  resulted  from  this  reaction,  and  using  the 
expression  [C3H3'"](t)  =  [C3H3^](0)  -  [CgHg^Kt)  = 
[CgHg^] (0)e""^*,  where  n  Is  the  CgHg  number  density. 
Overall  results  for  the  production  of  C5H5*  and  CgHg"^  Ions 
in  mixtures  of  acetylene  and  various  neutrals  used  as  the 
precursors  of  C3H3*  are  summarized  in  Table  4.1. 

Because  propargyl  Iodide  was  shown  to  produce  the 
highest  reactlve/unreactlve  ratio  of  C3H3'^  Ions  both  In 
earlier  (Baykut  et  al . ,  1986;  Holmes  and  Losslng,  1979)  and 
the  present  work  (see  Table  4.2),  It  was  used  as  a  precursor 
for  03113"^  Ions  In  these  reaction  kinetics  studies.   Since 
the  precursor  neutral  molecule  was  always  present  in  the 
FTICR  analyzer  cell,  It  was  a  competitor  with  the  reactant 
neutral  of  Interest  in  ion/molecule  reactions  involving 
C3H3'^.   In  order  to  determine  the  rate  coefficient  of 
reactions  of  CgHg"*^  with  the  neutral  reactant,  It  was 
necessary  first  to  monitor  the  reactions  of  this  ion  with 
C3H3I  and  then  to  subtract  the  rate  coefficient  for  this 
reaction  from  the  total  rate  coefficient  observed  in  the 
presence  of  both  the  precursor  neutral  and  the  reactant  of 
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TABLE  4.1 


Production  of  C^Hc* 


;;Mous'Se"„?J.Ji-i;/-.?5;3;jo„3  .„  „.«„e3  o. 


Propargyl 
todlde 


Neutral 

Ionic  sources^ 

of  CgHg"^  after 

Xe*  charge 

transfer 

Ionization  of 

a  mixture  of 

the  neutral 

and  acetylene 

Ratio  of 
C5H5*  prod, 
relative  to 

that  in 
allene  case*^ 

Ionic  sources'* 

Intensity 

°'  C5H5* 
vs.  C5H3* 

Percent 
reactive 
C3H3* 
ions 

Allene 

C2H2*(405) 
C3H4+(605) 

1.0 

- 

- 

<5 

Propyne 

C2H2*(40S) 
C3H4*(605) 

0.75 

- 

- 

30 

Propargyl      C2H2*( 40-505) 
chloride      C3H3C1*(50-60S) 


0.25 


C3H3C1+(20X) 
C3H3+(20%) 
C2H2*(605) 


3.0 


Propargyl  C3H3Br*(90-100%)  0.O8 

broffllde 


C2H2"'(70X) 

[C3H3*  * 

C3H3Br*-l(305) 


2.0 


C2H2*(40i) 
C3H3*"(60%) 


<0.02 


C3H3*(50%) 


1.7 


IS 


85 


'.)0 


^   ?erri^?f    •'ere    produced    by    chemical    Ionization    charno    transfer    r  T ' ' 

Percentages    show    the    relative    contributions    to   C   1.    '    .nnf-*^ '"''"'    ^<^    • 
deternxned    by   double    resonance   expe^J^^lts    and    Sv^   an    est^l^il   ed'"''""^^"    - 

c   „      -10-.  estimated    uncertainty 

the'loi.!!,?'''"'^    ^^^    '^'^'^    '''°    ^-""^    pressure    (7    X    10-     torr,     ,- 


w 
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f^.-^- 
# 


TABLE  4.2 


Percentages*  of  reactive  C3H3"*'  found  from  various  precursors 
by  various  Ionization  techniques  (monitored  by  observing 
reaction  with  the  precursor  neutral). 


precursor 


Ionizing  Propargyl    Propargyl    Propargyl 

Technique  iodide      bromide     chloride 


Electron  impact  (15eV)        90  40  10 

Chemical  ionization  charge    90  85  15 

transfer  with  Xe"*" 


^Estimated  error  is  +5%. 
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Interest.   Reactions  of  03113'*^  with  propargyl  Iodide  were 
monitored  as  a  function  of  time  following  charge  transfer 
chemical  ionization  of  C3H3I  by  Xe"^  and  ejection  of  all  ions 
but  C3H3'^  from  the  analyzer  cell.   Results  obtained  were 
identical  to  C3H3'^  reaction  channels  with  propargyl  iodide 
which  have  been  reported  elsewhere  (Baykut  et  al . ,  1986). 


Isomerization  of  Linear  C^^^'*' .       In  addition  to  the 
absence  of  any  significant  C5H3'*^  and  CgHg"*^  formed  by 
reaction  of  linear  €3113'^  with  C2H2.  it   was  also  observed 
that  C2H2  led  to  the  isomerization  of  linear  C3H3*  ions  t 
their  cyclic  form,  thus  rendering  them  unreactive  toward 
their  parent  neutral  (C3H3I)  as  well  as  toward  C2H2.   This 
Isomerization  was  followed  as  a  function  of  C2H2  pressure 
and  a  direct  pressure  dependence  was  found,  as  can  be  seen 
in  Figure  4.2.  .  * 


o 


■  • '  ".if-"',  ■- 


C2H2   Reactions  with  0203-   To  achieve  a  better 
understanding  of  the   isomerization  of  linear  C3H3''' ,  C2D2 
instead  of  CgHg  was  used  as  the  neutral  reactant.   The 
following  isotope  exchange  reactions  were  observed: 


CgHg"^  +  CgDg >  C3H2D"^  +  C2DH  (4.8) 

^Z^Z*    +  C2D2  >  C3HD2^  +  CgHg  (4.9) 

C3H2D'^  +  CgDg >  C3HD2*  +  C2DH  (4.10) 

C3H2D*  *  C2D2  >  C3D3*  +  C2H2  (4.11) 
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Deactivation  of   I- 031-1^    by    acetylene  at  different    pressures 


Figure  4.2. 
Pressures  of 


^3^3 


Ions  at  Different 


■2"2 


•3"3 


Isomerizat ion  of  Linear 

CoHo.   CoH-i*  ions  were  produced  by  charge 

p(Xe)  was  adjusted  to  maintain  a  constant  total  pressure  of 
2.6  X  10~"  torr  as  measured  on  the  ionization  gauge.   (All 
pressures  are  capacitance-manometer  corrected.) 


transfer  reactions  with  Xe.   p(C.,H.,I)  =  1.1  x  10'  torr: 


h^^ 
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CgHDg"^  +  CgDg >  C3D3*  +  CgDH  (4.12) 

Using  FTICR  ejection  capabilities,  It  was  found  that 
reactions  (4.9)  and  (4.10)  contribute  equally  to  the 
production  of  C3HD2^  while  reaction  (4.12)  produces  more  of 
0303"*^  (80%)  than  reaction  (4.11)  (20%). 

Ion  intensity  vs.  time  curves  for  the  C3H3'*'/C2D2  "' 
reaction  are  shown  in  Figure  4.3.   The  overall  rate 
coefficient  for  the  disappearance  of  C3H3"*'  was  calculated  by 
subtracting  the  observed  rate  coefficient  for  the  reaction 
with  propargyl  Iodide  from  the  total  observed  rate 
coefficient  in  the  presence  of  C2D2.   This  observed  rate 
coefficient  was  then  converted  to  the  true  rate  coefficient 
using  the  corrected  pressure  of  C2D2.   A  value  of  (4.5+1.9) 
X  10"^   cm^/s  was  found  at  a  cell  temperature  of  373  K  for 
the  disappearance  of  03113'^  (reactions  (4.8)  and  (4.9)).   In 
Figure  4.4,  ion  intensity  vs.  time  curves  of  C3H3'*'  are 
compared  for  reactions  with  and  without  C2D2. 

CgHg"*"  Reactions  with  Dlacetylene.    After  ejection  of 
all  ions  except  C3H3'*'  following  charge  transfer  chemical 
ionization  by  Xe"*"  of  a  mixture  of  dlacetylene  and  propargyl 
iodide,  the  ion/molecule  reactions  as  a  function  of  time 
were  monitored.   Consecutive  C2  and  C^H2  addition  reactions 
were  observed: 

CgHg""   +   C4H2 >  CrjH^^  (4.13) 
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Reactions  of  I-C3H3  with   C2D2 


80 


+ 

•     C3H; 

D 

C3H2D* 

0 

C3HD2 

A    • 

C^D^ 

1.5 


Time  /s 


C2D2 


Figure  4.3.   Isotope  Exchange  Reactions  of  03113"^  with 
Disappearance  of  03113'^  ion  includes  reactions  with  propargyl 
iodide.   Note  that  the  sum  of  all  isotooic  forms  of  C3H3"*" 


equals  the  total  unreactive  C3H3 


sum  of  all  isotopic  forms  of 
on  w 
when 


remaining  at  the  end  of  the  reaction  with  CoDo  approximately 


neutral  reactant  at 
p(C3H3l)  =1.1  X  10" 
=6.2x10^  torr. 


the  same 

7 

torr  ; 


pressure 
p(C2D2)  = 


CgHg 


2'^2 
is  used  as 


(see  Fig.  4.2). 
1.2  x  10"^  torr: 


p(Xe) 
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Reactions    of    L-C,H^  with    C,H,I      and    CoQ 


3  "3 


0 


'2^2 


+    :  C3H3  with    C3H3I 

D  :  C3H3  with  C3H3I  and  C2D2 


J — I — I I I  '  I  I   I  I 


0 


0.5 


1^ 


I.O 
Time/s 


-1 1 1 1 I     r    I     I 

1.5         2.0 


Figure  4.4.   CgHg^  Ion  Decay  Curves  for  Reaction  with  C^H  I 
and  CgDg.   (Pressures  are  the  same  as  piven  for  Fig.  1.3.5 
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C3H3       +    C^Hg    >    C5H3       +    CgHj 

CcHo  +      C/iHo >      CnHi 


'5"3 


'4"2 


'9"5 


CrHo       +    C^Ho ^    C17H0       +    CoH 


'5"3 


'4"2 


7"3 


'2"2 


C7H3    +   C4H2 >   CjjHg 

CgHg^   +   C4H2   >   CjiHg^   +   C2H2 


(4.14) 
(4.15) 
(4.16) 
(4.17) 
(4.18) 
(4.19) 


Some  of  these  product  ions  were  seen  to  react  further 
with  propargyl  Iodide  by  displacement  of  atomic  Iodine: 


^5^3  *  ^3*^3^  ~ 

^7^3  "^  ^3^3^  " 

C7H5  +  C3H3I  - 

^9^5  "^  ^3^3^  ~ 


->  CqHq*      +   I 

— >  Cl0"6^  *  I 


-->  CjoHs   *  I 


--->  Ci2H8^  +  I 


(4.20) 
(4.21) 
(4.22) 
(4.23) 
(4.24) 


Ion  Intensity  vs.  time  curves  for  the  C3H3'^/C4H2 
reaction  are  shown  in  Figure  4.5.   The  rate  coefficient  for 
the  disappearance  of  C3H3^  (reactions  (4.13)  and  (4.14), 
Figure  4.6)  was  calculated  as  described  earlier,  and  a  value 
of  k  =  (1.4  +  0.7)  X  10"^  cm^/s  was  found.   Propargyl 
iodide,  bromide,  and  chloride  were  all  used  as  precursors  of 
C3H3*  in  studying  its  reactions  with  diacetylene.   For  each 
precursor,  both  electron  impact  and  charge  transfer  chemical 
ionization  techniques  were  used.   The  percentages  of 
reactive  isomer  in  the  reaction  with  diacetylene  are  shown 
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Reactions    of    l-C,H^with   C^Ho  < 


3"  '3 


t^:i. 


+    : 

C3H3 

a   : 

C5H5        >: 

O  : 

C7H3 

0    : 

C9H5 

X    : 

C11H5 

Figure  4.5.   Reactions  of  C3H3*  with  C^Hg.   Disappearance 
CgHg"*"  and  product  ions  include  reactions  with  propargyl 
iodide.   p(C3H3l)  =  1.1  x  10"*^  torr:  p(C4H2)  =  4.8  x  10"'^ 
torr;  p(Xe)  =  6.2  x  10~"  torr.   (All  pressures  are 
capacitance-manometer  corrected. ) 


of 


IV 
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Reactionsof  I-C3H3  with  C3H3I  and  C4H 


f  rO 


+  :  C3H3  with  C3H3I 

□  :  CjHjWith  C3H3I  and  C4H2 


—J — I — ' — I — I — I — I — I — I — I — I 1 I I      1     I      I     ■      I 

0  1.0  2.0 

Time/s 


Figure  4.6.   C3H3'"  Decay  Curves  for  the  Reactions  with  C3H3I 
and  C4H2.   (Pressures  are  the  same  as  given  for  Fig.  4.5.) 
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TABLE  4.3 


Percentages®  of  reactive  C3H3'*"  observed  in  the  reaction  with 


C4H2 


precursor 


Ionizing  Propargyl    Propargyl    Propargyl 

Technique  iodide      bromide     chloride 


Electron  impact  (15eV)        75  30  5 

Chemical  ionization  charge    75  65  5 

transfer  with  Xe* 


'Estimated  error  is  +5%      ^PCC^Hg)  =  4.8  x  10"''  torr 
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In  Table  4.3.   When  these  percentages  of  reactive  Isomer 
were  compared  to  those  in  the  absence  of  C^H,  (see  Table 
4.2),  It  was  clear  that  some  Isomerlzatlon  of  the  reactive 
linear  C3H3'^  Ion.  as  well  as  reactions  (4.13)  and  (4.14). 
had  taken  place  (see  also  Fig.  4.6).   This  isomerlzatlon  was 
followed  as  a  function  of  C^Hg  pressure  and  a  direct 
pressure  dependence  was  found,  as  can  be  seen  In  Table  4.4. 

The  reactions  of  l-CgHg"^  with  propargyl  iodide  and  with 
both  acetylene  and  dlacetylene  have  also  been  followed  at 
several  elevated  temperatures  up  to  500  K.   All  the  rate 
constants  were  found  to  be  similar  to  their  room  temperature 
value  within  experimental  error. 

Discussion 
Effect  of  Different  Precursors.   Different  percentages 
of  reactive  CgHg"^  were  found  from  three  different 
precursors,  propargyl  iodide,  propargyl  bromide  and 
propargyl  chloride  as  shown  In  Table  4.2.   To  explain  the 
differences  observed  in  reactivity,  schematic  potential 
energy  surfaces  for  these  precursors  are  shown  in  Figure 
4.7.   Experimental  thermochemlcal  data  reported  by  Holmes 
and  Losslng  (1978)  were  used  in  the  generation  of  the 
potential  surfaces.   Reverse  activation  energies  for  c-CgHg"*" 
formation  from  propargyl  bromide  and  chloride  were 
determined  by  the  difference  between  the  experimental  and 
the  calculated  appearance  potentials.   Since  the  appearance 


I'^f 
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TABLE  4.4 


hanges  In  CoHo*  reactivity*  at  different  pressures  of 
iacetylene . " 


-7 


C^Hg  pressure/10  '  Torr 


*  unreactive  €3113'*^ 


0.8 
1.6 
4.8 
7.2 
8.0 
9.6 


16 
17 
25 
32 
35 
40 


^1-03113   Ions  were  produced  from  propargyl  Iodide  by 


chemical  Ionization  charge  transfer  with  Xe"*" .   (p(C3H3l)  = 

1.1  X  10"'^  torr;  Pj^g  was  adjusted  to  maintain  a  constant 

total  pressure  of  2.6  X  10~°  toj 

gauge) . 

°A11  pressures  are  capacitance-manometer  corrected 


irr  as  read  on  the  Ionization 
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I*.63  e/ 
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12.13  ev  z;^  ' 


ik 
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Figure  4.7.  Schematic  Potential  Energy  Surfaces  for 

C3H3  /C3H3X  System  from  Different  Precursors.   Propargyl 

Iodide  (a),  Propargyl  Bromide  (b).  and  Propargyl  Chloride 
(c). 


V  >. 


A-v 
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potential  of  C3H3'*"  produced  from  propargyl  iodide  very 
closely  corresponds  to  the  calculated  threshold  for  l-CgHg* 
rather  than  c-CgHg"^,  the  dissociation  to  the  latter  is 
assumed  to  have  a  significant  energy  barrier.   Thus,  reverse 
activation  energy  for  the  dissociation  channel  giving 
unreactive  c-CgHg"^  decreases  in  the  order  Ej^^jq  >  E{jj.qj,q  > 
^chloro  ^^  shown  in  the  figure.   Production  of  almost 
exclusively  reactive  l-CgHg^  by  both  electron  impact  and  Xe* 
charge  transfer  ionization  from  propargyl  iodide  suggests 
that  EjqjJq  is  so  large  that  the  fragmentation  channel 
leading  to  l-CgHg"*"  becomes  the  lowest  energy  channel.   In 
the  case  of  propargyl  bromide,  there  is  enough  excess  energy 
to  dissociate  to  both  l-CgHg"^  and  c-CgHg"^.   Production  of 
80*  l-CgHg"^  by  Xe"^  chemical  ionization  suggests  no 
significant  energy  barrier  for  the  1-CgHg"*^  channel.   It  is 
interesting  to  note  that  dissociation  to  l-CgHg'*'  reduces  by 
a  factor  of  two  when  electron  impact  ionization  is  used, 
which  demonstrates  the  effect  of  a  large  distribution  of 
electron  energies  from  an  electron  impact  ionization  source 
on  the  relative  abundances  of  two  channels  almost  equally 
accessible  energetically.   The  very  small  reactive  l-CgHg"^ 
percentage  produced  from  propargyl  chloride  by  both  electron 
impact  and  Xe"*^  charge  transfer  ionization  suggests  at  least 
a  small  energy  barrier  for  this  dissociation  channel,  as 
Indicated  schematically  in  Figure  4.7c. 
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Internal  Energy  of  the  Ions  in  Relation  to   Rate 
Coefficient  Measurements.   In  studying  reactions  it  is 
desirable  to  have  knowledge  of  the  internal  energy 
distribution  of  the  reactant  ions.   The  reactions  studied 
here  are  bimolecular  addition  reactions  followed  by 
unlBolecular  decomposition.   As  shown  in  Figure  4.7a,  when 
I'^S^S*  is  formed  from  propargyl  iodide  by  Xe"^  charge 
transfer  chemical  ionization,  1.5  eV  of  excess  energy  is 
available.   Much  of  this  excess  energy  will  be  converted 
into  translational  motion  of  the  heavy  Xe  and  I  neutrals 
resulting  from  the  charge  transfer  in  the  collision  process. 

—  6 

Under  typical  experimental  conditions  (ca.  2  -  3  x  10    torr 

total  pressure),  about  125  ms  was  allowed  for  the  charge 
transfer  process  and  for  ejecting  intermediates.   For  these 
conditions  the  C3H3*  Ions  collided  a  number  (10-15)  of  times 
with  the  excess  Xe  present  in  the  FTICR  cell,  leading  to 
near  thermal Ization  of  internal  energy  before  the 
ion/molecule  reactions  were  monitored.   Since  the  reaction 
time  scale  was  on  the  order  of  seconds,  it  can  be  assumed 
that  any  slight  initial  deviation  from  Boltzmann  behavior 
presents  no  serious  error.         On  the  other  hand,  the 
observation  that  the  rate  constants  for  the  reactions  of  1- 
CgHg"^  with  propargyl  iodide  and  with  both  diacetylene  and 
acetylene  are  temperature  Independent  implies  that 
thermalization  of  the  ions  is  not  complete  under  the 
conditions  reported  and  there  is  still  some  internal  energy 
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In  1-C3H3''",  which  Is  comparable  In  nagnltude  to  that 
contributed  from  the  range  of  temperatures  studied. 


Reactivity  of  l-CgHg"^  with  Acetylene.   Although  the 
results  of  03113'^  +  C2H2  reaction  are  not  In  agreement  with 
the  earlier  report  (Smyth  et  al . .  1982)  of  C3H3*/C2H2 
reactivity,  the  discrepancy  is  most  likely  due  to 
limitations  of  the  older  pulsed  ICR  (Smyth  et  al . ,  1982) 
Instrumentation  for  studying  ion/molecule  reaction  pathways 
in  complicated  systems  when  compared  to  newer  FTICR 
capabilities.   Facile  ejection  of  all  ions  except  the  one 
whose  ion/molecule  reactions  are  being  investigated  offers  a 
very  clean  monitoring  opportunity  for  product-parent 
relationships  even  in  complicated  consecutive  and 
competitive  reaction  systems.   Various  alternative  pathways 
for  the  production  of  C5H3'^  and  CgHg^  which  have  been 
described  above  probably  contributed  significantly  to  the 
Intensities  of  these  ions  seen  in  the  earlier  work. 
Additional  support  for  the  low  reactivity  of  C3H3'*^  with  CgHg 
is  found  in  a  recent  report  (Anlclch  et  al . ,  1986)  of  the 
rate  coefficient  for  this  reaction  as  less  than  0.01  X  10~® 
cm'^/s,  although  the  Isomeric  form  of  03113"^  was  not  given. 
It  Is  also  possible  that  the  highest  pressures  used  in  this 
work  did  not  reach  those  of  the  earlier  study  due  to 
differences  in  the  location  of  capacitance  manometers, 
Ionization  gauges,  etc.   Thus  third  body  stabilization  of 
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^5^5   collision  complexes  might  have  been  occurring  to  some 
extent  In  the  earlier  work  and  not  In  that  reported  here. 
In  fact,  such  colllslonal  stabilization  of  the  association 
complexes  for  the  reactions  of  C3H3"^  and  C4H4*  with  CgHg  has 
been  shown  to  occur  In  higher  pressure  SIFT  studies  (Smith 
and  Adams,  1987;  Knight  et  al . ,  1987). 

The  most  likely  mechanism  of  the  observed  Isomerlzatlon 
of  C3H3   Ions  by  collisions  with  acetylene  Is  a  "reactive" 
rather  than  a  "non-reactive"  one.   That  Is,  It  results  from 
an  Intimate  encounter  of  the  Ion  and  neutral  in  the   CrHk"^ 
collision  complex.   This  hypothesis  Is  confirmed  by  the  fact 
that  deuterated  forms  of  C3H3*  were  produced  when  CgDg  was 
the  neutral  reactant  (see   Figure  4.4).   Kinetic  modeling 
studies  (discussed  In  Chapter  5)  indicate  that  in  some  cases 
the  CgHg"^  collision  complex  dissociates  to  give  the  cyclic, 
unreactlve,  C3H3'^  isomer.  Instead  of  the  reactive,  linear 
form  which  reacted  initially.   The  possibility  of 
non-reactive  colllslonal  isomerlzatlon  of  linear  C3H3'^  to 
the  cyclic  Isomer  has  been  ruled  out  because  experiments  at 
elevated  pressures  of  xenon  (to  ca .  1  X  10"^  torr)  showed  no 
Interconversion.   Similar  interconversion  of  C^H^"^  ions  from 
a  linear  to  cyclic  form  has  also  been  reported  (Jarrold  et 
al.,  1984)  in  the  reaction  with  CgHg  and  has  also  been  shown 
to  take   place  via  complex  formation  by  using  Isotopically 
labeled  CgHg.   To  confirm  the  hypothesis  that  energetically 
less  stable,  reactive,  (linear)  C3H3''  ions  interconvert  to 


71 


more  stable,  unreactlve  ones,  cyclic  C3H3'*'  Ions  were  also 
reacted  with  C2D2  and  no  Isotope  exchange  reactions  were 
observed. 


Reactivity  of  l-CgHg"^  with  Dlacetylene.   Plots  of  C3H3* 
ion  Intensity  vs.  time  for  reaction  with  dlacetylene  (C4H2) 
(Figure  4.5)  Indicate  a  10-12%  Increase  In  the  Intensity  of 
the  unreactlve  Isomer  relative  to  the  reaction  when  the   ^ 
parent  precursor  only  Is  present.   Isomerlzatlon  of  reactive 
C3H3   was  also  seen  when  different  precursors  were  used 
(compare  Tables  4.2  and  4.3).   A  similar  mechanism  Involving 
complex  formation  may  be  responsible  for  this  Isomerlzatlon 
as  well,  although  It  was  not  investigated  in  any  detail. 


CHAPTER  5 
KINETIC  MODELING  OF  THE  REACTIONS  OF  CgHg"" 

Introduction 


As  reported  in  the  last  chapter,  binolecular  reactions 
of  the  propargyllum  form  of  C3H3'''  with  acetylene  aost  often 
result  In  an  Isomerlzatlon  to  the  cyclopropenyllum  Isomer. 
To  help  understand  this  Isomerlzatlon  process,  C3H3''' 
reactions  with  deuterated  acetylene  were  Investigated. 
These  studies  showed  that  the  Isomerlzatlon  proceeds  via  the 
CgHg"*"  ion/molecule  reaction  complex,  which  Is  sufficiently 
long-lived  under  the  experimental  conditions  employed  that 
deuterium  exchange,  as  well  as  isomerlzatlon,  takes  place. 
Thus  with  time  the  reactive  propargyllum  CgHg"*"  isomer  is 
converted  to  both  reactive  and  unreactlve  species  containing 
one,  two,  and  three  deuterium  atoms.   There  is  no  evidence, 
either  experimental  or  theoretical,  that  the  propargyllum 
cation  converts  into  the  cyclopropenyllum  cation  in  the 
absence  of  the  C5H5"*"  reaction  complex.  In  order  to   better 
understand  the  Isomerlzatlon  which  converts  the  reactive 
to  the  unreactlve  form  of  C^Hg''' .  kinetic  modeling  studies  of 
the  Ion  intensity  vs.  time  curves  reported  in  Chapter  4  were 
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73 
carried  out . ^   It  was  also  hoped  that  fitting  procedures 
would  produce  improved  Ion/molecule  reaction  rate 
coefficients.   Quantum  mechanical  calculations^  on  C3H3"*'  and 
CgHg^  structures  and  reactivity  were  used  to  guide  the 
modeling  effort. 


Experimental 
A  mixture  (predominantly  propargyl  lum )  of  the  CgHg"*" 
Isomers  was  in  most  cases  formed  by  charge  transfer  from 
propargyl  Iodide  to  Xe"*",  produced  by  15  eV  electron 
Ionization  of  Xe  (present  at  pressures  >  lOx  those  of  other 
gases).   Other  conditions  such  as  neutral  partial  pressures, 
pulse  sequences  and  reaction  times,  and  the  sources  of 
chemicals  were  kept  as  close  as  possible  to  those  reported 
In  Chapter  4  for  the  duplicate  kinetics  experiments  reported 
and  modeled  here.   Any  significant  deviations  are  given  In 
the  text,  table  headings,  or  figure  captions.   All  pressures 
reported  In  this  chapter  were  determined  by  a  capacitance 
manometer,  and  then  multiplied  by  a  "system  factor"  of  0.30 
which  corrects  for  the  fact  that  the  pressure  read  by  the 
capacitance  manometer  Is  not  the  same  as  that  In  the  FTICR 


^The  kinetic  modeling  studies  were  performed  In  the 
Envlronlcs  Division  of  Air  Force  Engineering  and  Services 
Center,  Tyndall  Air  Force  Base,  Florida  by  F.  Wiseman  using 
multiple  experimental  data  sets  produced  at  Identical 
conditions  to  those  reported  In  Chapter  4. 

^A.  Cameron,  J.  Leszczynskl,  M.  C.  Zerner  and  B  Welner, 
submitted.    J.  Feng,  J.  Leszczynskl  and  M.  C.  Zerner,  submitted 
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cell.   Non-linear  least-squares  fitting  routines  employing 
Marquardt's  algorithm  (Annino  and  Driver,  1986),  Implemented 
on  two  different  computers^,  were  used  for  kinetic  modeling. 
Complete  analytical  solutions  were  obtained  from  the 
chemical  models  developed  below  for  the  systems  CgHg^  +  CgHg 
and  CgHg^  +  C^Hg .   A  complete  analytical  solution  was  not 
possible  when  an  isotope  effect  was  included  in  the  chemical 
model  for  the  CgHg"*^  +  CgDg  system.   Numerical  integrations 
used  the  finite  difference  method  (Annino  and  Driver,  1986). 

Results 


Models  of  CgHg^  +    C2H2  Reaction.   As  reported  in 
Chapter  4,  collision  of  the  propargylium  cation,  I-C3H3*, 
with  acetylene  forms  the  cyclopropenylium  cation,  c-CgHg"*", 
which  is  unreactive  on  the  time  scale  of  the  experiments, 
given  the  pressures  attainable  in  the  FTICR  cell. 
Experiments  with  C2D2  showed  that  an  encounter  complex  which 
allows  for  Isotopic  scrambling  is  formed.   Hence,  whatever 
the  isomeric  form  of  this  complex,  a  structure  having  the 
chemical  formula  C5H5"*"  can  be  postulated.   Since  no  species 
of  m/z  65  is  observed  in  the  mass  spectrum,  the  (CgHg"*") 
species  must  be  In  steady-state  and  of  low  concentration. 
The  simplest  scheme  which  takes  into  account  this 
information  is  given  in  Figure  5.1. 


^Tektronix  Model  4054  and  Hewlett  Packard  Model  150. 
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i-C3H3-'  +  C2H2  >  CsHs"" 

k 


i-C^H-,   +  C^H-jI 


.C5H5 


3"3   "  "'3"3 

+    I 


->    sink 


>  l-C^H-i      +   C-,H 


'3  "3 


^S^S^   — ~^   ^'^^3^2*   ■*"  ^■'" 


2"2 
2"2 


Figure  5.1.   Reaction  Scheme  Postulated  for  the  Kinetic 
Modeling  of  the  Reaction  of  linear  CgHg"^  with  Acetylene 


fe"i  ■ 
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Applying  steady-state  kinetics  to  this  scheme  yields 


I(t)  =  jO  -  k  '11°  (1 


e~^*)/» 


(5.1) 


In  which  1°  and  I(t)  are  Ion  Intensities  Initially  and  as  a 
function  of  time,  respectively,  Ij  refers  only  to  the  linear 
form. 

'^p'  =  '^pPCgHgl  (5-2) 

and 


^    =  tkfk^Pc2H2/(ki  -  k^)]  +  kp- 


(5.3) 


Equation  (5.1)  was  fitted  to  several  kinetic  runs  reported 
in  Chapter  4.   Table  5.1  shows  results  of  these  fits.   A 
plot  of  8    vs.  Pq    H   should  be  linear,  as  Implied  by  equation 
(5.3),  and  this  Is  demonstrated  In  Figure  5.2.   Results 
yield  kfkg/(ki  +  k^)  =  2.3(.2)  x  10^  torr~^s"^  and  kp'  = 
1  .4(0.2)  s"^ . 


Models  of  CgHg^  -*■  C2D2  Reaction.   Reaction  of  C3H3'*^  + 
C2D2  Is  complicated  by  the  observation  that  Isotoplc 
scrambling  occurs  and  Isotope  effects  are  possible.   Several 
models  were  tried.  Including  those  which  allowed  for 
complete  Isotoplc  scrambling  and  those  which  allowed  only 
partial  scrambling.   The  simplest  model  allowed  for  complete 
scrambling  and  no  Isotope  effects  . 


^Lampe  and  Field  (1959)  studied  the  reaction  of  CD^"^  + 
H4  ,  and  observed  the  following  yields:  C^HD^"^  :  1/10, 


'2"4 


C3H3D2*  =  2/5.  CgHgDg"^  :  2/5,  and  C3HD4'^  :  1/10.  Statistical 
yields  with  no  Isotope  effect  would  have  been:  C3H4D*  : 
1/14,  C3H5D2*  :  3/7.  C3H2D3"^  :  3/7.  and  CgHD^"^  :  1/14,  very 
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TABLE  5.1 


Results  of  fits  of  equation  (5.3)  to  kinetic  data  for  the 
CgHg"^  +  CgHg  reaction.^ 


Pf>  u  /torr 


kp'I°/Arb.  units  s"^ 


fl/s 


1 


4.4  X  10 


1.0  X  10 


1.6  X  10 


-7 


-6 


-6 


2.0  X  10' 


1  .53( .06) 
1  .65( .20) 
1  .74( .12) 
1  .71( .08) 
1  .71 ( .11) 


1.58( .08) 
2.40( .31) 
3.60( .24) 
4.82( .23) 
6.33( .38) 


^The  standard  error  of  estimate  computed  by  the  fitting 
program  is  shown  In  parentheses. 


close  to  the  observed  values.   Hence,  since  the  CqLq"*" 
complex  exhibited  almost  complete  scrambling,  it  is 
reasonable  to  expect  the  CgLg"*"  complex  modeled  here  to 
undergo  complete,  or  nearly  complete  scrambling. 


;,nr*T 
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,13 


.26       .39 


52       .65 


;■-  » 


Figure  5.2.   The  Plot  of  5  as  a  Function  of  Acetylene 
Pressure  (see  Equation  (5.3)). 
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Conplete  scrambling  occurs  when  fragmentation  of  the 
complex,  CgLg"*"  (L  =  H,  D),  yields  precursors  having  a 
statistical  distribution  of  hydrogen  and  deuterium  atoms. 
CgHgOg"*"  will  then  yield  the  following  ratios  of  precursors: 
CgHg"^  =  1/10,  CgHgD"^  =  3/5,  and  CgHDg"^  =  3/10.  Using  these 
statistics  for  obtaining  the  Isotoplc  distribution  in  the 
CgLg*  precursors,  we  obtain  the  scheme  shown  in  Figure  5.3. 
With  the  assumption  again  that  all  four  CgL^"^  complexes  are 
in  steady-state,  a  full  analytical  solution  Is  possible  for 
the  set  of  kinetic  differential  rate  equations.   The 
solutions  to  the  set   of  equations  are  given  in  Appendix 
II. A.   Figure  5.4  shows  the  best  fit  curves  to  a  typical 
data  set.   Table  5.2  shows  the  fitted  parameters,  errors, 
and  residual  sum  of  squares  from  fits  of  a  typical  data  set. 

An  examination  of  Figure  5.4  shows  that  the  model  given 
by  Figure  5.3  does  not  adequately  explain  the  production  of 
CgHDg^.   In  an  attempt  to  examine  this,  incomplete,  or 
partial  Isotoplc  scrambling  was  next  assumed.   To  do  this 
correctly  requires  a  detailed  knowledge  of  the  chemistry  of 
the  system,  which  is  not  available.   A  somewhat  crude 
application  of  Isotope  effects  applies  multiplicative 
factors  to  the  Individual  rate  constants  and  this  procedure 
requires  but  a  nominal  knowledge  of  the  structure  of  the 
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l-CjH,*  ♦  CjOj  — ^>   CsHjDj*^ 


t-CjHj      +  C3H3I  — S->   sink 


S"3°2*  —->  To    ^-*^3«3*    *  To  ^2°2   *  T   ^^3"2°^  * 
I  S«°   *  To    ^^3«°2*   ^  To   ^2"2 


S"3°2*  -^>  lO    '=^3"3'   *  1¥  S°2    "  T  =-^3"2°*  * 


I  C^HD    .  ^3  C-C3HD,*   .  ^  C^H^ 


t-CjHjD*  +  C2D2  ^— >   C5H2D3* 


t-CjHjD*  +  C3H3I  — —>    sink 


CgHjD^*  _L->  i_   t^:n   +    ^  J^    ^   „       .3 


*  To  *-<^3°3     *  -^o  S"2  *  T 

f  *=2"''  *  Tn   t-C.ny  +  1^  CD 


1-C3HD2  * 


10   --3"2"      "To  ^2°2 
S"2°3*  — >  \o  =^3°3*  *  ^  S«2   *  r  <=-<^3«°/  - 


5-V°*T0=-<^3«2°      *foV2 


kf 
*"^3"°2*   *  *^2°2  >  ^5"°4* 


I-C3HD2*  ♦  C3H3I  — 2->   sink 


5      4 


'    -   I   ^■^3""2      *   5  '^2°2    *  f  *-<^3°3^   *  f  C^HD 


C.HD."  i->  ^   t-C   HD,*    +  i  C   D      ^  i    .^   -   *    .     2 


+        c    .    3 


••■    .    3 


CHT,.-    -^>  I  c^3„0^"  *  I  C^O^    .  1  0-0303-   .  I  .^., 


S      4 


+  —  C,HD 


^-^3°3      *S°2-^>V5 


I^^Oj       ♦   ^3^3^  ~^    ^'■'^ 


'3    3 


Figure  5.3.   Reaction  Scheme  Postulated  for  the  Kinetic 
Modeling  of  the  Reaction  of  linear  CgHg"^  with  Deuterated 
Acetylene  Assuming  Complete  Scrambling  and  No  Isotope 
Effects . 
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MODEL  FIT 

a       'C3H2D  + 

+       'C3H3  + 


1.5 


(a) 


h) 


Figure    5.4.       Model    Fit    (using    the    scheme    of    Figure    5.3)    to    a 
Typical    Data    Set    for    CgHg"^    +    CgDg    Reactions.       Ion    Intensity 

and 


vs 


C3D3 


time  curves  for  (a)  C3H3 


and  CgHgD" 


and  (b)  CgHDg 
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species  Involved^.   This  simple  scheme  was  applied  as 
described  below. 

Quantum  mechanical  calculations^  indicate  that  CgHg 
does  not  readily  react  with  the  cyclopropenyl ium  cation,  but 
does  react  with  the  propargylium  cation  (l-CgHg"^)  without 
barrier  with  formation  of  the  four  possible  products  shown 
in  Figure  5.5.   Thermodynamical ly ,  only  Structure  (5.1)  is 
stable  with  respect  to  decomposition  to  c-CgHg"^  +  CgHg. 

However,  since  the  l-CgHg^  +  CgHg  >  Structure  (5.1) 

reaction  is  at  least  60  kcal/mol  exothermic^  ,  In  the 
absence  of  stabilizing  collisions  this  energy  then  permits 
many  different  Isomeric  forms  to  be  sampled  before 
decomposition  back  to  CgHg^  +  CgHg.   If  the  sampling  of  all 
isomers  is  fast,  complete  isotopic  scrambling  is  expected. 
Given  the  uncertalnity  in  CgHg"^  structures,  total 
equivalence  of  carbon  atoms  in  the  complex  is  assumed.   For 
ease  of  understanding,  a  cyclic  CgHg"^  complex,  in  which  all 
C  atoms  are  sp^  hybridized,  might  be  visualized.   With  these 
assumptions,  it  is  a  simple  matter  to  determine  which  carbon 
atoms  undergo  hybridization  changes  during  the  reaction. 
For  instance,  in  the  attack  of  CgHg  on  I-C3H3*,  three  carbon 


Although  somewhat  crude,  applying  a  multiplicative  factor 
to  the  rate  constant  for  each  deuterated  site  is  in  keeping  with 
experimental  evidence.   For  instance,  in  the  acetolysls  of  some 
identical  tosylates,  each  deuterium  atom  substitution  changed  the 
rate  constant  by  ca.  0.84  ( Strei twei ser  et .  al . .  1958). 

Details  of  these  studies  are  given  in:  J.  Feng,  J. 
Leszczynski  and  M.  C.  Zerner,  submitted;  and  J.  Leszczinski,  M. 
C.  Zerner  and  F.  Wiseman,  submitted. 
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atoms  change  from  sp  to  sp^  hybridized.   Upon  fragmentation, 
some  carbon  atoms  remain  sp^  hybridized;  others  become  sp 
hybridized. 

In  general,  for  o-secondary  Isotope  effects,  an 
Increase  In  hybridization  In  going  from  the  reactant  state 
to  the  transition  state  yields  an  inverse  Isotope  effect  (kp 
<  Kq),  whereas  a  decrease  in  hybridization  yields  a  normal 
isotope  effect  (kjj  >  kp)  (Dreuth  and  Kwart,  1980).   In  this 
system  the  blmolecular  addition  reaction  will  have  an 
inverse  Isotope  effect,  and  the  fragmentation  reaction 
should  have  a  normal  Isotope  effect.  Figure  5.6  shows  how 
the  Isotope  effects  for  the  reaction  of  l-CgHg"*"  +  C2L2  can 
arise.   For  the  forward  addition  reaction,  Epjj  >  Epp  and  for 
the  fragmentation  reaction,  Ej^jj  <  Ej^p ,  in  which  the 
subscripts  "F"  and  "R"  refer  to  forward  and  reverse, 
respectively,  and  "H"  and  "D"  refer  to  protonated  and 
deuterated  species,  respectively.  Since  the  blmolecular 
addition  is  very  exothermic  for  the  formation  of  most  CgLg"*^ 
Isomers,  it  might  be  expected  that  the  "average"  transition 
state  structure  might  closely  resemble  the  reactants  and  not 
any  of  the  C5L5*  isomers.   This  in  turn  implies  that  E^p  - 
Epjj  >  Epjj  -  Epjj  (zero-point  effects).   However,  the  excess 
energy  in  the  reaction  will  allow  longer  sampling  times  for 
the  more  energetic  CgLg"^  Isomers.   The  less  energetic 
Isomers  which  are  sampled  will  undoubtedly  be  in  higher 
rotational  and  vibrational  levels.   The  overall  effect  is  to 
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lessen  the  normal  Isotope  effects  expected  In  the 
fragmentation  of  the  C^Lg"^  Isomers,  unless  there  are 
sufficient  collisions  to  stabilize  the  Isomers  prior  to 
fragmentation. 

In  the  model  development  outlined  below,  only 
a-secondary  Isotope  effects  will  be  considered  important. 
P-secondary  Isotope  effects,  arising  predominantly  from 
hyperconjugatlon,  can  sometimes  be  important  (Melander, 
1960),  but  will  be  assumed  here  to  be  minor  compared  to  the 
a-effects.   It  is  also  possible  that  hydrlde/deuterlde 
shifts  may  be  occurring  in  the  transition-  state.   However, 
hydride  transfers  often  exhibit  small  Isotope   effects 
(Melander,  and  Saunders.  1980)  and  hydride  shifts,  if  they 
occur  at  all,  will  be  assumed  here  to  give  negligible 
contributions  to  the  isotope  effects. 

Even  though  the  different  reactions  in  Figure  5.3  will 
have   different  Isotope  effects,  the  Introduction  of  an 
Independent  fitting  parameter  for  each  type  of  reaction  is 
not  justified.   Though  crude,  only  one  additional  parameter 
was  Introduced  into  the  model  to  account  for  all  potential 
o-secondary  Isotope  effects.   The  method  for  introducing 
this  parameter  Is  outlined  as  follows. 

If  a  carbon  atom  bearing  a  deuterium  atom  undergoes  a 
hybridization  change  from  sp  to  sp^  (force  field  becoming 
stronger),  the  "Isotope  effect  factor",  F,  Is  Introduced  as 
a  multiplicative  factor  in  the  rate  constant.   For  two 
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deuterium  atoms,  F^  is  the  multiplicative  factor,  etc.   If 
the  deuterium  atom  Is  attached  to  an  atom  changing  from  sp^ 
to  sp  hybridized,  the  rate  constant  Is  divided  by  F,  for  two 
deuteriums,  F*^,  etc.   Introducing  the  same  factor  for  both 
addition  and  fragmentation  reactions  Implies  a  constraint 
which  Is  at  best  only  qualitatively  correct.   Applied  In  the 
numerator  the  factor  corrects  for  a  single  deuterium  atom 
(F   for  two,  etc.)  attached  to  a  center  undergoing 
hybridization  change  from  sp  to  sp^  in  the  transition-state 
complex.   Applied  in  the  denominator,  it  corrects  for  a 
change  from  sp^  to  sp  hybridization.   Since  the 
transition-state  complex  has  a  stronger  force  field  than  the 
reactant  state  (l-CgLg"^  +  CgLg)  at  these  centers,  the 
a-secondary  isotope  effect  will  be  "inverse"  and  F  should 
therefore  be  greater  than  unity. 

Using  the  structural  notation. 


l-HDCgH^  =  [   C-C=C-H]^,  l-HgCgD^  =  [  ^C-C=C-D]^, 


H' 


Hv 


l-HDCgD^  =  [   C-C  =  C-D]-',  and  l-DgCgH^  =  [   C-C  =  C-H]*. 


D 


the  full  kinetic  scheme  is  shown  in  Figure  5.7.   The  neutral 
fragmentation  products,  CgHg.  CgHD,  and  CgDg.  have  not  been 
included  for  brevity.   It  is  assumed  that  I-C3L3*  undergoes 
the  same  kind  of  hybridization  changes  when  reacting  with 
C3H3I  as  it  does  with  CgLg.   As  shown  in  Chapter  4,  the 
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Figure  5.7.   Reaction  Scheme  Postulated  for  the  Kinetic 

Modeling  of  the  Reaction  of  linear  r  h  +  ..<..u  r,   *    ^  ^ 
,  ,  '^ii'itJarugHg   withDeuterated 

Acetylene  Assuming  Complete  Scrambling  and  a-secondary 
Isotope  Effects. 
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products  of  the  1-C3L3'^  +  C3H3I  reaction  are  of  higher  mass 
and  do  not  enter  further  Into  the  kinetic  schemes  modeled 
here  . 

Steady-state  conditions  are  applied  for  all  four  CgLg* 
Isotoplc  species  as  before  (without  regard  for  Isomeric 
differentiation)  and  the  differential  rate  equations  for 
scheme  shown  In  Figure  5.7  are  given  In  Appendix  II. B. 
Table  5.2  shows  fitting  results  for  several  data  sets  using 
the  same  kinetic  scheme.   Figure  5.8  shows  plots  of  the  best 
fit  of  this  model  to  the  same  data  set  as  fitted  In  Figure 
5.4. 

Some  experiments  were  conducted  where  certain  Ions  were 
ejected  from  the  analyzer  cell  using  FTICR  double  resonance 
techniques  (Comlsarow  et.  al.,  1978)  as  they  formed.   Among 
the  ions  ejected  were  C3H2D'^  and  C3HD2*.   Without  further 
fitting,  this  model  was  used  to  predict  the  behavior  of  the 
kinetic  system  If  these  ions  were  ejected.   Figure  5.9  shows 
predicted  results  and  data  points. 

Models  for  CgHg^  +  C^Hg.   The  reaction  of  C3H3*  with 
C4H2  (dlacetylene )  is  klnetlcally  more  complicated  than  that 
of  C3H3"*"  with  C2H2.   There  are  more  isomeric  possibilities, 
and  ion/molecule  reaction  products  of  higher  m/z  are 
detected.   Several  models  were  tried  in  attempts  to  fit  the 
experimental  data,  and  the  best  of  these  made  the 
assumptions  that  CjYl^*    and  the  excited  forms  of  CgHg**^  and 
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«^       ^11*^5^  ^^^    ^"  steady  state,  while  C5H3'^  and  CjHq'*^    are  not. 

'       In  order  to  account  for  the  build-up  of  CgHg"^  and  CjiHc"^ 
V:..      ions,  it  was  necessary  to  include  stabilization  steps  for 
the  excited  forms  of  these  Ions.   A  further  assumption  was 
made  that  the  stabilized  forms  of  CgHg"^  and  CjjHg*  are  not 
reactive  within  the  time  frame  (2s)  of  the  experiments.   The 
resulting  scheme  is  shown  in  Figure  5.10.   Reactions 
,^^     involving  CgHgl   are  included  because  CgHg*  is  detected  in 
It       the  ICR  experiments.   Applying  steady-state  assumptions  to 

the  C7H5*,  CgHg"^,  CijHg"^,  and  CgHg^  ions  yields  a  model  with 
f       a  full  analytical  solution  (Appendix  II. C).   Figure  5.11 

shows  results  of  the  model  fit  to  a  typical  data  set.   The 
.^i/;      kinetic  parameters  have  been  grouped  as  exponential  and 

pre-exponential  terms  in  the  equations,  yielding  9  fitting 
parameters.   Table  5.3  shows  results  of  model  fits  of  two 
data  sets.   The  term  k^gPc  H  j'^^ 

-(^2  +  k-2102  "^  '^sl^2^'  ^"*^  hence  can  be  calculated  from  the 
parameters.   The  results  of  the  first  data  set  fit  in  Table 
5.3  yield  k  2Pc„HoI  "  1.2(0.3)  s"^;  from  the  second  data  set 
fl^'  *^P2PC3H3I  =  1-3  (0.2)  s'^. 

Discussion 
Q^^l—±-£^K2l^R2  •       Several  insights  into  the  mechanism 
f.      of  the  C3H3*  +  C2H2/C2D2  reactions  are  obtained  from  these 
b       modeling  studies.   First,  a  value  of  1.4  (0.2)  s~^  is 
■Iv      obtained  for  kp '  from  the  intercept  of  Figure  5.2  (values 
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for  k  '  are  a  little  higher  from  model  fits  shown  In  Table 
5.2).   Assuming  a  cell  temperature  of  363K  and  given  that 
Pq  jj  J  =  3.3  X  10~®  torr  (after  system  factor  correction), 
k   =  1.6  (0.2)  X  10"^  cm^  s~^,  in  good  agreement  with  an 
experimentally  determined  value''.   Second,  the  model  gives 
branching  ratios  for  the  complex  CgLg"^.   The  ratio,  k^/kj, 
is  ca.  0.34  (.03)  (Table  5.2),  implying  that  CgLg"*^  fragments 
to  l-CgHg"*"  at  a  rate  three  times  that  of  fragmentation  to 
C-C3H3'*".   Since  C-C3H3'''  is  thermodynamical  ly  more  stable 
than  l-CgHg"^  by  ca.  25  kcal  mol~^  (Lossing,  1972),  any 
factor  affecting  the  internal  energy  of  CgL^"*",  such  as  a 
collisional  stabilization,  will  undoubtedly  affect  the  ratio 
kg/kj.   Hence,  k^/kj  is  probably  very  sensitive  to 
experimental  conditions  (method  of  ionization,  total 
pressure,  and  temperature).   Finally  the  "direction"  and 
magnitude  of  F  (1.67  (.10))  suggest  not  only  that  the  model 
is  reasonable,  but  that  the  isotope  effect  plays  an 
important  part  in  the  overall  kinetics. 

It  should  be  remembered  that  F  is  not  an  overall 
isotope  effect;  isotope  effects  are  ratios  of  rate 
constants.   The  magnitude  of  F  does,  however,  indicate  the 
kinetic  contribution  of  replacing  one  hydrogen  atom  with  a 
deuterium  atom  on  a  hydrogenic  site  undergoing  a 
hybridization  change  in  the  transition-state  structure.  F  is 
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"'k   =  (3.0  +  1.2)  X  10"^  cm^  mol"^  s~^  was  determined 
experimentally  using  a  typical  reaction  conditions  given  in 
Chapter  4 . 
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the  a-secondary  kinetic  Isotope  effect  for  the  elementary 
binolecular  addition  reactions  per  hydrogenlc  site  (only 
those  sites  giving  rise  to  a-secondary  effects  included). 
1/F  Is  the  Isotope  effect  for  the  elementary 
fragmentation  reactions  (to  I-C3H3''"  and  C-C3H3"*')  per 
hydrogenlc  site  (again  only  those  sites  giving  rise  to 
o-secondary  effects  Included).   Because  of  the  complexity  of 
the  scheme  In  Figure  5.7,  overall  Isotope  effects  cannot  be 
obtained  by  determining  product  ratios,  as  Is  often  done  In 
Isotoplc  studies  of  more  elementary  ion/molecule  reactions. 
However,  F  can  be  compared  to  Isotope  effects  found  In 
elementary  Ion/molecule  reactions. 

o-secondary  kinetic  Isotope  effects  can  be  quite  large 
(Mead,  et.al.,  1980;  Tumas ,  et.  al . ,  1987).   The  larger 
effects  found  In  Ion/molecule  reactions  are  due  to  a  narrow 
(non-  Boltzmann)  energy  distribution  centered  close  to  the 
threshold  energy  (Derrick,  1983).   The  magnitude  of  the 
Isotope  effect  factor,  F,  Indicates  a  fair  amount  of  bond 
rearrangement  In  the  transition  state  structures.  Due  to  the 
number  of  assumptions  and  the  indirect  method  for  obtaining 
F,  it  is  not  attempted  to  calculate  isotope  effects  from 
molecular  models  and  calculated  vibrational  frequencies. 

C3H3   +  C^H2-   Kinetic  data  and  modeling  of  this  system 
reveal  some  insight  into  the  mechanism  for  the  reaction  of 
C3H3"*'  and  C^Hg.   Although  the  scheme  shown  in  Figure  5.10 
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does  not  account  for  specific  isonerlc  reactivities,  the 
scheme  does  account  reasonably  well  for  the  data.   The 
signal-to-noise  ratios  of  the  available  data  do  not  warrant 
further  refinements  of  the  model.  Several  crude  models  were 
attempted,  and  the  model  reported  here  may  not  be  the  best 
or  only  model  which  can  adequately  explain  the  results. 
However,  of  the  variations  of  the  scheme  in  Figure  5.10 
which  were  tried,  the  one  reported  gave  the  best  fitting 
results.   In  particular,  a  scheme  was  tried  in  which  stable 
forms  of  CgHg"*"  and  CjjHg"^  were  formed  from  the  corresponding 
bimolecular  additions  reactions  rather  than  requiring 
stabilization  steps  of  "hot"  CgHg"^  and  CiiHg"^  ions.   This 
scheme  did  not  fit  the  data.  Apparently,  only  in  the  ions 
larger  than  C^Hg*  are  there  enough  internal  vibrational 
modes  to  give  rise  to  an  ion/molecule  collision  complex  of 
sufficient  lifetime  to  allow  colllsional  stabilization  at 
these  pressures. 

Interestingly,  the  value  for  kpg  in  the  scheme  of 
Figure  5.10  is  the  same,  within  error,  as  the  value  for  k^ 
in  the  scheme  of  Figure  5.7.   The  ratio  of  the  reduced 
masses  of  the  two  reacting  systems  (C3H3'^  +  C3H3I )  /  (C5H3*  + 
C3H3I)  Is  0.83,  which  Is  well  within  the  error  for  the  value 
of  kpg/kp.   Hence,  It  cannot  be  determined  from  this 
comparison  whether  the  reactions  with  C3H3I  are 
collision-controlled  or  not.   However,  if  the  reactions  of 
C3L3   with  C3H3I  and  C2L2  are  collision-controlled,  there 
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would  only  be  a  small  kinetic  isotopes  effect  arising  from 
the  rotational  partition  functions.   If  the  scheme  in  Figure 
5.7  is  correct,  there  is  a  transition-state  structure  for 
the  reaction  of  CgHg"^  and  CgHg ,  and  undoubtedly  with  C3H3'^ 
and  C^Hg  as  well,  which  does  have  bond  changes.   This 
implies  that  the  rate  is  predicted  to  be  lower  than 
collision-controlled. 

Figure  5.10  does  not  show  C5H5*  and  C^Hg*  fragmenting 
to  any  products.   Undoubtedly  any  isomers  of  05113"^  and  CrjH^'^ 
which  form  are  all  resonance-stabilized.   Any  fragmentation 
product  would  be  highly  energetic  and  unlikely  to  form,  even 
at  the  energy  levels  available  in  these  reacting  systems. 
Hence,  once  C^Hg"^  and  C,7H3'^  are  formed,  they  will  react  only 
in  bimolecular  addition  reactions,  which  are  orders  of 
magnitude  slower  than  the  unimolecular  fragmentation 
reactions  of  CrjH^'^ ,    CgH^"*^,  and  CjjHg^.  Hence,  there  is  a 
non-steady-state  build  up  of  C5H3'^  and  C7H3'^  ions  at  the 
first  part  of  the  reaction. 


CHAPTER  6 

REACTIONS  OF  CgHa""  AND  CgHg""  IONS  WITH  ACETYLENE  AND 

DIACETYLENE 


Introduction 
According  to  the  proposed  lon/nolecule  mechanism  of 
soot  formation  (Calcote.  1981;  Olson  and  Calcote,  1981a), 
C3H3   forms  CgHg"^  and  CgHg"^  in  reactions  with  acetylene,  and 
^1^5*   when  reacting  with  diacetylene.   However,  as  seen  in 
Chapter  4,  investigation  of  the  reactions  of  CgHg"^  with 
acetylene  and  diacetylene  did  not  reveal  facile  formation  of 
*'5^3*  °^    ^5^5*  '^y  reaction  of  this  ion  with  acetylene, 
although  CgHg"^  did  react  readily  with  diacetylene,  yielding 
C5H3   and  CyHg*.   All  three  of  the  postulated  products  of 
CgHg"^  reactions  with  acetylene  and  diacetylene  (CgHg"^, 
C5H5  ,  Cr^Hg  )  have  been  identified  by  mass  spectrometry 
(Michaud  et  al . ,  1981;  Olson  and  Calcote,  1981b)  in  flames, 
but  it  remains  to  be  determined  which,  if  any,  of  them  may 
be  important  in  soot  formation  mechanisms. 

Heats  of  formation  of  the  C5H3*  ion  from  different 
precursors  have  been  reported  in  three  experimental  studies 
(Franklin  and  Carroll,  1969;  Dannacher  et  al . ,  1979;  Baer  et 
al.,  1979).  Dannacher,  et  al  (1979)  suggested  the  presence 
of  two  different  linear  structures,  CH3-C-C-C=C*  and 
HC'C-C^C-CHg"^  for  the  C^H^*    ions  produced  by  .  CH3  and  H. 
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loss  from  2 , 4-hexadlyne  and  1 , 3-pentadlyne  parent  Ions, 
respectively. 

There  are  «any  possible  structures  for  the  C5H5*  Ion, 
and  despite  a  number  of  theoretical  and  experimental  studies 
involving  it,  few  definitive  results  exist  regarding  the 
specific  relative  energies  of  various  isomeric  forms. 
Experimental  studies  concentrated  on  determining  the 
appearance  potential  and  heats  of  formation  of  CgHg"*^  ions 
from  different  sources  by  mass  spectrometric  methods  (Pottle 
and  Lossing  1963;  Dorman,  1965;  Harrison  et  al.,  1965; 
Occolowltz  and  White.  1968;  Franklin  and  Carroll,  1969; 
Lossing  and  Traeger,  1975;  McCreary  and  Preiser,  1978). 
These  resulted  in  heat  of  formation  values  ranging  from  239 
to  309  kcal/mole  depending  on  the  source  and  technique  of 
CgHg"^  production.   Probably  the  best  estimate  of  heat  of 
formation  is  for  the  cyclopentadienyl  Ion  (AH^  =  255 
kcal/mole),  obtained  by  Lossing  and  Traeger  (1975)  from 
measurement  of  the  ionization  potential  of  the  ^ 

cyclopentadienyl  radical  with  an  electron  monochromator-mass 
spectrometer.   Other  than  this,  all  the  AH^'s  measured  from 
different  sources  using  electron  Impact  techniques  fall  into 
two  categories:   270-290  kcal/mole  and  300-309  kcal/mole. 
Evaluating  these  differences  in  terms  of  different  isomeric 
structures  cannot  be  Justified  considering  the  large 
uncertainty  in  the  internal  energies  of  ions  which  may  be 
present  following  electron  Impact  ionization.   Early  ICR 
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experiments  were  carried  out  to  Identify  structures  of  CgH^"^ 
according  to  their  reactivity  with  different  neutrals  (Brill 
and  Eyler,  1981;  Buckley,  1982;  Brill.  1983,  Eyler,  1984). 
Several  precursors  were  used  in  the  formation  of  C^Hg   ions 
and  the  results  of  the  reactivity  studies  indicated  the 
possibility  of  four  different  isomers  of  CgHg^  (Brill, 
1983).   Acetylene  reacted  quite  slowly  with  the  CgHg"*"  ions 
(Buckley,  1982;  Brill,  1983)  while  dlacetylene  and  aromatics 
with  side  chains  reacted  at  an  appreciable  rate  (Buckley. 
1982).   As  a  result,  it  has  been  suggested  that  soot 
nucleatlon  may  proceed  by  adding  a  few  large  molecules 
rather  than  through  addition  of  many  smaller  ones  (Buckley, 
1982;  Baykut  et  al . ,  1986).   No  definitive  assignment  of 
CgHg^  isomeric  structure  was  possible  in  the  earlier  studies 
(Buckley,  1982;  Brill,  1983).   Proton-transfer  reactions 
involving  one  relatively  unreactlve  CgHg"*"  Isomer  gave 
(Buckley,  1982)  a  proton  affinity  of  227.9  +  0.3  kcal/mole 
for  the  CgH^  neutral  which  remained  after  proton  transfer. 
When  combined  with  estimates  of  the  heats  of  formation  of 
possible  CgH^  species,  the  results  were  consistent  with  (but 
did  not  conclusively  prove)  a  vlnylcyclopropenyl lum  form  for 
the  CgHg"*"  isomer. 

A  number  of  theoretical  studies  have  examined  CgHg"^ 
structures  (Stohrer  and  Hoffman,  1972;  Hehre  and  Schleyer, 
1973;   Dewar  and  Haddon ,  1973;  Kollmar  et  al ,  1973;  Borden 
and  Haddon,  1979).  Schleyer  and  co-workers  located  two 
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minimum  energy  forms  on  the  C^Hg"^  potential  surface,  the 
more  stable  one  corresponding  to  planar  cyclopentadlenyl , 
the  other  one  to  a  square-based  pyramid  structure  (Kollmar 
et  al.  1973;  Hehre  and  Schleyer,  1973).   Similar  results 
were  reported  by  Stohrer  and  Hoffman  (1972)  although  they 
proposed  the  pyramidal  structure  as  the  more  stable  form 
compared  to  the  planar  cyclopentadlenyl.   Recent 
calculations®  have  shown  that  the  vlnylcyclopropenyl lum  Ion 
has  the  lowest  energy  (AH^  =  256.7  kcal/mole)  among  a  number 
of  other  possible  structures  as  shown  In  Figure  6.1.   All  of 
these  structures  were  found  to  be  no  more  than  30  kcal/mole 
higher  In  energy  than  the  lowest  energy  structure''. 

In  this  work,  the  kinetics  and  reaction  mechanisms  of 
CgHg   and  G5H5   Ions  produced  from  different  precursors  and 
reacting  with  acetylene  and  dlacetylene  have  been  studied  In 
order  to  Identify  Isomeric  structures  and  to  obtain  rate 
coefficients.   Also,  the  energetics  of  CgHg"^  formation  from 
norbornadlene  and  cycloheptatrlene  were  followed  to 
Investigate  two  possible  competing  pathways  of  CgHg"^ 
formation  which  could  lead  to  different  Isomeric  structures. 


f 


Experimental 
Reaction  rate  coefficients  were  determined  by 
monitoring  the  intensity  of  the  CgHg^  or  C-Hg"^  ions, 
respectively,  as  a  function  of  time  after  ejection  of  all 

"J.  Feng,  J.  Leszczynskl  and  M.  C.  Zerner,  submitted. 
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other  Ions  from  the  analyzer  cell.   Excitation  amplitude  and 
total  pressure  were  kept  constant  at  optimized  values  for 
all  the  kinetic  runs  In  order  to  minimize  the  unwanted 
effects  of  such  factors  on  ion  intensities,  which  have  been 
discussed  in  detail  elsewhere.^ 

C5H3'*'  ions  were  produced  by  50  eV  electron  ionization 
of  2 , 4-hexadiyne  and  by  the  reaction  of  CgHg"^  with 
diacetylene.   CgHg"*^  ions  used  to  produce  CgHg"^  were  formed 
by  Xe*  charge  transfer  reactions  with  propargyl  iodide  at  an 
ionizing  energy  of  15  eV.   After  30  ms  reaction  time  with 
C4H2.  all  the  other  ions  except  CgHg*  were  ejected  from  the 
cell  to  follow  the  reactions  of  this  Ion  with  diacetylene  as 
a  function  of  time.   CgHg"*"  ions  were  produced  by  charge 
transfer  reactions  of  various  precursors  ( dicyclopentadiene , 
cyclopentadiene ,  norbornadiene ,  l-penten-3-yne , 
cycloheptatrlene )  with  different  charge  transfer  agents 
(Xe"*",  Ng  .  Ar  )  formed  with  an  electron  beam  pulse  of  5  ms 
duration  at  an  ionizing  energy  of  13,  16.5,  and  20  eV, 
respectively. 

Cyclopentadiene  was  prepared  by  cracking 
dicyclopentadiene  and  was  kept  in  dry  ice  when  not  used  to 
prevent  the  dlmerization  process.   All  the  other  compounds 
used  were  obtained  commercially  and  their  purity  was  checked 
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M.  Noini  and  J.  R.  Eyler,  to  be  published 
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by  obtaining  wide  mass  range  spectra.   All  the  samples  were 
used  after  multiple  f reeze-pump-thaw  cycles. 

All  reactions  were  followed  at  a  cell  temperature  of 
363  K.   Some  C5H3'*'  and  all  65115''"  ions  used  in  rate  constant 
determination  studies  were  produced  by  chemical  reactions  in 
order  to  minimize  the  internal  energy  imparted  to  the  ions. 
Also,  since  the  total  pressure  in  the  reaction  cell  was 
almost  1  X  10"   torr  and  the  reactant  ion  formation  time  was 
>30  ms ,  substantial  collisional  relaxation  of  the  ions  took 
place  before  kinetic  data  were  collected. 


Results 
Reactions  of  C5H2"'".    A  very  low  number  of  C5H3"''  ions 
were  produced  from  2 , 4-hexadiyne  at  electron  energies  above 
30  eV.   Reactivity  of  these  ions  was  monitored  at  an 
electron  energy  of  50  eV  although  the  ion  Intensity  was 
still  very  low.   Although  almost  all  of  the  C5H3"*"  ions 
produced  reacted  with  the  2 , 4-hexadiyne  precursor,  no 
reaction  was  observed  with  either  C2H2  or  C4H2-   Reactions 
of  C5H3'''  with  2 , 4-hexadiyne  were: 


r^v' 


C5'*3   *  ^6^6 
^5^3   *  CgHg 

^5^3   *  ^6^6 


-->  CgHe*  +  [C5H3] 

■    y  CqH»7  +  l-^22-l 

-->  CgHg^  +  [C5H4] 

-->  C^H^^  +  [C4H2] 


(6.1) 
(6.2) 
(6.3) 
(6.4) 
(6.5) 
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CgHy    +   CgHg >   C25HJ3 


(6.6) 


CgH3'''  ions  produced  as  products  of  the  reaction  of 
C3H3''"  +  C4H2  were  100%  reactive  with  both  propargyl  Iodide 
and  with  dlacetylene  (the  neutrals  present  in  the  reaction 
medium) . 
C5H3'''  reactions  with  C4H2  were: 


C^Hg   +  C4H2 
^9^5   "^  ^4'^2 


>  CnHo    ■*"  LConpJ 

>  C9H5* 

-  — >  CllHg^ 

>  C23H7 


(6.7) 

(6.8) 

(6.9) 

(6.10) 


Some  of  the  product  Ions  were  observed  to  react  further  with 
propargyl  Iodide  by  displacing  atomic  iodine: 


CjUq      +  C3H3I      >  CjqHq   +  I 
C9H5   +  C3H3I      >  CjgHg   +  I 


(6.11) 
(6.12) 


CgHg"*"  reactions  with  propargyl  iodide  were 


^e^6       *    ^3^3^       ^  ^11^9  *    ^ 


(6.13) 
(6.14) 


The  CgHg"*"  ion  reacted  further  with  C4H2: 


109 


CgHg   +  C4H2      >  CjgHg 


(6.15) 


t:.i 


Ion  Intensity  vs.  time  curves  for  the  CgHg"*"/ ( C^Hg  + 
C3H3I)  system  are  shown  In  Figure  6.2.   The  decay  of  C5H3'*' 
Ions  Involved  reaction  with  both  C^Hg  and  propargyl  iodide. 

The  procedure  for  rate  coefficient  determination  used 
in  previous  studies  of  C3H3*  and  CgHg"^  ion/molecule 
reactions  required  subtraction  of  the  observed  rate  constant 
for  the  reaction  of  the  ion  with  the  precursor  neutral  from 
the  observed  rate  constant  for  the  sum  of  the  reactions  with 
precursor  and  reactant  neutrals.   In  this  case,  however,  the 
reactant  neutral  (C^Hg)  was  also  the  precursor  of  the  ion  of 
Interest  (C5H3'*'),  and  thus  the  subtraction  procedure  could 
not  be  used.   An  alternative  method  for  rate  constant 
determination  was  thus  required.   Following  ejection  of  all 
other  ions  from  the  FTICR  cell,  the  decay  of  C5H3'^  as  a 
function  of  time  is  given  by: 


[CsHg-]^  =  [C5H3Moe 


("p'^p  ■"  "d'^d)^ 


(6.16) 


where  Hpkp  and  n^jk^j  refer  to  the  products  of  the  number 
densities  and  ion/molecule  rate  coefficients  for  propargyl 
iodide  and  diacetylene.  respectively  (Reactions  (6.7),  (6.8) 
and  (6.13)).   The  quantity  n^jk^j  +  n  k   can  thus  be 
determined  from  the  slope  of  a  semilog  plot  of  05113"^  decay 
as  a  function  of  time.   At  short  reaction  times,  the 
following  expressions  hold  true: 
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C5H3+  +  C4H2  Reactions 
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C4H2. 


Disappearance  of 


Figure  6.2.   Reactions  of  CgHg"^  with 

C5H3   and  product  ions  includes  reactions  with  propargyl 
iodide.   C^Hg"^  ions  were  produced  from  the  reaction  of  CgHg" 
with  diacetylene  within  30  ms  reaction  time.   CgHg"*^  ions 


were  produced  from  propargyl  iodide  by  charge  transfer 


reactions  with  Xe"^ 


1.3  X  10  °  torr ,  p(Xe) 


P^CgHg 


I)  = 
4  X 


1  .9 


10 


-6 


X  10 

torr  . 


-7 


torr  , 


P(C4H2 


)  = 


^., 


Ill 


dD/dt  =  n^kd  [CgHg^loe 


-("pkp  +  n^kj)! 


and 


dP/dt  =  npkp  [C5H3^]oe 


-("pkp  +  njjkjj)t 


(6.17) 


(6.18) 


where  D  and  P  refer  to  product  Ions  of  the  reaction  with 
dlacetylene  (C^Hg"^  and  CgHg^,  Reactions  (6.7)  and  (6.8))  and 
with  propargyl  iodide  (CgHg"^,  Reaction  6.13)  respectively. 
Thus,  the  ratio  n^jk^j/npkp  was  calculated  from  the  ratio  of 
slopes  of  product  formation  as  a  function  of  time.   Next 
n^jkjj  was  obtained  using  the  calculated  sum  and  ratio  of  the 
two  rate  constants.   Finally  the  absolute  rate  constant  was 
determined  following  the  usual  procedure.   It  was  found  that 
CgHg"^  Ions  reacted  with  C4H2  with  a  rate  constant  of  (5.6  + 
1.7)  X  10"^°  cm^/s. 

Kinetic  Modeling  of  CgHg'*'  reactions  with  dlacetylene. 
In  order  to  better  understand  the  reaction  mechanisms 
Involved  in  the  CgHg"*^  +  C^Hg  system,  kinetic  modeling 
studies  of  the  ion  intensity  vs.  time  curves  shown  in  Figure 
6.2  were  carried  out.^°   Details  of  the  modeling  methodology 
are   given  in  Chapter  5.    A  kinetic  model  involving  the 


^^The  kinetic  modeling  studies  were  performed  in  the 
Envlronlcs  Division  of  Air  Force  Engineering  and  Services 
Center,  Tyndall  Air  Force  Base,  Florida  by  F.  Wiseman  using 
multiple  experimental  data  sets  produced  at  identical 
conditions  to  those  reported  for  Figure  6.2. 
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reactions  above  was  fitted  to  the  experlfflental  data  assuming 
steady-state  concentrations  for  the  excited  forms  of 
(^g'^S  )  ^^^    ^^11^5*)   Intermediate  complexes.   The  complete 
reaction  scheme  used  for  modeling  Is  shown  In  Figure  6.3. 

Comparison  of  the  model  fit  with  the  experimental  data 
is  shown  in  Figures  6.4a  and  6.4b. 


Effect  of  Precursor  Neutrals.   CgHg"^  ions  produced  from 
different  precursors  exhibited  behavior  indicative  of  both 
reactive  and  unreactive  populations  toward  both  the 
precursors  and  the  reactant  neutrals.   Figures  6.5a  and  6.5b 
show  C5H5   normalized  ion  intensity  vs.  time  curves  for  two 
(cycloheptatrlene  and  l-penten-3-yne )  of  the  five  precursors 
used  to  form  ions  before  reaction  with  dlacetylene  and 
acetylene.   Following  an  exponential  decay  Indicative  of 
pseudo-first  order  kinetics,  a  substantial  fraction  of 
unreactive  ions  remains  at  long  reaction  times,  particularly 
in  Figure  6.5b.   Table  6.1  shows  the  method  of  preparation 
of  C5H5   ions  from  different  precursors  and  the  percent  of 
unreactive  ions  remaining  at  long  reaction  times. 

Two  (norbornadlene  and  cycloheptatrlene)  of  the  five 
precursors  mentioned  above  produce  CjKq^    rather  than  CgHg"^ 
parent  ions.   The  C^Hg"^  ions  produced  from  norbornadlene 
were  reported  (Davidson  and  Skell,  1973;  Ausloos,  1982) 
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^9^5*  "^  ^4^2      ^  ^13^7 


f3 


CvHs"  *  C4H2  >  C11H5+ 


P3 


CtHs"  +  C3H3I >  CioHei"- 


-32 


(C11H5+)*  -=-->  C7H3*  +  C4H2 


-31 


(CjiHg*)*   --  — >   CgH3*   +   C2H2 


's2 


'  (CiiHg-^)* >  C11H5*  (stabilized) 


Figure  6.3.   Reaction  Scheme  Postulated  for  the  Kinetic 
Modeling  of  the  Reaction  of  CgHg"^  with  Dlacetylene. 
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Reactions.^  (a)  Ion  Intensity  vs  time  curves  for  CgHg"^  and 
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The  best  fit  values  for  the  fitting  parameters  were 
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TABLE  6.1 

'•■*■ 

Percentages^  of  unreactlve  CgH^"^  found  from  various 
precursors  by  charge  transfer  chemical  ionization  monitored 
by  observing  reaction  with  the  precursor  neutral. 


Precursor 
Neutral 

Charge  Transfer           *  unreactlv 
Agent                     CgHg^ 

dicyclopentadiene 

Ng^/Ar""                      20 

cyclopentadiene 

Xe""                        17 

l-penten-3-yne 

Xe"^                       65-70 

cycloheptatriene 

Ar*                        0-5 

norbornadiene 

Ng^/Ar*                     18-20 
Kr"^                        5-10 

.  r 

■■■  ,  .   s 


a 


Estimated  error  is  5%. 
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earlier  to  lead  to  two  different  channels  (Reactions  (6.19) 
and  (6.20))  for  CgHg^  Ion  formation. 


or 


C^Hg 


^5^6^ 


>  C5H5   +  CgHg 


>   CgHg     +   CgHg 

>    C5H5*  ^  H 


(6.19a) 
(6.19b) 

(6.20a) 
(6.20b) 


Thus,  CgHg"*"  Ion  formation  pathways  were  studied  In  more 
detail  for  norbornadlene  and  cycloheptatrlene .   When 
different  charge  transfer  gases  were  used  for  Ionization  of 
norbornadlene,  different  percentages  of  reactive  CgHg*  were 
observed,  as  Is  noted  In  Table  6.1.   Also  the  abundances  of 
CyHfy"*" ,  CgHg"*"  and  CgHg"*"  Ions  were  measured  following  charge 
transfer  Ionization,  and  different  behavior  was  observed  for 
compounds  with  Ionization  potentials  In  the  range  of  14  to 
16  eV.   Table  6.2  shows  this  effect  for  two  reagent  gases 
(Kr*  and  Ng"^)  which  have  Ionization  energies  of  14.0  and 
15.7  eV  respectively.  ■• 

In  order  to  further  investigate  this  behavior,  the 
relative  ion  abundances  vs.  electron  energy  for 
cycloheptatrlene  and  norbornadlene  were  obtained  (see 


Figures  6.6  and  6.7) 


11 


^^Electron  Impact  ionization  was  used  for  this  study  and 
therefore  the  energy  scale  in  Figures  6.6  and  6.7  should  be 
considered  as  approximate  with  at  least  +1  eV  uncertainty. 
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TABLE  6.2 

Changes  in  ion  abundances  at  two  different  ionization 
energies  for  norbornadiene . 


Ionizing  energy/eV 
(charge  transfer  agent 

Abundances  as  a  fraction  of 
)           total  ion  signal^ 

C^Hq         C7H7         CgHg         C5H5 

14.0  (Kr*) 

0.07      0.45       0.32      0.16 

15.7  (Ng-^) 

0.09      0.28       0.16      0.47 

a 


Estimated  error  is  +0.02 
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Figure  6.6.   Relative  Ion  Intensities  Produced  from 
Cycloheptatriene  as  a  Function  of  Electron  Impact  Energy 
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Figure  6.7.   Relative  Ion  Intensities  Produced  from 
Norbornadiene  as  a  Function  of  Electron  Impact  Ener 
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CgHg"*"  reactions  with  dlacetylene.   After  ejection  of 
all  Ions  except  CgHg"*"  following  charge  transfer  chemical 
Ionization  of  a  mixture  of  dlacetylene  and  a  precursor 
compound,  the  Ion/molecule  reactions  as  a  function  of  time 
were  monitored.   Independent  of  the  precursor  used,  the  main 
reaction  was  the  addition  of  C4H2  to  C^Hg"^  to  produce  the 
CgHy"^  Ion. 


^5^5   *  ^4*^2      ^  ^9^7 


(6.21) 


Other  minor  reactions  observed  were: 


^5^5   "^  ^4*^2 
^5^5   *  ^4^2 

9  7   ^  ^4" 2 


-->   C^Hg^   +   [C2H2] 

— ->  CyH^"^  +  [Cg] 


->   CjjHg"^ 


>  C,oH 


13"9 


(6.22) 
(6.23) 
(6.24) 
(6.25) 


For  CgHg"^  Ions  produced  from  cyclopentadlene ,  the  C, 
addition  reaction  (6.22)  was  not  observed  with  C^Ho. 

The  rate  coefficient  for  the  disappearance  of  CgHc^ 
was  calculated  as  described  earlier  In  Chapter  3  and  the 
values  found  for  CgHg"*^  Ions  produced  from  different 
precursors  are  given  In  Table  6.3. 


CgHg"^  reactions  with  acetylene.   CgHg"*"  Ions  formed  from 
four  of  the  five  precursors  (cyclopentadlene. 
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TABLE  6.3 


Rate  coefficients  for  the  reaction  of  CgHc"*^  Ions  from 
different  precursors®  with  dlacetylene  and  acetylene. 


Precursor  Absolute  rate  coefficient  with 

C^Hg/ClO"^^  cm^/s)      CgHg/ClO"^^  cm^/s) 

cyclopentadiene 

l-penten-3-yne 

norbornadiene 

cycloheptatrlene 

dicyclopentadlene*' 


(1.0  +  0.5) 

(2.0  +  1.0) 

(2.9  +  1.4) 

(3.3  +  1.9) 

(1.6  +  0.8) 


(4.8  +  1.9) 

(2.7  +  1.0) 

(3.1  +  1.3) 

(1.8  +  0.7) 


^CgHg"*"  ions  were  produced  by  charge  transfer  chemical 
ionization  using  different  gases  as  shown  in  Table  6.1. 
"Ion  signal  for  dicyclopentadiene  was  too  small  to  produce 
reproducible  results  for  the  rate  coefficient  of  the 
reaction  with  acetylene. 


i;.*,< 


u 
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l-penten-3-yne ,  norbornadlene ,  cycloheptatrlene )  reacted 
very  slowly  with  C2H2,  producing  very  small  amounts  of  C^Hi^"*" 
and  CgHg"^. 

CgHg"^  +  CgHg >  C^H^""  (6.26) 

CjH^*    +  CgHg  >  C9H9+  (6.27) 


Rate  constants  for  the  CgHg"*^  +  CgHg  reaction  were  about  one 
order  of  magnitude  less  than  those  for  reactions  with  C^Hg 


(see  Table  6.3). 


Discussion 

CgHg^.   Two  different  CgHg"*^  structures  were  postulated 
(Dannacher  et  al.,  1979)  based  on  PIPECO  measurements  of 
appearance  energies  of  C5H3*  Ions  from  1 , 3-pentadlyne  and 
2 ,4-hexadlyne.   In  the  work  reported  here,  CgHg"^  Ions  were 
produced  from  2 , 4-hexadlyne  and  by  reaction  of  C3H3*  with 
dlacetylene.   Although  CgHg"*"  Ions  produced  either  way  were 
reactive  toward  their  precursor  neutrals,  CgHg"^  Ions  from 
2, 4-hexadlyne  were  unreactlve  with  both  C^Hg  and  CgHg.  while 
those  formed  as  Ion/molecule  reaction  products  reacted  with 
C4H2  with  the  relatively  high  rate  constant  of  (5.6  +  1.7)  x 
10"^°  cm^/s.   Dannacher  et  al .  (1979)  suggested  the 
structures  CH3-C  =  C-C  =  C'^  and  HC  =  C-C  =  C-CH2'^  for  the  C^H^*    Ions 
produced  from  2 , 4-hexadlyne  and  1 , 3-pentadlyne  respectively. 
Experimental  appearance  potentials  for  these  Ions  suggested 


'■'■'i 

'■■■■■  i 


<«<- 


I-.;, 

3^' 
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that  the  former  probably  had  a  higher  heat  of  fornatlon, 

<1431  kJ/nole  (Baer  et  al . ,  1979),  compared  to  the  latter, 

1317  kJ/mole  (Dannacher  et  al . ,  1979).   In  this  study, 

i 
enhanced  reactivity  of  C5H3*  ions  produced  by  the  C3H3*  + 

C4H2  reaction  compared  to  those  produced  by  electron 

ionization  of  2 , 4-hexadlyne  suggests  that  the  former  have  a 

higher  heat  of  formation  with  a  stable  structure  different 

from  those  reported  earlier  for  C5H3'*'  ions  from  both 

2 , 4-hexadiyne  and  1 , 3-pentadiyne .   Enhanced  reactivity  does 

not  always  correlate  with  a  higher  heat  of  formation,  but 

this  has  often  been  found  true  for  isomers  of  other  small 

hydrocarbon  ions.   Our  observation  that  CgHg*  is  formed  in 

an  exothermic  reaction  of  linear  C3H3'''  with  diacetylene 

leads  to  an  upper  limit  of  1375.9  kJ/mole  for  AH^  of  the 

C5H3*  ion  formed  in  this  manner. ^^   This  value  is  not 

Inconsistent  with  either  of  those  reported  earlier. 

Although  no  structural  assignment  can  be  made  for  the 

C5H3*  ion  produced  by  reacting  CgHg"^  with  diacetylene  based 

on  the  results  of  this  study,  several  structures  can  be 

considered  possible  on  the  basis  of  initial  charge 

distribution  on  C3H3'^  leading  to  the  formation  of  three 

product  ion  structures  by  reaction  mechanisms  shown  in 

Figure  6.8.   The  three  possible  C5H3*  product  ions  shown  in 


12r 


'This  calculation  assumes  no  significant  internal 
energy  in  the  C3H3'*"  ions  prior  to  reaction  and  heats  of 
formation  of  281  kcal  mol"^  (Losslng,  1972),  for  I-C3H3*. 
54.2  kcal  mol""^  (Wagman  et  al  .  ,  1968),  and  102  kcal  mol~* 
for  C4H2  (Coats  and  Anderson,  1957). 


■:m 
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C  =  C  =  C-H 


H-C=C-C=C-H 


>   ^C  =  C  =  C'^   + 

W^  ^  =  C-C  =  C-H 


^C  =  C  =  C-^     C  =  C-H 


> 


> 


c=c  =  c 


^C=C-H 


H-C=C-H 


(6.10) 


Hv 
H^ 


C-C=C-H   +   H-C=C-C=C-H 


■>       H-C=C-C^  ^C^ 


"C  =  C-H 


> 


H-C  =  C-C' 


H 


■H 


C< 


~C  =  C-H 


> 


H-C  =  C'^  ^C  =  C-H    +    H-C  =  C-H 


and/or 


(6.11) 


H-C=C-C^I 


./ 


-c+ 


■4-f^ 


H-C=C-H 


(6.12) 


Figure  6.8.   Proposed  Reaction  Mechanisms  for  the  Reaction 
of  the  Propargyllum  Cation  with  Diacetylene,  Leading  to 
Three  Different  C5H3*  Product  Ions  with  Resultant  Loss  of 
Acetylene . 
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Figure  6.8  have  structures  which  differ  from  those 
postulated  (Dannacher  et  al , ,  1979)  earlier.   While  no 
definitive  theoretical  or  experimental  evidence  as  to  the 
relative  stability  of  these  (or  other)  CgHg"*"  Isomers  Is 
available.  Structure  (6.12)  Is  most  similar  to  the 
vlnylcyclopropenyllum  ion  found  most  stable  among  the  CcHr"*^ 
Isomers.   Thus  Structure  (6.10)  or  (6.11)  might  be  the 
unstable,  reactive  structure  formed  in  these  experiments. 
When  deuterated  dlacetylene  was  reacted  with  CgHg^  ion,  the 
deuterated  product  ratios.  CgHgD"^ :  CgHDg"^ :  CgHg*  were  found  to 
be  6:3:1.   This  result  is  also  consistent  with  the  reaction 
mechanisms  of  formation  of  the  three  C5H3*  structures  shown 
in  Figure  6.8  assuming  complete  scrambling  of  hydrogens  in 
the  intermediate  CyHg"^  complex. 


C^H^   ion  structures  produced  from  different 
precursors .   CgHg"^  ions  produced  from  different  precursors 
all  exhibited  at  least  two  populations,  one  reactive,  other 
unreactive.   This  behavior  for  numerous  other  ionic 
reactants  has  been  used  previously  (Ausloos  and  Lias,  1981; 
Smyth  et  al . ,  1982;  Wagner-Redeker  et  al . .  1983;  Smith  and 
Adams.  1987)  to  argue  for  the  existence  of  at  least  two 
different  structural  Isomers;  one  reactive  and  one 
non-reactive.   These  results  are  in  agreement  with  the 
earlier  reactivity  study  of  C^H^*    ions  from  different 
precursors  (Brill,  1983).   To  distinguish  the  reactive  C^H^* 
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structures,  their  reaction  mechanisms  and  reaction  rates 
with  acetylene  and  dlacetylene  were  studied.   Although  they 
reacted  with  their  precursor  neutrals  with  different  rates, 
they  all  reacted  with  acetylene  and  dlacetylene  at  similar 
rates  within  the  experimental  error  limits,  as  shown  in 
Table  6.3.   Thus  no  isomeric  differentiation  based  on 
reactivity  can  be  made.   Rate  coefficients  for  the  reaction 
with  acetylene  were  similar  to  those  determined  earlier 
(Buckley,  1982;  Brill,  1983),  i.e.  in  the  range  of  IQ-ll 
cm  /s.   Reaction  mechanisms  were  similar  for  all  CgHg*  ions 
with  the  exception  of  cyclopentadiene ,  which  didn't  give  C, 
addition  reactions  with  C^H,. 

Although  no  definitive  structure  assignment  can  be  made 
for  the  CgHg"^  ions  produced,  several  suggestions  are 
possible  from  the  results  of  this  study.   Similar  reaction 
rates  within  experimental  error  suggest  the  formation  of  the 
same  reactive  CgHg"^  structure  from  all  precursors.   Assuming 
that  the  latest  theoretical  Information^^  on  the  energies  of 
different  CgHg"^  structures  is  fairly  accurate,  we  can 
identify  the  possible  structures  which  are  energetically 
accessible  by  the  initial  internal  energy  transferred  to 
neutral  parents  by  chemical  ionization.   As  shown  in  Table 
6.4  ,  simple  calculation  using  the  thermochemical  data^^ 


13 


14 


J.  Peng.  J.  Leszczynskl  and  M.  C.  Zerner.  submitted. 


kca 
kca 


This  calculation  assumed  heats  of  formation  of  31.8 
1  »ol  for  cyclopentadiene  (Harrison  et  al . .  1960),  43 
1  mol    for  cycloheptatriene  (Finke  et  al . .  1956),  and 
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TABLE  6.4 

CoBparlson  of  the  energy  transferred  to  neutral  precursors 
by  chemical  ionization  with  theoretical  appearance  energies 
(assuming  no  reverse  activation  energy)  of  different  C5H5'*' 
structures  from  the  corresponding  precursors. 


Neutral 
precursor 


charge  transfer 
agent  and  (IP/eV) 


structure* 


theoretical 
AP/eV 


Cyclopentadiene 


Xe"^  (12.1) 


6.  1 
6.2&6 
6.4 
6 
6 
6 
6 


5 
6 
7 
8 


6.9 


12 

12 

12 

12 

12, 

12 

13, 

13 


Norbornadiene 


Kr* 

Ar* 


(14.0) 
(15.8) 


6.1 
6.2&6 
6.4 
6.5 
6.6 
6.7 
6.8 
6.9 


13 
13 
13 
13 
13 
14 
14 
14 


Cycloheptatriene 


N 


2^    (15.6) 


6.  1 
6.2&6 
6.4 
6 
6 
6 
6 
6 


13 
14 
14 
14 
14 
14 
14 
15 


Numbers  refer  to  the  structures  shown  in  Figure  6.1 


62.3  kcal  mol  ^  for  norbornadiene  (Lifshitz  and  Bauer,  1963) 
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available  for  cyclopentadlene ,  norbornadlene  and 
cycloheptatrlene  indicates  that  CgHg"*^  formed  from 
cyclopentadlene  by  Xe"^  charge  transfer  chemical  ionization 
can  only  have  the  lowest  energy  structure  (Structure  6.1). 
On  the  other  hand,  CgHg"^  ions  produced  from  norbornadlene 
and  cycloheptatrlene  by  N2"*^/Ar^  charge  transfer  chemical 
ionization  have  enough  excess  energy  to  sample  almost  all 
the  Isomeric  structures  shown  in  Figure  6.1.   Since  the 
energy  barriers  between  these  structures  are  not  yet  known, 
nothing  can  be  said  concerning  the  possibility  of 
isomerizatlon  between  these  structures.   Experimental  data 
indicates  that  at  least  two  stable  (reactive  and  unreactive) 
structures  are  formed  in  all  the  cases  studied.   More 
definitive  identification  of  CgHg"*^  structures  may  be 
forthcoming  when  better  theoretical  calculations  are 
reported,  although  identification  of  the 
vlnylcyclopropenylium  ion  as  (one  of)  the  unreactive 
structure(s)  seems  reasonable  in  the  light  of  the  recent 
calculations  and  earlier  ICR  studies. 

To  compare  the  behavior  of  ions  produced  from 
cycloheptatrlene  and  norbornadlene  precursors  (both  having 
the  formula  C^Hg),  relative  abundances  as  a  function  of 
electron  energy  shown  in  Figures  6.6  and  6.7  were  obtained. 
As  shown,  the  curves  for  CyH,^*  (Figures  6.6a  and  6.7a)  and 
C5H5*  and  03113"^  (Figures  6.6b  and  6.7b)  are  very  similar 
both  in  terms  of  energetics  and  of  general  shape. 


p--y- 
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Differences  In  CgHg"*"  relative  Intensities  are  probably  due 
to  the  normalization  procedure  used  in  each  case.   An 
increase  in  the  relative  intensity  of  C5H5'*'  ions  in  both 
cases  in  the  energy  range  of  15-25  eV  along  with  the  fact 
that  the  C^E.j'*'    curves  no  longer  Increase  in  this  energy 
range  implies  that  most  of  the  C5H5'''  ions  are  produced  from 
CjHrj'*' .      Similar  behavior  is  seen  for  €3113'^  ions,  i.e..  the 
intensity  of  these  ions  increases  at  about  22  eV  (on  the 
energy  scale  shown  on  Figures  6.6b  and  6.7b)  where  the  rate 
of  increase  in  the  intensity  of  CgHg"^  is  decreased.   This 
similar  behavior  implies  that  C^Hi^*^  ions  formed  from 
cycloheptatr iene  and  norbornadiene  follow  the  same       ■  >;  ^ 
fragmentation  pathways  in  the  same  energy  range.   Similar 
behavior  for  C^H.^'''  ions  from  these  two  compounds  has  also 
been  reported  earlier  in  collisional  activation  studies 
(Winkler  and  McLafferty.  1973). 

As  seen  in  Figures  6.6b  and  6.7b,  formation  of  CgHg"^, 
both  in  terms  of  energetics  and  relative  intensity,  is 
different  for  cycloheptatriene  and  norbornadiene,  indicating 
that  Ci^Hg"*"  ions  formed  from  these  compounds  by  electron 
impact  have  different  structures  (according  to  Reaction 
(6.20a).   Based  on  the  significantly  higher  relative 
intensity  of  CgHg"^  for  norbornadiene  along  with  the  results 
shown  in  Tables  6.1  and  6.2,  there  exists  the  possibility  of 
formation  of  CgHg"^  ions  from  the  CgHg"*^  channel  by  Reaction 
(6.20b)  for  norbornadiene.   As  shown  in  Table  6.1,  the 
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percent  of  unreactlve  CgHg"*^  increases  from  (5-10)%  to  (18- 
20)X  with  changing  ionizing  charge  transfer  energy  for 
norbornadiene  while  no  such  change  in  CgHc'*'  reactivity  was 
observed  for  cycloheptatriene  under  the  same  conditions. 
Relative  ion  intensities  at  two  different  energies  of  charge 
transfer  chemical  ionization  are  shown  in  Table  6.2  for 
norbornadiene.   The  results  indicate  that  when  Ar*/N2'*" 
replaces  Kr"*^  as  the  charge  donor,  the  relative  intensity  of 
the  observed  CgHg"*"  increases  while  there  is  a  corresponding 
decrease  in  the  yield  of  the  CgHg"*"  ion.   Thus,  evaluation  of 
the  results  on  CgHg^  reactivity  (Table  6.1)  together  with 
the  results  on  relative  ion  intensity  (Table  6.2)  implies 
that  CgHg*  formed  from  CgHg*  channel  has  a  higher  percent  of 
unreactlve  isomer  compared  to  that  formed  from  C-H,,''' 
channel . 
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CHAPTER  7 
REACTIONS  OF  GASEOUS  C^H^*  IONS 


.  _.  ■:         Introduction 
The  CyHy   Ion  has  been  proposed  as  one  species 
involved  in  an  ion/molecule  mechanism  of  soot  formation 
(Olson  and  Calcote,  1981a;  Calcote,  1981).   Mass 
spectrometrlc  measurements  have  shown  that  as  premlxed 
flames  approach  a  sooting  condition,  the  predominant  C^jHo"^ 
ion  is  replaced  by  larger  ions  (Olson  and  Calcote,  1981b) 
Including  C^H^^.   These  ions  in  turn  have  been  postulated 
to  react  with  neutrals  such  as  acetylene  and  dlacetylene, 
forming  still  larger  ions  (Olson  and  Calcote,  1981a). 
According  to  the  proposed  ion/molecule  mechanism.  CrHr"^ 
ferns  C,jH,j^    in  reactions  with  acetylene.   However,  the 
reactions  of  CgHg-^  with  acetylene  did  not  reveal  significant 
formation  of  C^H^"^  (see  Chapter  6).   Much  higher 
concentrations  of  C^H^^  than  observed  experimentally  were 
predicted  (Olson  and  Calcote,  1981a)  by  a  simple  model  which 
assumed  a  benzyl  structure  for  the  Ion.   To  account  for  this 
difference  between  experiment  and  theory,  it  was  suggested 
that  the  actual  ion  in  the  flame,  presumably  formed  by 
Ion/molecule  reactions  of  CgH^^  and  C^H^^    ions,  is  some 
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Other  CjHj'*'    Isoffler  less  stable  than  benzyl  (Olson  and 
Calcote,  1981a).       | 

There  have  been  a  significant  number  of  both 
theoretical  and   experimental  studies  on  structural 
characterization  of  CyH^"^   Ions.   Theoretical  studies 
(Abbaund  et  al . ,  1976;  Cone  et  al..l977)  on  the 
determination  of  heats  of  formation  for  C^H^*  Ions  Indicated 
that  tropyllum,  with  AHf=207.9  kcal/mol  (Abbaund  et  al . , 
1976),  and  195.6  kcal/mol  (Cone  et  al.,1977)  and  benzyl, 
with  AHf=217.1  kcal/mol  (Abbaund  et  al . .  1976),  and  220,4 
kcal/mol  (Cone  et  al.,1977)  are  the  most  stable  structures 
with  an  activation  energy  of  32.7  kcal/mol  (Cone  et  al . , 

1977)  for  benzyl  >  tropyllum  Isomerlzatlon .   Experimental 

measurements  of  ionization  potentials  for  tropyllum  and 
benzyl  radicals  by  monoenergetlc  electron  Impact  gave 
calculated  values  of  AHf=209  and  211  kcal/mol  (Thrush  and 
Zwolenlk,  1963;  Elder  and  Parr,  1969;  Losslng,  1971) 
respectively.   Studies  of  metastable  Ions  (Cooks  et  al  .  , 
1973;  Mclafferty  and  Winkler.  1974;  Grotemeyer  and 
Grueltzmacher,  1981)  colllslonal  activation  (Winkler  and 
McLafferty.  1973;  Koeppel  et  al . ,  1978McLaf f erty  and 
Bockhoff,  1979).  gas  phase  radlolysls  (Yamamoto  et  al . , 
1969;  Takamuku  et  al . ,  1972;  Takamuku  et  al . ,  1973;  Sagl  et 
al . ,  1974),  photolonlzatlon/photodlssoclatlon  (  Dunbar, 
1975;  Traeger  and  McLoughlln.  1977;  McCreary  and  Preiser, 
1978:  Traeger  and  McLoughlln,  1978;  McLoughlln  et  al . ,  1979; 
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Yaroslavtsev  et  al.,  1984).  PEPICO  (Boabach  et  al . ,  1983a, 

and  b),  and  Ion/molecule  reactions  (Shen  et  al,,  1974; 

Jackson.  1977;  Jackson  et  al . ,  1977;  Ausloos  et  al . ,  1980; 

Sharma  and  Kebarle.  1981)  of  C^H^"^  ions  from  different 

precursors  were  used  to  characterize  the  different 

structures.   General  experimental  evidence  including  the 

determination  of  appearance  potentials  and  heats   of 

formation  (Traeger  and  McLoughlin.  1977;  Traeger  and 

McLoughlin.  1978;  McLoughlin  et  al..  1979;  Bombach  et  al., 

1983a. and  b)  supported  that  these  ions  exist  in  the  gas   ,,. 

phase  with  at  least  two  structures;  tropylium.  and  benzyl. 

Unreactive  and  reactive  C^H^"^  observed  in  ion/molecule 

i 
reactions  have  been  attributed  to  the  tropylium  and  benzyl 

structures,  respectively  (Shen  et  al.,  1974;  Jackson,  1977; 

Jackson  et  al.,  1977;  Ausloos  et  al.,  1980;  Sharma  and 

Kebarle.  1981).   Furthermore,  it  was  proposed  that  there  is 

an  equilibrium  between  these  two  structures  (McLafferty  and 

Bockhoff,  1979;  Andrews  and  Keelan,  1981)  similar  to  that 

proposed  between  the  parent  ions  toluene  and 

cycloheptatriene  (Mclafferty  and  Bockhoff,  1979;  Ausloos, 

1982)  (see  Figure  7.1).    Several  groups  studied  the  effect 

of  the  energy  above  CyH,^*  formation  threshold  on  the 

dominance  of  one  structure  over  the  other  (McLafferty  and 

Winkler,  1974;  Shen  and  Dunbar,  1974;  Dunbar,  1975;  Ausloos. 

1982).   Their  results  in  general  were  supportive  of  the 

equilibria  between  different  structures  mentioned  above. 
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Toluene    ion 
A  Hf  =  215.4  kcal/mol 


0.9  ev 


0.5  ev 


-H' 


Benzyl 

AH,  =  2IlkcaI/moI 


1.8  ev 


2.0. 


■^ 


1.0  ev 


Cycloheptafriene 


ion 


AHj=    235.1     kcol/mol 


Tropylium 
AHf  =  209    kcal/mol 


Figure  7.1.   Proposed  (Mclafferty  and  Bockhoff   1979. 
Ausloos.  1982)  Equllbrlums  between  Different  C*Ho^  and  C  H  + 
Structures.  ^  8    "°  ^7^7 
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Although   McLafferty  et  al  (Winkler  and  McLafferty,  1973; 
McLafferty  and  Winkler,  1974)  reported  evidence  for  stable 
tolyl  and   norbornadlenyl  structures  In  their  earlier 
collisional  activation  studies  no  definitive  indication  for 
stable  structures  other  than  benzyl  and  tropylium  was 
provided  by  later  studies  using  different  precursors 
(Jackson  et  al.;  1977;  Koeppel  et  al.,  1978;  McLafferty  and 
Bockhoff,  1979;  Ausloos,  1982).   In  the  work  reported  here 
CrjHj      ions  have  been  produced  from  three  different         .-- 
precursors  using  different  charge  transfer  gases.   Reaction 
pathways  and  the  rate  coefficients  for  the  reaction  of  CrjH,,* 
ions   with  precursor  neutrals,  acetylene  and  diacetylene 
near  room   temperature  were  investigated. 

Experimental 
CyH-^   ions  were  produced  from  toluene,  cycloheptatriene 
and  norbornadiene  by  dissociative  charge  transfer  . 
Percentages  of  reactive/unreacti ve  isomers  were  determined 
when  the  three  neutrals  above  reacted  with  different  charge 
transfer  agents  (C4H2''.  Xe^,  Kr"^  and  Ar"")  formed  with  an 
electron  beam  pulse  of  5ms  duration  at  an  ionizing  energy  of 
11.5,  13,  15,  and  20eV,  respectively.   To  do  this,  all  ions 
except  CrjHrj'^    were  ejected  from  the  cell  and  enough  reaction 
time  (typically  3s)  at  the  pressures  used  (typically  10"^- 


10 


torr)  was  allowed  so  that  the  reactive  structures  could 
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be  lost  via  Ion/molecule  reactions,  thus  leaving  an 
unreactlve  population. | 

However,  In  the  case  of  cycloheptatrlene ,  as  reported 
earlier  (Jackson  et  al . ,  1977;  Ausloos,  1982),  because  of 
the  growth  of  unreactlve  Cp^H^"^  by  time  as  a  result  of  the 
ion/molecule  reactions  (7.1a)  and  (7.1b),  this  approach 
required  some   adjustments  In  the  experimental  sequence  in 
order  to  obtain  accurate   percentages  of  Ci^Hy"*"  populations. 


CfrHfT     +   CfrH 


'7"7 


7"8 


C^Hg     +   CiyHg 


>  C^Hg""  +  C7H7  (7.1a) 

>  CrjUrj*    (unreactlve)  +  C7H9     (7.1b) 


Thus,  at  long  reaction  times  C^Hg"^  was  ejected   contlnously 
during  reaction  time  In  order  to  avoid  reaction  (7.1b). 

All  reactions  were  followed  at  a  cell  temperature  of 
363K.   In  order  to  avoid  any  possible  effect  of  differences 
In  Initial  internal  energy  distribution  of  C^H^*  ions  on 
rate  coefficient   determination,  charge  transfer  agents 
which  produce  CrjE^'*'    ions   with  similar  internal  energies 
from  different  precursors  were   used  for  kinetic  runs. 
Also,  since  the  total  pressure  in  the   reaction  cell  was 
almost  1  X  10"^  torr  and  the  charge  transfer  time   was 
100ms,  substantial  colllsional  relaxation  of  the  ions  took 
place  before  kinetic  data  were  collected. 
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Results 

Effect  of  Precursor  Neutrals:  CyHy"^  ions  produced  from 
different  precursors  exhibited  behavior  Indicative  of  both 
reactive  and  unreactive  populations  toward  the  precursors. 
Figures  7.2a.  7.2b.  and  7.2c  show  normalized  ion  intensities 
of  CrjHrf      and  product  ions  vs  time  for  reactions  with 
precursor  neutrals;  toluene,  norbornadiene  and 
cycloheptatrlene  respectively.   Following  an  exponential 
decay  indicative  of  first  order  kinetics,  a  substantial 
fraction  of  unreactive  ions  remains  at  long  reaction  times. 
Table  7.1  shows  the  charge  transfer  agents  used  to  produce 
CyH^   ions  from  different  precursors,  the  initial  internal 
energy  imparted  to  the  ions,  and  the  percent  of  unreactive 
ions  remaining  at  long  reaction  times. 

The  major  reaction  of  C^Hy'*'  produced  by  charge  transfer 
from  toluene  with  the  toluene  neutral  gives  CgHg"*"  (reaction 
(7.2));  when  norbornadiene  is  used  as  a  precursor  and 
neutral  reactant,  reaction  (7.3)  is  observed. 


CrjHrj  +       CnHo 

Crjnrj  +       CnHo 


>  CgHg     +   CgHg 

>  CgHg     +   CgHg 


(7.2) 
(7.3) 


p:-i-. 


CiyH.^   from  norbornadiene  also  reacted  to  give  small  amounts 
of   higher  molecular  weight  products.   On  the  other  hand, 
CrjU'j'*'    produced  from   cycloheptatrlene  reacted  with  its 
precursor   neutral  to  give  significant  amounts  of  higher 
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TABLE  7.1 


Percentages  of  unreactlve  CrjHrj'^    found  from  various  precursors  by 
different  ionization  methods  monitored  by  observing  reaction  with 
the  precursor  neutral.^  .   , 


Precursor 
Neutral 


Ionizing  charge 
transfer  agents 


Maximum  Internal 
Energy  (eV)^ 


Percent 
Unreactlve 


Toluene 


Xe  + 


Kr^ 


Ar 


1  .4 


3.3 


5.0 


45-50 


25-30 


25-30 


Norbornadiene 


^4^2^ 


Xe^ 


Kr^ 


0.9 


2.5 


4.4 


70-75 


30-35 


25-30 


Cycloheptatrlene 


C^Hg"" 


Xe^ 


Kr^ 


1.1 


2.8 


4.6 


70-75 


35-40 


35-40 


For  those  ions  produced  by  C^Hg*  charge  transfer,  reactions  with 
both^precursor  neutral  and  diacetylene  were  monitored.  ^Based  on 
CrjHrj      appearance  energies  given  Traeger  and  McLoughlin,  (1978). 
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molecular  weight   product  Ions  along  with  smaller  amounts  of 
CgHg"*"  and  CgHg"*"   product  Ions: 


CyHy  +  C^Hg  >  CgHg           +      CgHg 

C7H/  +  C7H8 >  CgHg*      +      C5H6 

CrjErj  +  C^Hg  >  CjQHg          +      C4He 

C7H/  +  C7Hg >  CiiHg-'  +  CgHg 

C^H^"^  +  CyHg >  CjgHii*  +  C2H4 

CrjErj*  +  CyHg >  C  j  3  H  J  j  "^       +       CH4 

C^Hy"^  +  C7H8 >  Cj4Hii"*"  +  2H2 

CyHy  +  Cr^Hg >  C24Hjg     +   Ho 


(7.4) 
(7.5) 
(7.6) 
(7.7) 
(7.8) 
(7.9) 
(7.10) 
(7.11) 


Deuteratlon  Studies:   In  order  to  further  Investigate 
the  mechanisms  of  some  of  the  reactions  shown  above, 
deuterated  toluene  was  used  both  to  produce  C^D^"^  for 
reactions   with  neutral  precursors  and  as  a  neutral  reactant 
for  reactions   of  C^H^*  produced  from  different  precursors. 
Reactions  of  Cr^Drj'*'   with  toluene  and  norbornadiene  were 
respectively : 


CrjDrj*    +  C^Hg  (toluene) >  CgH^Dg"^  +  CgHDg     (7.12) 

CyD^""  +  C^Hg  (NBD) >  CgD^Hg*  +  CrHg  (7.13) 


When  CrfDrj      reacted  with  cycloheptatr lene ,  the  major  product 
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observed  was  unreactlve  Ci^Hy"^  due  to  reactions  (7.14)  and 
(7.15).  I 


CyD,""  +  C^Hg >  C^Hg''  +  C7D7  (7.14) 

C^Hg"^  +  C^Hq >  CyH^"^  (unreactlve)  +  C^Hg      (7.15) 

Other  reactions  observed  were: 

I 

C7D7*  +  C^Hg  >  CgDgHg^  +  C5H5D  (7.16) 

C7D/  +  C^Hg  >  CjoDgHg*  +  C4H5D  (7.17) 

C^D^^  +  C^Hg  >  CjiDgHg^  +  C3H5D  (7.18) 

C7D7*  +  C^Hg >  C12D7H4-'  +  C2H4  (7.19) 

C.yH.j,"'"  Ions  produced  from  toluene  reacted  with   deuterated 
toluene  as  shown  by  equation  (7.20): 

C7H/  +  C^Dg  >  CgD7H2^  +  CgHgD  (7.20) 


Very  low  intensities  of  this  product  ion  were  also  observed 
in  the  mass  spectra  of  CyHy"*"  ions  produced  from 
cycloheptatriene  and  norbornadlene .   The  rate  coefficients 
for  the  disappearance  of  CjUrj* /CrfDrj*    by   reactions  with 
precursor  neutrals  were  calculated  as  described  earlier  and 
the  values  are  shown  in  Table  7.2. 
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TABLE  7,2 


Rate  coefficients  for  different  C^H^^^/C^Hq  reaction  systems 


jPrecursor   j        k 

/lO-lO  cm3s-l                       1 

1      \  Ion  1 

1       \     1  Toluene- 

|Neutral\    j 

1 \__| 

dg      1      NBD       1      CHT      1 

I  Toluene-dgj  (3.5  +  i 

.0)3   1  (0.8  +  0.3)  1  (1.7  +  0.5)  1 

1    NBD       1(3.2+1.0)     1(5.4+1.7)1               | 

1    CHT       1(4.6+1.5)     1                1(8.2+2.6)1 

^Previously  reported  as  4.0  x  lO'^^  cm^/s  for  the  reaction 
of  benzyl  ions  with  toluene  (Dunbar,  1975). 
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CyHy   Reactions  with  Acetylene:  CiyHy*  Ions  produced  by 
Xe    dissociative  charge  transfer  of  norbornadlene  and 
cycloheptatrlene  reacted  slowly  with  acetylene  to  produce 
very   small  amounts  of  CgHg"*^. 


C7H7   +  CgHg 


>   CgHg 


(7.21) 


*"■'.■ 


The  following  reaction  was  also  observed  with  the  ions 
produced  from  cycloheptatrlene: 


CyHn    +   CoH 


2"2 


^  ^7^*0   **"  C ' 


(7.22) 


The  rate  coefficients  for  the  disappearance  of  C^H^*  ions 
produced  from  norbornadlene  and  cycloheptatrlene  were  found 
to  be  (1.0  +  1.0)  *  10~11  and  (3.6  +  1.3)  x  10"^^  cm^/s 
respectively.   However,  when  toluene  was  used  as  a  precursor 
for  C^Hy"^  production  by  Kr^  dissociative  charge  transfer,  no 
reaction  with   acetylene  was  observed. 


CyHy-t-  Reactions  with  Diacetylene:  C7H7*  ions  from 
different   precursors  were  reacted  with  both  diacetylene  and 
deuterated   diacetylene  in  order  to  identify  the  reaction 
products  from  those  with  the  precursor  neutral  (for  instance 
the  reactions  with  both  diacetylene  and  cycloheptatrlene 
give  CjjHg  ).   The  following  reactions  of  CyH^^"^  with 
diacetylene  were  observed. 
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CirHiT  +       ^4"0 


'7"7 


'4"2 


CitHit       +    C^Do 

CgHy  +        C4H2 


>       CgH^  +       C2H2 


C^H^"^    +    C^Do >    CqHcDo"^    +    CoH 


9"5"2 


■2"2 


^7^7       ^    ^4^2 ^    ^11^9 


>       CjgHg 


CgHgDg^  +  C4D2  ->  Ci3H5D4^ 


(7.23a) 
(7.23b) 
(7.24a) 
(7.24b) 
(7.25a) 
(7.25b) 


In  the  toluene  case,  the  following  reaction  was  also 
observed: 


CgHg     +   C4H2 
CgHg    +   C4D2 


■>  CjgHg    +  H2 

>  C22H7D2   +  H2 


(7.26a) 
(7.26b) 


Ion  Intensity  vs  time  curves  of  C^H-y^  (produced  by  Kr"^ 
charge  transfer  from  toluene)  and  product  ions  for  the 
reaction  with  diacetylene  are  shown  in  Figure  7.3.   The  rate 
coefficients  for  the  disappearance  of  C^H^"^  by  the  sum  of 
reactions  (7.23)  and   (7.24)  for  different  precursors  are 
given  in  Table  7.3. 

When  CrjHrj  ions  were  produced  in  a  100ms  reaction  time 
by  Xe  dissociative  charge  transfer  from  norbornadiene  and 
cycloheptatriene  in  the  presence  of  C4H2.  production  of  - 
unreactive  C^Hy"^  was  observed  by  the  following   reactions. 


Xe*  +  C4H2 


>  C4H2   +  Xe 


(7.27) 


o 

u 


o 
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C7H7    +   C4H- 


Figure  7.3.   Reactions  of  C^H^"'  Ions  Formed  from  Toluene  by 
Kr   Chemical  Ionization  with  Dlacetylene. 
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TABLE  7.3 
Reaction  rate  coefficients  for  C^Hy"^  +  C4H2  reaction. 


Precursor 


k/10-^0  cm^s"^ 


Toluene 


(  1  .50  +  0.50  ) 


Norbornadlene 


(  0.25  +  0.19  ) 


Cycloheptatrlene 


(  0.75  +  0.25  ) 
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C^H,   +  C^Hg >  C,7H7"^(unreactlve)  +  C4H2  +  H   (7.28) 


'4"2 


Thus  In  order  to  obtain  accurate  information  on 
reaction  kinetics  of  the  C^h^''  +  C4H2  reaction,  the  C^K^^ 
ion  was  continuously  ejected  during  C^H^^  ion  formation  by 
dissociative  charge  transfer  using  Xe"^ . 

Also  in  this  study,  the  relative  abundances  of  the 
ions  produced  from  toluene  as  a  function  of  electron  Impact 
energy  were  determined  for  the  purpose  of  comparison  with 
results  obtained  for  norbornadlene  and  cycloheptatrlene  (see 
Chapter  6).   Relative  abundances  of  C5H5*  and  CgHg"^  as  a 
function  of  energy  for  toluene  are  shown  in  Figure  7.4. 
Results  of  the  earlier  study  given  in  Chapter  6  for 
norbornadlene  and  cycloheptatrlene  are  also  included  in  the 
figure  for  comparison. 


Discussion 


Since  numerous  earlier  investigations  ^Dunbar.  1975; 
Jackson  et  al . .  1977;  Traeger  and  McLoughlln.  1978; 
McLafferty  and  Bockhoff,  1979;  Sharma  and  Kebarle,  1981) 
have  demonstrated  that  C^H^""  ions  have  at  least  two  stable 
structures,  one  reactive  and  the  other  unreactlve,  the 
results  of  this  study  are  interpreted  in  terms  of  the 
possibility  of  other  stable  structure ( s ) . 
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The  nature  of  the  reaction  products  observed  for  the 
reactions  of  C^Hy"^  Ions  with  their  precursor  neutrals, 
toluene,  norbornadlene  and  cycloheptatrlene .  were 
significantly  different,  as  shown  in  Figures  7.1a,  7.1b  and 
7.1c  respectively.   In  order  to  determine  whether  these 
differences  were  due  to  the  neutral  or  ion  structures 
deuterated  toluene  was  used  to  monitor  the  reaction 
mechanisms  with  different  neutral  precursors.   As  seen  in 
equations  (  7  .  12  )  -  (  7  .  19  )  ,  C^D^"^  ions  reacted  with  the 
precursor  neutrals,  producing  deuterated  analogues  of  the 
same  product  Ions   observed  for  each  non-deuterated 
Ion/precursor  reaction  system.   This  result  showed  that  the 
differences  in  reaction  products  are  the  result  of  the 
different  structures  of  neutral  precursors,  and  not 
necessarily  of  different  ion  structures.   Then  deuterated 
toluene  was  used  as  a  reactant  neutral  for  C^H^*  ions 
produced  from  different  precursors.   As  seen  from  equation 
(7.20),  CyHy   ions  produced  from  all  three  precursors 
reacted  with  C^Dq  to  produce  CgDyHg"^,  although  very  low  ion 
intensities  were  produced  with  cycloheptatrlene  and 
norbornadlene  precursorsors .   The  reaction  mechanism  which 
gives  this  product  ion  is  also  operative  in  reaction  (7.12) 
and  is  believed  (Shen  and  Dunbar,  1974;  Dunbar,  1975; 
Jackson  et  al.,  1977;  McCreary  and  Preiser,  1978)  to  be 
indicative  of  the  benzyl  structure  for  C^H^"^  ions  (see 
Figure  7.5).   This  mechanism  has  also  been  studied  earlier 
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Igure  7.5.  Proposed  (Shen  and  Dunbar.  1974;  Dunbar  1975. 
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by  similar  Isotoplc  labeling  experiments  (Jackson  et  al . , 

1977;  Ausloos  et  al . .  1980).   The  only  Important  difference 

in  the  behavior  of  reactive  C^H^"*"  Ions  produced  from 

I 
different  precursors  was  their  reaction  rates  with  precursor 

neutrals.   C^h^'^/C^D^'^  Ions  produced  from  toluene/deuterated 

toluene  reacted  at  similar  rates  with  toluene,  norbornadlene 

and  cycloheptatrlene  whereas  C^H^"^  Ions  produced  either  from 

norbornadlene  and  cycloheptatrlene  reacted  at  much  slower 

rates  with  toluene  than  with  their  own  precursor  neutrals 

(see  Table  7.2). 

Reactions  of  Ci^Hy"^  and  product  ions  with  dlacetylene 

are  represented  by  equations  ( 7 . 23 ) - ( 7 . 25 )  Independent  of 

the  precursor  used.   An  interesting  mechanism  is  observed 

for  Cg  and  C^  addition  reactions  shown  by  equations  (7.23b) 

and  (7.26b)  which  is  indicative  of  an  ion/molecule  complex 

from  which  CgHg  and  Hg.  respectively,  rather  than  CgDg  and 

Dg  are  the  apparent  neutral  products.   Since  complete 

scrambling  of  product  ions  is  not  observed,  preferential 

retention  of  CgDg  and  Dg  is  probably  due  to  a  specific 

structure  of  the  Intermediate  complex  which  does  not 

rearrange  fast  enough  to  give  complete  scrambling  before 

falling  apart  to  give  the  products.   Reaction  rates  with 

dlacetylene   followed  a  different  trend  than  those  with 

precursor  neutrals  for   Cp^H^^  ions  from  different 

precursors.   In  this  case,  CyHi^"^  ions  from  toluene  reacted 
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faster  than  those  from  norbornadlene  and  cycloheptatrlene 
(see  Table  7.3). 

Different  behavior  was  observed  when  C2H2  was  used  as  a 
reactant  neutral.  CjHrj'^    from  norbornadlene  and 
cycloheptatrlene  reacted  slowly  with  acetylene  as  shown  by 
equations  (7.21)  and  (7.22).   On  the  other  hand,  no  reaction 
products  were  observed  with  C^Hy"^  from  toluene. 

Table  7.1  shows  the  percent  of  unreactlve  C-^Hy"*"  as  a 
function  of  the  Internal  energy  of  the  Ion  formed  from 
different   precursors.   At  the  lowest  Internal  energies 
above  the  CyH^"*^  formation  threshold,  C7H7*  from  toluene  had 
a  much  smaller   percentage  of  unreactlve  Isomer  than  C^Hy"^ 
from  norbornadlene  and  cycloheptatrlene.   Formation  of 
unreactlve  CrjErj'*'    from  all  three   precursors  decreased  with 
increasing  internal  energy  and  stayed  constant  at  high 
enough  energy.   This  behavior  is  different  from   that 
observed  previously  for  toluene  and  norbornadlene  (Dunbar, 
1975;  Ausloos,  1982)  in  which  the  unreactlve  percent  first 
decreased  and  then  increased   with  Internal  energy.   The 
previous  results  were  interpreted  In  terms  of  the  interplay 
between  the  two  equilibria  involving  different  C^Hg"^  and 
C^Hy"^  structures  (Sagl  et  al  .  ,  1974;  McLafferty  and 
Bockhoff.  1979;  Andrews  and  Keelan,  1981;  Ausloos,  1982). 
The  different  behavior  observed  in  an  earlier  ICR  study 
(Jackson  et  al . ,  1977),  which  is  similar  to  that  observed  in 
this  work,  was  suggested  (McLafferty  and  Bockhoff,  1974; 
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Ausloos,  1980)  to  be  the  result  of  the  displacement  of  the 

equilibrium  between  two  C^H^"^  structures  by  the  removal  of 

the  reactive  CyHy"^  Ions  through  chemical  reaction. 

In  Chapter  6,  the  relative  abundances  of  ions 

produced  from  norbornadiene  and  cycloheptatriene  as  a 

function  of  ionizing  electron  energy  were  reported. 

Similarities  of  C^H-y"^.  CgHg"^  and  CgHg"^  appareance  curves 

were  interpreted  to  imply  that  C-^H^"^  formed  from  these  two 

compounds  behaves  similarly,  i.e..  follows  the  same 

fragmentation  pathways  in  the  same  energy  range. 

Differences  in  the  appareance  curve  of  CgHg*  for     '  '■ 

norbornadiene  and  cycloheptatriene  were  Interpreted  as 

suggesting  that  CyHg*  formed  from  these  compounds  by 

electron  impact  may  have  different  structures.   Extensive 

fragmentation  of  C^Hg"^  to  give  C^Hg"^  was  also  reported 

earlier  for  norbornadiene  (Ausloos,  1982).   In  this  study, 

the  relative  abundances  of  the  ions  produced  from  toluene  as 

a  function  of  energy  were  determined  for  the  purpose  of 

comparison  with  earlier  results  obtained  for  norbornadiene 

and  cycloheptatriene.   In  contrast  to  the  observed 

similarity  mentioned  above  between  norbornadiene  and 

cycloheptatriene  fragment  ion  appareance  curves,  those  of 

toluene  were  significantly  different  both  in  terms  of 

I 
energetics  and  of  general  shape.   For  example,  the  onsets  of 

formation  for  both  CgHg-"  and  CgHg-"  were  different  than  those 

values  observed  for  the  other  two  precursors.   Furthermore, 
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the  relative  Intensities  for  CgHg"^  and  CgHg"^  Increased  much 
more  slowly  with  energy  relative  to  those  from  norbornadlene 
and  cycloheptatrlene.   The   maximum  relative  Intensities  for 
CgHg"^  and  03113'^  Ions  were  respectively  0.017  and  0.005  for 
toluene  and  0.10  and  0.035  for  norbornadlene  and 
cycloheptatrlene  (see  Figure  7.4).   This  different  behavior 
In  the  fragmentation  of  C^Hy"*^  produced  from  toluene  suggests 
the  Involvement  of  a  different  structure  (benzyl"^)  for  C^H^"^ 
In  the  case  of  toluene. 

A  possible  potential  energy  surface  for  CyH^^/CYHg"^ 
system  Is  shown  In  Figure  7.6.   Based  on  the  differences  In 
behavior  observed  In  this  study  for  CjH,^'*^    Ions  produced  from 
three  precursors,  the  possibility  of  stable  8tructure(s) 
other  than  benzyl  and  tropyllum  must  be  considered.   The 
unreactlve  structure  which  was  observed  In  all  three  cases 
has  been  shown  to  be  cyclic  tropyllum  by  a  number  of 
previous  experimental  (see  for  example,  Jackson  et  al . , 
1977)   and  theoretical  (Abbaund  and  Hehre,  1976;  Cone  et 
al.,  1977)  studies.   Similarly,  the  reactive  Isomer  produced 
from  toluene  was  Identified  as  the  benzyl  structure  numerous 
times  previously  (see  for  example,  Ausloos  et  al . ,  1980; 
Bombach  et  al.,  1983).   However,  the  results  of  this  study 
have  shown  that  it  behaves  somewhat  differently  from  the 
reactive  C^Hy"^  isomers  produced  from  norbornadlene  and 
cycloheptatrlene.   Although  no  definitive  assignment  for 
reactive  CjHrj'*^    ion  structure(s)  produced  from  norbornadlene 
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and  cycloheptatriene  can  be  made  at  this  point,  the  most 
probable  structures  will  be  discussed  considering  the 
energetics  of  CyHr^"^  Ions. 


Cycloheptatriene.   AMI  calculations  for  the  C^Hg'^/CyH-y'*' 
surface^^  suggest  the  structure  produced  by  hydrogen  loss 
from  the  blcyclo  [3  .  2  .  0]heptadlene  radical  cation  (BCH"*^),  is 
a  possible  stable  C^Hy"^  isomer,  accessible  at  about  2.8  eV 
above  the  appearance  potential  of  tropylium  from 
cycloheptatriene  (see  Figure  7.6).   This  structure  is  in  a 
sufficiently  deep  potential  well  that  it  might  reasonably  be 
expected  to  live  for  an  extended  time.   It  was  also 
suggested    that  this  structure  is  most  probably  formed  from 
cycloheptatriene  by  isomerization  on  the  radical  cation 
surface  followed  by  hydrogen  loss  from  (BCH*). 
Isomerization  of  (BCH-H)"^,  ( AHf  =  274  .  4kcal/mol  )  ^^  to 
tropylium  is  considered  unlikely  due  to  the  high  energy 
barrier  between  two  structures  as  seen  in  Figure  7.6. 
Although  the  barrier  to  form  the  toluene  Ion  is  lower  than 
that  for  (BCH""),  it  was  suggested^^  that  the  simplicity  of 
the  rearrangement  to  the  latter,  compared  to  the  convoluted 
surface  leading  to  benzyl  or  toluene  radical  ion  may  make 
formation  of  BCH"*^ ,  a  faster  rearrangement  at  high  enough 
energy.   Consistent  with  these  calculations,  our 
experimental  results  show  that  the  percent  of  unreactive 
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C^Hy   from  cycloheptatrlene  decreases  very  rapidly  with 
Increasing  Internal  energy.   This  can  be  Interpreted  as  a 
result  of  increased  Isomerlzation  of  the  cycloheptatrlene 
cation  to  (BCH"^).   At  high  energy,  when  the  equilibrium 
between  these  two  parent  ions  Is  established,  the  percent  of 
unreactive  Isomer  stays  constant  due  to  the  absence  of 
isomerlzation  between  two  C^H^^  structures  with  an  energy 
barrier  of  14.1  kcal/mol^^.   Other  experimental  evidence  for 
the  involvement  of  different  structures  for  C^Hg"^  ions 
formed  from  toluene  and  cycloheptatrlene  comes  from  the 
strikingly  different  photodissoclation  curves  of  these  ions 
(Dunbar  and  Fu,  1973). 

Norbornadlene.   In  spite  of  the  general  similarity  of 
the  behavior  of  CjHj'*'    ions  from  cycloheptatrlene  and 
norbornadlene,  differences  in  reactivity  (see  Tables  7.2  and 
7.3)  suggest  the  formation  of  a  C^H^"^  structure  different 
from  that  discussed  above  for  cycloheptatrlene.   A  simple 
bond  breakage  in  norbornadlene  can  directly  form  an 
Ipso-protonated  benzyl  radical.  C6H5(H'')CH2  ion.  in  which  a 
1.2-hydrogen  shift  leads  to  the  Intermediate.  5-methylene- 
1  ,3-cyclohexadiene  ion.  (MCH"^)  with  a  theoretical   AHj  =  219 
kcal/mol  (Dewar  and  Landman,  1977)  (see  Figure  7.6).   This 
ion  was  proposed  to  account  for  the  ring  expansion  of  the 
.,  toluene  ion  in  the  toluene"^ >  cycloheptatrlene* 
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isomerization  process  (Dewar  and  Landman,  1977).   The 
formation  of  (MCH^)  from  other  precursor  cations  at  least 
with  a  transitory  existence  has  been  indicated  by  previous 
experimental  evidence  (Leusen  et  al . ,  1973).   After  the 
availability  (Gajewskl  and  Gortva,  1982)  of  5-methylene-l , 3- 
cyclohexadlene,  experimental  determination  of  its  ionization 
potential  indicated  a  heat  of  formation  value  18  kcal/mol 
higher  than  that  of  the  toluene  ion  (Bartmess,  1982). 

Once  this  intermediate  ion  is  formed,  it  can  rearrange 
to  another  structure  or  decompose  directly  to  the  C^H™"^  ion. 
In  fact,  it  was  suggested  that  direct  decomposition  of  this 
Ion  is  in  competition  with  rearrangement  to  other 
structures,  cycloheptatr lene  cation  being  the  most  probable 
rearrangement  product  (Bursey  et  al . ,  1973).   The 
competition  between  these  two  processes  is  believed  to 
depend  on  the  internal  energy  with  rearrangement  being 
favored  at  lower  energies  (Dunbar  and  Klein,  1977).   Thus  at 
low  energies,  fragmentation  to  C^H,^^  after  rearrangement  to 
cycloheptatrlene  but  not  to  toluene  ion  is  expected  due  to 
the  lower  activation  barrier  for  fragmentation  of 
cycloheptatrlene  ion  relative  to  toluene  parent  ion  (Dunbar, 
1975).   This  is  consistent  with  our  results  which  indicate  a 
high  percent  of  unreactive  C^H^*    ( tropyl  ium"^ )  formed  from 
norbornadlene  at  low  energy  (see  Table  7.1).  Furthermore,  as 
a  result  of  the  photodlssoclatlon  study  of  Ci^Hg*  structures 
(Dunbar  and  Fu,  1973),  possibility  for  conversion  of  about 
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half  the  norbornadlene  cations  to  the  cycloheptatrlene 
structure  was  not  inconsistent  with  the  data.   On  the  other 
hand,  It  has  been  reported  that  there  is  no  rearrangement  of 
this  ion  to  toluene  ion  on  a  time  scale  of  10"^  s  (Bursey  et 
al. ,  1971) . 

Photodissociatlon  spectra  of  n-butylbenzene  and  2- 
phenylethanol  ions  which  are  believed  to  form  (MCH"^)  were 
studied  in  order  to  identify  any  possible  rearrangement 
products  (Dunbar  and  Klein,  1977).   The  results  were 
consistent  with  the  retention  of  the  (MCH*)  structure,  with 
less  than  20  percent  rearrangement  to  cycloheptatrlene  ion. 
Although  there  was  not  any  evidence  for  ready 
interconverslon  among  C-^Hg*  structures  for  nondecomposlng 
methylenecyclohexadlene  ions.  It  was  suggested  the  that  more 
energetic,  decomposing  ions  may  interconvert  to  other 
structures.   In  this  study,  the  decrease  in  the  percent  of 
tropyllum  structure  with  the  increase  in  energy  is 
interpreted  as  a  result  of  decomposition  of  (MCH"^)  becoming 
competitive  with  rearrangement  to  the  cycloheptatrlene  ion. 
Thus,  the  CyHi^"^  structure  formed  by  hydrogen  loss  from 
(MCH  )  can  be  suggested  as  the  reactive  Isomer  observed  for 
norbornadlene.   The  overall  results  of  this  study  in  terms 
of  the  proposed  structures  for  CyHi^"^  ions  formed  from 
different  precursors  are  given  in  Table  7.4. 


161 


TABLE  7.4 


Proposed  CyH,^*  structures  from  different  precursors 


1           1  Unreactive 
IPrecursorj  structure 

1  Reactivel    AH^     J 
1  structure} (kcal/mol) | 

JToluene   j  Tropylium"^ 

1  Benzyl'*'  |    211     | 

1   NBD     1  Tropyllum"*^ 

1  (BCH-H)"^^!    271      1 

1   CHT     1  Tropylium"' 

|(MCH-H)''^|   ?(236)   1 

^Structure  produced  by  hydrogen  loss  from 
bicyclo[3.2.0]heptadlene. 

Structure  produced  by  hydrogen  loss  from  5-methylene-l , 3- 
cyclohexadlene  Ion. 


^. 
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An  ICR  study  (Bartmess.  1982)  Involving  the  formation 
of  (MCH)"^  from  Its  precursor  neutral  at  30eV  electron  energy 
indicated  Ion/molecule  chemistry  different  from  that  of 
toluene.   C^H^"^  produced  from  (MCH"^)  was  reported  to  have  a 
different  rate  coefficient  and  reaction  efficiency  for  the 
formation  of  CgHg"".   It  Is  Interesting  to  note  that  the  rate 
coefficient  reported  (0.6  x  10~^°  cm^/s)  (Bartmess,  1982) 
for  this  reaction  agrees  very  closely  with  the  value,  (0.8  + 
0.3)  X  10  -^^  cm^/s  found  In  this  work  for  the  reaction  of 
C^H^   formed  from  norbornadlene  with  toluene  (see  Table 
7.2).   Although  a  definitive  structure  can  not  be  assigned 
to  CyHy   formed  from  (MCH'*')  yet,  one  possible  structure  is 
the  norcaradlenyl  Ion,  which  was  calculated  to  be  stable 
with  AHf  =  235.9  kcal/mol  (Cone  et  al  .  ,  1977)  (see  Figure 
7.6).  I 

Supporting  experimental  evidence  for  the  existence  of 
stable  CrjHrj'*'    structures  other  than  benzyl  and  tropyllum 
comes  from  a  UV  absorption  study  of  C-yH^"^  Ions  In  solid 
argon  (Andrews  and  Keelan,  1981).   In  relation  to  rate 
coefficient  data  for  reactions  with  precursor  neutrals, 
acetylene,  and  diacetylene  the  same  trend  was  observed  for 
norbornadlene  and  cycloheptatrlene .   In  all  three  cases,  the 
reactive  structure  formed  from  norbornadlene  reacted  more 
slowly  than  that  formed  from  cycloheptatrlene.   This  Is 
consistent  with  the  general  observation  that  less  reactive 
structures  of  small  hydrocarbon  ions  have  lower  heats  of 
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formation  since  the  C.^H^'^  structure  formed  from  (MCH*)  has  a 
lower  heat  of  formation  than  that  formed  from  (BCH"^). 
Finally,  In  relation  to  the  Ionic  soot  formation 
mechanism,  the  existence  of  a  different  Cf^H^"^  Ion  less 
stable  than  benzyl  In  flames  would  account  for  the 
disagreement  between  experiment  and  theory  In  relation  to 
CfjErj'*'    concentration  in  flames  (Olson  and  Calcote,  1981a). 


CHAPTER  8 
CONCLUSIONS  AND  RECOMMENDATIONS 


The  experimental  work  reported  here  contributes 
significantly  to  an  understanding  of  the  reactivity, 
energetics,  and  structures  of  small  hydrocarbon  ions. 
However,  it  does  not  prove  or  disprove  an  ionic  mechanism  of 
soot  formation.   It  does  lead  to  several  clarifications  of 
the  ionic  mechanism,  which  will  be  discussed  below. 

The  extensive  ion/molecule  condensation  reactions  at 
rates  approaching  the  Langevln  limit  observed  when  CgHg"^, 
CsHs"^  and  C^H^"^  react  with  C^Hg  suggest  that  this  aspect  of 
the  proposed  (Calcote,  1981;  Olson  and  Calcote,  1981a)  ionic 
path  to  soot  formation  is  quite  credible.   The  general 
reaction 


C^H/    +  C4H2  — >  Cx^2«y^  ^  ^2"^ 


(8.1) 


appeares  to  occur  in  a  facile  manner  and  could  be  quite 
Important  in  soot  nucleation.   Some  of  the  product  ions 
formed  (e.g.  CgHg*  and  CyHg"")  have  been  seen  to  be  abundant 
In  both  nonsootlng  and  sooting  flames  (Olson  and  Calcote, 
1981b).   Furthermore,  product  ions  with  Increased  C:H  ratios 
having  as  many  as  13  carbon  atoms  have  been  observed  in  the 
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reaction  of  C^Hy"^  with  dlacetylene.   On  the  other  hand,  the 
observation  of  C3H3*  Isomerlzatlon  and  not  condensation  with 
acetylene  as  well  as  low  blmolecular  reactivity  of  acetylene 
with  C5H5   and  C^Hy"^  suggests  that  the  proposed  sequential 
acetylene  addition  reactions  to  C^h^^    In  the  Ionic  soot 
formation  mechanism  be  reconsidered.   Other  channels  such  as 
direct  reaction  of  neutral  aromatlcs  with  CgHg"^  (Mlchaud  et 
al.,  1981;  Baykut  et  al . ,  1986)  may  be  as  Important  In  the 
formation  of  small  polycycllc  Ions. 

These  experiments  were  carried  out  at  significantly 
lower  temperatures  (363K)  and  pressures  (p  <  5  X  10"^  torr) 
than  those  found  In  most  combustion  environments.   It  Is 
thus  quite  possible  that  third-body  collisions  In 
atmospheric  pressure  flames  can  stabilize  a  fraction  of  the 
reaction  complexes  leading  to  significant  termolecular 
reactivity  of  the  Ions  studied  with  acetylene.   In  fact, 
such  colllslonal  stabilization  of  the  association  complexes 
for  the  reactions  of  C3H3^  and  C4H4''  with  CgHg  has  been 
shown  to  occur  In  higher  pressure  SIFT  studies  (Smith  and 
Adams.  1987;  Knight  et  al.,  1987). 

Continuation  of  this  work  In  several  directions  could 
have  a  major  Impact  on  assessing  the  Importance  of  an  Ionic 
mechanism  of  soot  formation.   First,  experiments  should  be 
carried  out  at  higher  pressures,  using  a  high  pressure 
source  for  the  FTICR  mass  spectrometer,  to  assess  the 
Importance  of  third-body  stabilization  of  Ion/molecule 
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collision  complexes.   Actual  sampling  of  ions  from  flames 
into  the  FTICR  mass  spectrometer,  In  order  to  probe  their 
structures  and  reactivity,  would  extend  the  isolated  studies 
done  here  to  an  actual  combustion  environment. 

A  crucial  argument  in  favor  of  an  ion/molecule  soot 
nucleation  route  is  the  ease  with  which  cyclic  ions,  which 
could  lead  to  polycyclic  aromatic  hydrocarbons,  are  formed. 
Yet  the  work  reported  here  did  not  identify  definitively 
when  the  first  cyclic  ions  were  formed.   Predicted  spectra 
from  theoretical  calculations  on  C^Hg^,  C5H3'*^  and  C^Hj^"^  ions 
should  be  combined  with  laser  spectroscopic  studies  on  these 
ions  to  determine  which,  if  any,  have  cyclic  structures. 
Also,  theoretical  studies  which  probe  the  geometries  of  the 
ion/molecule  collision  complexes,  such  as  those  reported  in 
Chapter  5  for  CgHg"^  +  acetylene,  should  be  extended  to  other 
important  reactions,  e.g.  with  diacetylene,  and  the  nature 
of  the  proposed  chemi-ionization  reaction  which  forms  the 
C3H3'*"  ion  (see  Chapter  1)  should  be  investigated. 


APPENDIX  I 

PROGRAM  TO  CALCULATE  ABSOLUTE  RATE  CONSTANTS  AND  THEIR  95% 
CONFIDENCE  LIMITS  FROM  RAW  OR  NORMALIZED  INTEGRATED  PEAK 
AREAS  OF  THE  REACTING  ION  AS  A  FUNCTION  OF  TIME  IN  FOURIER 
TRANSFORM  ION  CYCLOTRON  RESONANCE  MASS  SPECTROSCOPY. 

WRITTEN  JULY  AND  AUGUST  1986  BY  FEZA  OZTURK  AND  BRYAN  HEARN 
Adapted  for  use  with  the  IBM  pc  on  9-17-87  by  CLIFF  WATSON. 

UNDER  THE  DIRECTION  OF  DR.  JOHN  EYLER 
DEPARTMENT  OF  CHEMISTRY 
UNIVERSITY  OF  FLORIDA 
GAINESVILLE,  FLORIDA   32611 


TYPICAL  ORDER  OF  EXECUTION: 

A.  Calculation  of  the  system  factor  using  the  literature 
value  of  a  well-known  rate  constant  for  a  reference  reaction 
(this  Involves  calculating  the  observed  rate  constant  for 
the  reference  reaction,  #3,  the  actual  pressure,  #1,  then 
the  system  factor,  #4,  which  uses  the  literature  value  for 
the  absolute  rate  constant  along  with  your  experimental 
observed  rate  constant  and  actual  pressure  Just  calculated 
using  #3  and  *1).   The  system  factor  will  be  used  later  in 
the  absolute  rate  constant  calculation,  so  it  should  be 
noted.   If  this  calculation  is  not  desired,  a  value  of  1  may 
be  input  for  the  system  factor. 

B.  Running  the  program  a  second  time,  calculate  the 
observed  rate  constant  for  the  reaction  of  the  ion  with  its 
precursor  neutral,  *3 . 

Step  3  should  be  followed  by  step  9  here  if  the  no 
baseline  correction  option  is  chosen  in  the  calculation  of 
k. 

C.  Calculation  of  the  observed  rate  constant  for  the 
reaction  of  the  ion  with  reactant  neutral  (gives  the  sum  and 
individual  rate  constants  with  precursor  and  reactant 
neutrals,  *  5).   This  should  be  followed  by  step  10  if  the  . 
no  baseline  correction  option  is  chosen  in  the  calculation 


of  k 


2- 


D.  Determination  of  the  actual  pressure  for  the  reactant 
neutral  by  entering  the  ionization  gauge  and  capacitance 
manometer  readings  manually  (Baratron  factor  is  determined 
from  the  slope  by  linear  least  squares,  #  1). 

E,  Calculation  of  the  absolute  rate  constant  using  the 
actual  pressure  and  system  factor  (#  6). 
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For  steps  B  and  C,  raw  or  normalized  Integrated  peak  areas 
of  the  reacting  Ion  as  a  function  of  time  should  be  entered 
manual ly . 

System  factor  Is  a  correction  factor  for  pressure 
calibration  due  to  the  fact  that  Ionization  gauge  and 
capacitance  manometer  are  located  at  different  points  of  the 
instrument . 

Options  2,  8,  11  and  12  can  be  used  to  Input  the 
appropriate  values  once  they  are  calculated.   This  allows 
one  to  calculate  absolute  rate  constants  without  going 
through  the  entire  process  by  entering  the  precalculated 
values . 

I 

Baseline  correction  is  needed  to  correct  for  the 
unreactive  population  of  the  reacting  Ion  when  both  reactive 
and  unreactive  isomers  exist  (#'s  9  and  10  are  used  when 
there  is  no  unreactive  isomer  ) . 


te: 


DIMENSION 
RDATA{50,50) .TIMES (100) .TIC (50) .TMPDAT(50) .TMP2DT(50) 

DIMENSION 
YP(IOO) .XP(IOO) .TEMP (50) .ETEMP(50) ,CMA(100) ,TXGRAF(52) . 
&TYGRAF(52) 
REAL  K,IG.LITK.K2 
INTEGER  FLAG 
CHARACTER*1  YESNO 
10    WRITE(*.*) '   • 

«  )  I  I    i 

*)•  1  PRESSURE  CALC  (MANUAL)  2  ENTER  BF  ' 
*)'  3  K  AND  95*  CONFIDENCE    4  SYSTEM  FACTOR 
*)'  5  K2  AND  95%  CONFIDENCE   6  ABSOLUTE  RATE 


WRITE(* 
WRITE(* 
WRITE(* 
WRITE(* 

CONSTANT    ' 
WRITE(* 
WRITE(* 

C0RR(K2)  • 

WRITE(* 

95*  CL 

WRITE(* 
WRITE(* 


*)'  7  STOP  8  INPUT  SYSTEM  FACTOR 

*)'  9  NO  BASELINE  CORR(FOR  K)  10  NO  BASELINE 


*) '  11  ENTER  K  AND  95*  CL 


12  ENTER  K2  AND 


*)• 
*)' 


INPUT  CHOICE 


READ(*,*)I 

GOTO (100. 200. 300. 4 00. 500. 600. 700. 800. 900. 1000,  1300. 
&1400)I 

GOTO  10 
100   WRITE(*,*)*  NUMBER  OF  POINTS  TO  ENTER  MANUALLY  FOR 
PRESSURE  ' 

READ(*.*)NPRESS 

DO  123  I=1,NPRESS 

WRITE(*. 120)1 
120   FORMAT('  INPUT  IG  READING  #  '.15) 

READ(*.*)XP(I) 
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WRITE(*, 129)1 
129   FORMATC  INPUT  BARATRON  READING  #  '.15) 
123   READ(*.*)YP(I) 

CALL  GRAF2D(XP,YP,NPRESS) 

WRITE(*.*)'  INPUT  THE  NUMBER  OF  POINTS  FOR  LLSQ  ' 

READ(*,*)NPRESS 

FLAG=3 

CALL  LLSQ(XP.YP,EYP.NPRESS.BF,EBF, FLAG. INT) 

DO  196  I=1,NPRESS 

TYGRAF(I)=YP(I) 
196   TXGRAF(I)=XP(I) 

TXGRAF(NPRESS+1)=XP(NPRESS) 

TYGRAF{NPRESS+1)=BF*XP(NPRESS) 

TXGRAF(NPRESS+2)=XP(1 ) 

TYGRAF(NPRESS+2)=BF*XP(1) 

CALL  GRAF2D(TXGRAF,TYGRAF, (NPRESS+2) ) 

WRITE(*,*)'  THE  BARATRON  FACTOR  IS' 

WRITE(*.*)BF 

WRITE(*,*)'  THE  95*  CONFIDENCE  LIMIT  FOR  BARATRON 
FACTOR  IS  • 

WRITE(*,*)EBF 
150   WRITE(*,*)'  INPUT  THE  IONIZATION  GAUGE  READING  ' 

READ(*,*)IG 

P=IG*BF  I 

WRITE(*,*)'  THE  PRESSURE  IS  ' 

WRITE(*,*)P 

WRITE(*.*)'  DEFAULT  95*  CON. LIMIT  FOR  IG  IS  0.1*GAUGE: 
CHANGE=Y  ' 

READ(*,88)YESN0 

IF(YESNO.EQ. ' Y' )THEN  ^ 

WRITE(*,*)'  INPUT  THE  *  MULTIPLIER  DESIRED  95*  CON. 
LIMIT  FOR  IG  • 

READ(*.*)EI 
EI=EI*IG 
GOTO  152 

END  IF 

EI=0.1*IG 
C  DEFAULT  95*  CONF. LIMIT  FOR  ION.  GAUGE  READING  IS  EST ' ED  AS 
0. 1*I0N. GAUGE 
152   EP=SQRT(BF*BF*EI*EI+IG*IG*EBF*EBF) 

WRITE(*.*)'  THE  95*  CONFIDENCE  LIMIT  FOR  PRESSURE  IS' 

WRITE(*,*)EP 

GOTO  10 
200   WRITE(*.*)'  ENTER  THE  BARATRON  FACTOR  ' 

READ(*,*)BF 

WRITE(*,*)'  ENTER  THE  95*  CONF.  LIMIT  OF  BF  ' 

READ(*.*)EBF 

GOTO  150 
300    WRITE(*.*)'   NUMBER  OF  POINTS  TO  ENTER  MANUALLY  =' 

READ(*.*)NTIM 

DO  366  1=1 ,NTIM 
WRITE(*, 364)1 
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364   FORMAT('  ENTER  TIME  #'.I5) 

READ(*.*)TIMES(I) 

WRITE(*.*)'  ENTER  THE  CIN  OR  CIN  NORMALIZED  TO  TOTAL 
ION  CIN  ' 

READ(*,»)TEMP(I) 
366   CONTINUE 
345   WRITE(*,*)'  DO  YOU  WISH  BASELINE  CORRECTION?  YES=1  ' 

READ(*.*)NANS 

IF(NANS.NE.l)  GOTO  10 

CALL  CNIC(NTIM,TEMP,ETEMP,XMAX.EXMAX) 

390  YLSTPT=TEMP(1 ) 
KOUNT  =1 

DO  391  1=2, (NTIM-1 ) 

IF(TEMP(I ) .GT. YLSTPT)GOTO  392 

YLSTPT=TEMP(I)   | 

391  KOUNT=KOUNT  +  1 

392  CALL  GRAF2D(TIMES, TEMP, KOUNT) 

WRITE(*,*)'  INPUT  THE  NUMBER  OF  POINTS  TO  USE  ' 
READ(*.*)NPOINT 

395  CALL  LLSQ(TIMES, TEMP. ETEMP.NPOINT.K.ERRK, FLAG. YINT) 
DO  396  I=1.NP0INT 

TYGRAF(I)=TEMP(I) 

396  TXGRAF(I)=TIMES(I) 
TXGRAF(NP0INT+1)=TIMES(NP0INT) 
TYGRAF(NP0INT+1)=K*TIMES(NP0INT)+YINT 
TXGRAF(NP0INT+2)=TIMES(1  ) 
TYGRAF(NP0INT+2)=K*TIMES(1)+YINT 
CALL  GRAF2D(TXGRAF.TYGRAF, (NPOINT+2) ) 
K  =  -K 

WRITE(*.*)'   K  AND  95%  CONFIDENCE  LIMIT  IN  K= ' 
WRITE(«.*)K.ERRK 
GOTO  10 
400   WRITE(*,*)'  HAVE  YOU  CALCULATED  THE  PRESSURE?  YES=1 ' 
READ(*,*)NANS 
IF  (NANS.NE. 1 )GOTO  10 

WRITE(*.*)'  INPUT  THE  LITERATURE  VALUE  K  ' 
READ(*.*)LITK 

WRITE(*,*)'  INPUT  THE  LITERATURE  ERROR  VALUE   ' 
READ(*.*)ELIT 

WRITE(*,*)'  INPUT  THE  TEMPERATURE  ' 
READ(*,*)T 

WRITE(*.*)'  THE  DEFAULT  ERROR  IN  TEMPIS  2  DEGREES. 
CHANGE=Y  ' 

READ(*.88)YESN0 
IF(YESNO.EQ. ' Y' )THEN 

WRITE(*,*)'  INPUT  THE  DESIRED  ERROR  IN  TEMP.   ' 
READ(*,*)ET 
GOTO  415 
^.]'  END  IF 

'"''  ET=2 

C  DEFAULT  ERROR  IN  TEMPERATURE  IS  ESTIMATED  AS  2  DEGREES 
|f^  '    415   K=l  .036E-19*T*K/P 
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EK=SQRT(ET*ET/(T*T)+EK*EK/(K*K)+EP*EP/(P*P) )*K 

SF=K/LITK 

ESF=SQRT(EK*EK/(K*K)+ELIT*ELIT/(LITK*LITK) )*SF 

WRITE(*,*)'  THE  SYSTEM  FACTOR  AND  ITS  95%  CONFID 
LIMIT  IS  • 

WRITE(*.*)SF.ESF 

GOTO  10 
500    WRITE(*,*)'   NUMBER  OF  POINTS  TO  ENTER  MANUALLY  = 
.   READ(*,*)NTIM 

DO  566  I=1,NTIM 

WRITE(*,564) I 
564   FORMAT('  ENTER  TIME  #',I5) 

READ(*.*)TIMES(I) 

WRITE(*,*)'  ENTER  THE  CIN  OR  CIN  NORMALIZED  TO  TOTAL 
CIN  ' 

READ(*.*)TEMP(I) 
566   CONTINUE 
545   WRITE(*,*)'  DO  YOU  WANT  BASELINE  CORRECTION?  YES=1  ' 

READ (*.*) NANS 

IF(NANS.NE. 1)G0T0  10 

CALL  CNIC(NTIM.TEMP.ETEMP.XMAX.EXMAX) 

590  YLSTPT=TEMP(1 ) 
KOUNT  =1 

DO  591  1=2, (NTIM-1) 
IF(TEMP(I) .GT.YLSTPT)GOTO  592 
YLSTPT=TEMP(I) 

591  KOUNT=KOUNT  +  1  i 

592  CALL  GRAF2D(TIMES, TEMP, KOUNT) 

WRITE(*,*)'  INPUT  THE  NUMBER  OF  POINTS  TO  USE  ' 
READ(*,*)NPOINT 

CALL  LLSQ(TIMES,TEMP,ETEMP,NP0INT,K2,ERRK2.FLAG.YINT) 
DO  596  I=1,NP0INT 
TYGRAF(I)=TEMP(I) 
596   TXGRAF(I)=TIMES(I ) 

TXGRAF(NP0INT+1 ) =TIMES ( NPOINT ) 
TYGRAF(NP0INT+l)=K2*TIMES(NP0INT)-t-YINT 
TXGRAF(NP0INT  +  2)=TIMES(1  ) 
TYGRAF(NP0INT+2)=K2*TIMES(1 )+YINT 
CALL  GRAF2D(TXGRAF,TYGRAF,  (NPOINT+2)  ) 
K2=-K2 

WRITE(*,*)'  THE  OBSERVED  SLOPES  (Kl  AND  K2 )  ARE  ' 
WRITE(*,*)K,K2 

WRITE(*,*)'  THE  ERRORS  IN  Kl  AND  K2  ARE  ' 
WRITE(*, *)ERRK,ERRK2 
598   K=K2-K 

EK=SQRT(ERRK*ERRK+ERRK2*ERRK2) 
WRITE(*,*)'  THE  NET  RATE  CONSTANT  IS' 
WRITE(*.*)K 

WRITE(*,*)'  THE  95*  CON.  LIMIT  IN  THE  OBSERVED  RATE 
CONSTANT  IS' 

WRITE('*,*)EK 
GOTO  10 
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600      WRITE(*.*) 'CALCULATED    THE    SYSTEM    FACTOR    AND    PRESSURE? 
Y=l     • 

READ(*.»)NANS 
IF(NANS.NE.1)G0T0    10 

WRITE(*,*)'     INPUT    THE    TEMPERATURE' 
READ(*,*)T 

TRUEK=1.036E-19*K*T/{SF*P) 
WRITE(*,*)'     THE    TRUE    K    IS     ' 
WRITE(*.*)TRUEK 

WRITE(*.*)'     THE    DEFAULT    ERROR    IN    TEMP    IS    2    DEGREES. 
CHANGE=Y    • 

READ(*,88)YESN0   | 
IF(YESN0.EQ. 'Y' )THEN 

WRITEC*,*)'     INPUT    THE    DESIRED    ERROR    IN    TEMP.      ' 
READ(*,*)ET 
GOTO    615 
ENDIF 
ET=2 

C    DEFAULT    ERROR    IN    TEMPERATURE    IS    ESTIMATED    AS    2    DEGREES 
615 

ETRUEK=SQRT(EK*EK/(K*K)+ET*ET/(T«T)+ESF«ESF/{SF*SF)+EP*EP/ 
&(P*P) )*TRUEK 
WRITE(*.*)'     THE    95%    CONFIDENCE    LIMIT    FOR    TRUE    K    IS     ' 
WRITE(*.*)ETRUEK 
GOTO    10 
800       WRITE(*,*)'     INPUT    SYSTEM    FACTOR' 
READ(*,*)SF 

WRITE(*.*)'     INPUT    THE    ERROR    IN    THE    SYSTEM    FACTOR     ' 
READ(«.*)ESF 
GOTO    10  I 

900       WRITE(*,*) '****N0    BASELINE    CORRECTION****' 

WRITE(*.*) 'DEFAULT    ERROR    MULTIPLIER    IN    CIN    IS    0.1 
CHANGE=Y' 

READ(*.88)YESN0 
IF(YESNO.EQ. 'Y' )THEN 

WRITE(*,*)'     INPUT    THE    DESIRED    MULTIPLIER    OF    ERROR 
IN    CIN     ' 

READ(*.*)ECIN 
GOTO    915 
ENDIF 
ECIN=0. 1 
915       DO    910    1=1 ,NTIM 

ETEMP(I)=ECIN*TEMP(I  ) 
ETEMP(I)=ETEMP(I)/TEMP(I) 
910       CONTINUE 

DO    920    I=1,NTIM 
920       TEMP(I )=LOG(TEMP(I  )  ) 
990      YLSTPT  =  TEMP(1  ) 
KOUNT    =1 
DO    991     I=2.NTIM 

IF(TEMP( I ) .GT. YLSTPT)GOTO    992 
YLSTPT=TEMP(I) 


991 
992 

CALC 


173 


KOUNT=KOUNT  +  1 

CALL  GRAF2D(TIMES,TEMP,K0UNT) 

WRITE(*,*)'  INPUT  THE  NUMBER  OF 


POINTS  TO  USE  IN  LLSQ 


READ(*,*)NPOINT  ' 
FLAG=2 

CALL  LLSQ (TIMES, TEMP. ETEMP.NPO I NT. K.ERRK. FLAG, YI NT) 
DO  996  I=1,NP0INT 
TYGRAF(I)=TEMP(I) 
996   TXGRAF(I)=TIMES(I) 

TXGRAF(NP0INT+1)=TIMES(NP0INT) 
TYGRAF(NP0INT+1)=K*TIMES(NP0INT)+YINT 
TXGRAF(NP0INT+2)=TIMES(1 ) 
TYGRAF(NP0INT+2)=K*TIMES(1)+YINT 
CALL  GRAF2D(TXGRAF.TYGRAF,  (NPOINT  +  2)  ) 
K  =  -K 

WRITE(*,*)'  OBSERVED  RATE  CONSTANT  AND  95%  CONFIDENCE 
LIMIT  IS  • 

WRITE(*,*)K,ERRK 

GOTO  10 
1000  WRITE(*,*) •****N0  BASELINE  CORRECTION****' 

WRITE(*,*) 'DEFAULT  ERROR  MULTIPLIER  IN  CIN  IS  0.1 
CHANGE=Y' 

READ(*,88)YESN0 
88    format (a) 

IF(YESNO.EQ. ' Y' ) THEN 

WRITE(*,*)'  INPUT  THE  DESIRED  MULTIPLIER  OF  ERROR 
IN  CIN  ' 

READ(*,*)ECIN 
GOTO  1015 

ENDIF 

ECIN=0.1 
1015  DO  1060  I=1,NTIM 

ETEMP(I)=ECIN*ETEMP(I) 

ETEMP ( I ) =ETEMP ( I ) /TEMP ( I ) 
1060  CONTINUE 

DO  1070  I=1,NTIM 
1070  TEMP(I)=LOG(TEMP(I)  ) 

1090  YLSTPT=TEMP(1) 
KOUNT  =1 

DO  1091  I=2,NTIM 

IF(TEMP( I ) .GT. YLSTPT)GOTO  1092 
YLSTPT=TEMP(I )    | 

1091  KOUNT=KOUNT  +  1 

1092  CALL  GRAF2D(TIMES. TEMP, KOUNT) 
WRITE(*,*)'  INPUT  THE  NUMBER  OF 


CALC 


POINTS  TO  USE  IN  LLSQ 

READ(*,*)NPOINT 
FLAG=2 

CALL  LLSQ (TIMES. TEMP, ETEMP, NP0INT,K2,ERRK2. FLAG. YINT) 
DO  1096  I=1,NP0INT 

TYGRAF(I)=TEMP(I ) 
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1096   TXGRAFd  )=TIMES(I  ) 

TXGRAF(NP0INT+1 ) =TIMES ( NPOINT ) 
TYGRAF(NP0INT+1)=K2*TIMES(NP0INT)+YINT 
TXGRAF(NP0INT+2)=TIMES(1) 
TYGRAF(NP0INT+2)=K2*TIMES(1 )+YINT 
CALL  GRAF2D(TXGRAF,TYGRAF,  (NPOINT  +  2)  ) 
K2=-K2 

WRITE(*,*) 'OBSERVED  RATE  CONSTANT  AND  95%  CONFIDENCE 
LIMITS  ARE  ' 

WRITE(*,*)K2.ERRK2 

GOTO  10 
1300  WRITE(*,*)'  ENTER  OBSERVED  RATE  CONSTANT  ' 

READ(*,*)K 

WRITE(*.*)'  ENTER  95%  CONF.  LIMIT  OF  K  ' 

READ(*,*)ERRK    | 

GOTO  10 
1400  WRITE(*.*)'  ENTER  K2  ' 

READ(*.*)K2 

WRITE(*,*)'  ENTER  95%  CONF.  LIMIT  IN  K2  ' 

READ(*,*)ERRK2 

GOTO  598 
700   STOP 

END 

SUBROUTINE  LLSQ (T , D . ERRD , NPOINT . SLOPE . ERRS . FLAG . YINT ) 
C  THIS  SUBROUTINE  CALCULATES  THE  SLOPE  AND  ITS  95% 
CONFIDENCE  LIMITS. 
C  EITHER  ESTIMATED  OR  "T  TEST"  GENERATED  ERR  Y  IS  USED. 

DIMENSION  D(IOO) ,T(100) , ERRD (100) 

INTEGER  FLAG 

CHARACTER*1  YESNO 

DIMENSION  CT(31) 

DOUBLE  PRECISION 
SUMX . SUMY . SUMXY . SUMYS , SUMXS . BOTTOM . SIGMAY , SUMRS 

DATA 

CT(1).CT(2) ,CT(3) ,CT{4) ,CT(5)/12.706.4.303.3.182.2.776. 
42.571/ 

DATA 

CT(6) .CT(7) .CT(8) .CT(9) ,CT(1 0)72.44  7,2. 365. 2. 306. 2. 262, 
&2.228/ 

DATA 

CT(11).CT{12),CT(13),CT(14),CT(15)/2.201.2.179.2.160.2.145. 
42.131/ 

DATA 

CT(16) ,CT(17),CT(18) ,CT(19) .CT(20)/2.120.2.110,2.101,2.093. 
42.086/ 

DATA 

CT(21) ,CT(22)  ,CT(23)  .CT(24)  ,CT(25)/2.080.2.074.2.069.2.064. 
42  .  060/ 

DATA 

CT(26).CT(27).CT(28) ,CT(29) ,CT(30)/2.056,2.052,2.048.2.045. 
42.042/ 

IP(FLAG.EQ.2)G0T0  600 


w 
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IF(FLAG.EQ.3)G0T0  100 

WRITE(*,*)'  0=ZERO  INTERCEPT      1=N0N-ZER0  INTERCEPT 

READ(*.*)N 
IF(N.EQ.O)GOTO  100 
C   ***«««***«*«NON  ZERO  INTERCEPT  LEAST  SQUARES********** 
600   SUMX=0 
SUMY=0 
SUMXY=0 
SUMXS=0 
SUMYS=0 
SUMRS=0 
DO  870  1=1 .NPGINT 

SUMX  =  SUMX  +  T(I  ) 

SUMY=SUMY+D(I) 

SUMXS=SUMXS+T(I )*T(I ) 

SUMYS=SUMYS+D(I)*D(I) 

SUMXY=SUMXY+D(I)*T(I) 
870   CONTINUE 

B0TT0M=FL0AT(NP0INT)*SUMXS-SUMX*SUMX 

SLOPE =( FLOAT (NPO I NT )*SUMXY-SUMX* SUMY) /BOTTOM 

YINT=(SUMY*SUMXS-SUMXY*SUMX) /BOTTOM 

DO  900  1=1 ,NPOINT 

SUMRS=SUMRS+( ( SLOPE*T ( I ) +YINT) -D ( I ) ) * ( { SL0PE*T ( I ) +YINT ) 

-D(I)) 
900   CONTINUE 

WRITE(*,*)  'DO  YOU  WISH  TO  USE  THE  ESTIMATED  ERRY? 
YES=1  ' 

READ(*.*)NANS 

IF(NANS .EQ. 1  )  GOTO  910 

SIGMAY=SUMRS/(NP0INT-2) 

SIGMAY=SQRT(SIGMAY) 

ERRY=CT(NP0INT-2)*SIGMAY/SQRT(FL0AT(NP0INT) ) 

ERRS=ERRY*(SQRT(NPOINT/BOTTOM) ) 

RETURN 
910   EDMAX=ERRD(1) 

DO  930  I=2,NP0INT 
930   EDMAX=MAX(ERRD{I) .EDMAX) 

ERRY=EDMAX 

ERRS  =  ERRY* SQRT( FLOAT (NPO I  NT) /BOTTOM) 
RETURN 
100   WRITE(*.*)'  ****ZERO  INTERCEPT  LEAST  SQUARES****  ' 
SUMXY=0 
YINT=0.0 
SUMXS=0 
SUMRS=0 

DO  200  I=1,NP0INT 
SUMXY=SUMXY+D(I )*T(I) 
SUMXS=SUMXS+T( I )*T(I) 
SLOPE=SUMXY/SUMXS 
SUMRS=SUMRS+(SLOPE*T(I)-D(I) ) * ( SLOPE*T ( I ) -D ( I ) ) 
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200   CONTINUE 

WRITE(*.*)'  DO  YOU  WISH  TO  USE  ESTIMATED  ERRY?  YES=1  ' 
READ (*,*) NANS 
IF  (NANS.EQ.l)  GOTO  300 
SIGMAY=SUMRS/FL0AT(NP0INT-1) 
SIGMAY=SQRT(SIGMAY) 

ERRY=CT(NP0INT-1 ) *SIGMAY/SQRT ( FLOAT ( NPO INT ) ) 
WRITE(*,*) 'ERRY  IS  ' 
WRITE(*,*)ERRY 
ERRS=ERRY/SQRT(SUMXS) 
RETURN 
300   WRITE(*,*)'   DETERMINING  THE  BARATRON  FACTOR?  YES=1  ' 
READ(*,*)NANS 
IF  (NANS.EQ.l)  GOTO  500 
EDMAX=ERRD(1) 
DO  400  I=2.NP0INT 
400   EDMAX=MAX(ERRD(I) .EDMAX) 
ERRY=EDMAX 

ERRS=ERRY/SQRT(SUMXS) 
RETURN 
500   WRITE(*,*) 'DEFAULT  ERROR  IN  BARATRON  IS  2E-6. 
CHANGE=Y' 
88    FORMAT(A)         ' 
READ(*.88)YESN0 
IF(YESNO.EQ. 'Y' )THEN 

WRITE(*,*)'  INPUT  THE  DESIRED  MULTIPLIER  OF  ERROR 
IN  CIN  ' 

READ(*,*)ERRY 
GOTO  515 
ENDIF 
ERRY=2E-6 
C  DEFAULT  ERROR  IN  BARATRON  READING  IS  ESTIMATED  AS  2E-6 
TORR 
515   ERRS=ERRY/SQRT(SUMXS) 
RETURN 
END 

SUBROUTINE  CNIC (NTIM . TEMP , ETEMP . XMAX . EXMAX ) 
C  THIS  SUBROUTINE  MAKES  BASELINE  CORRECTION  AND  FINDS  ERROR 
BARS  OF  Y  VALUES  FOR  BOTH  ZERO  AND  NONZERO  INTERCEPT  LLSQ 
CALC. 

DIMENSION   TEMP(IOO) ,ETEMP(100) 
CHARACTER*1  YESNO 

WRITE(*.*) 'DEFAULT  ERROR  MULTIPLIER  IN  CIN  IS  0.1 
CHANGE=Y' 

READ(*.88)YESN0 
IF(YESNO.EQ. 'Y' )THEN 

WRITE(*,*)'  INPUT  THE  DESIRED  MULTIPLIER  OF  ERROR 
IN  CIN  ' 
c.  READ(*  ,*)ECIN 

^J:  GOTO  915 

K  ENDIF 

p:,.  ECIN  =  0.1 


V  i 


915 
500 


600 


625 


700 


750 


900 
ETEMP 


910 
950 
960 

970 
88 
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DO  500  I=1,NTIM 

ETEMP(I)=ECIN*TEMP(I) 

XMIN=TEMP(1) 

DO  600  I=2,NTIM 

XMIN=MIN(TEMP(I) ,XMIN) 

EXMIN=ECIN*XMIN 

XMAX=TEMP(1)      1 

DO  625  I=2,NTIM  ' 

XMAX=MAX ( TEMP ( I ) , XMAX ) 

EXMAX=ECIN*XMAX  , 

DO  700  1=1 .NTIM 

TEMPd  )=TEMP(I  )-XMIN 

DO  750  1  =  1, NTIM 

ETEMP(I)=SQRT(ETEMP(I )*ETEMP(I ) +EXMIN*EXMI N ) 

WRITE(*,*)I 

WRITE(*,*)TEMP(I  )  .ETEMP(I) 

CONTINUE 

WRITE(*.*)  'LLSQ  CALC:   NON-ZERO  INTERCEPT=1 

READ (*,*) NANS 

IF(NANS.EQ. 1 )  GOTO  950 

XMAX=XMAX-XMIN 

EXMAX=SQRT(EXMAX*EXMAX+EXMIN*EXMIN) 
DO  900  1=1 , (NTIM-1 ) 


( I )=SQRT(ETEMP( I ) *ETEMP(I ) / (TEMP( I ) *TEMP( I ) ) +EXMAX* 
&EXMAX/(XMAX*XMAX) )*TEMP(I ) 
DO  910  1=1 , (NTIM-1) 
TEMP(I)=TEMP(I)/XMAX 
DO  960  1  =  1 ,  (NTIM-1  ) 
ETEMP(I)=ETEMP(I)/TEMP(I ) 
DO  970  1=1 , (NTIM-1 ) 
TEMP(I)=LOG(TEMP(I) ) 
FORMAT(A) 
RETURN 
END 


ZERO=0 


YR.NDP) 

for  cga  adaptor. 

OBJ  which  contains 


SUBROUTINE  GRAF2D(XR 
C  General  plotting  routine 
C  Must  be  linked  with  GRAF 
language 

C  subroutines  WRITDOT,  and  SETDIS. 
C  Written  by  Cliff  Watson  9-17-1987. 

DIMENSION  XR(1024)  ,YR(1024) 

CHARACTER  YESN0*1 

COMMON  /SCALE/  XMIN . XMAX , XS , YMIN . YMAX . YS 

YMAX  =  YR( 1 ) 

YMIN  =  YR(1) 

XMIN  =  XR(1 ) 

XMAX  =  XR(1) 


the  assembly 


C 
AND 


FIND  XAW 


MIN  VALUES 
DO  10  I  = 


1  .NDP 
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10 


IF  (XR(I)  .GT.  XMAX)  XMAX  =  XR(I) 

IF  (XR(I)  .LT.  XMIN)  XMIN  =  XR(I) 

IF  (YR(I)  .LT.  YMIN)  YMIN  =  YR ( I) 

IP  (YR(I)  .GT.  YMAX)  YMAX  =  YR (  I  ) 

CONTINUE 

WR I TE  (*.*)  XMAX,  XMIN,  YMAX,  YMi;^ 

PAUSE 


SET  DISPLAY  TO  HIGH 


SET  SCALE 

XS  =  639  /  (XMAX  -  XMIN) 
YS  =  199  /  (YMAX  -  YMIN) 
C      ENTER  GRAPHICS  MODE 

RESOLUTION  GRAPHICS 

CALL  SETDIS(6)   ' 
CALL  BOX(0,0,639 , 199) 


PLOT  DATA  I 

DO  20  I  =  l.NDP 

IX  =  (XR(I )  -  XMIN)  *  XS 
lY  =  (YR(I)  -  YMIN)  *  YS 
CALL  PSET(IX,IY,1 , lERR) 
CALL  BOX( IX-5, IY-5, IX+5. IY+5) 
20    CONTINUE 

C  I  PAUSE 

TO  DISPLAY  PLOT 

READ(*. • (A) • )  YESNO 
C 
DISPLAY  TO  TEXT 

CALL  SETDIS(2) 

RETURN 

END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
SUBROUTINE  PSET ( IX , I Y , I  COLOR ,  TERR ) 

C  Subroutine  to  set  pixel  on  high  resolution  graphics 
screen . 

C  X  =  0  ,  639 ;   y  =  0 , 199 
lERR  =  0 


SET 


85^' 


IP( (IX.LT. 0) .OR. (IX.GT.639) )  lERR  =  1 

IF( (lY.LT. 0) .OR. ( lY.GT. 199) )  lERR  =  1 

IF  ( (ICOLOR.LT.-l  )  .OR.  ( ICOLOR.GT. 1 )  )  lERR  =  1 

IF( lERR. EQ. 1 )  GO  TO  100 

INVY  =  199  -  lY 


W: 
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CALL  WRIDOT(IX,INVY,ICOLOR) 
100   CONTINUE 
RETURN 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  LRL(IXMIN, IXMAX. lY.IERR) 
C  Subroutine  to  draw  horizontal  line 
DO  100  I  =  IXMIN. IXMAX 
CALL  PSET(I , lY.l , lERR) 
100   CONTINUE 
RETURN 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  LUD ( I YBOT , I YTOP , IX , lERR ) 
C  Subroutine  to  draw  vertical  line 
DO  100  I  =  lYBOT.IYTOP 
CALL  PSETdX.I  ,1  ,  lERR) 
100   CONTINUE  ' 

RETURN 
END 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  BOX (IXMIN, lYMIN. IXMAX. lYMAX) 
C  Subroutine  to  draw  box 

C  Draw  top  of  box 

CALL  LRLdXMIN,  IXMAX.  lYMAX.  lERR) 
C  Draw  bottom 


'  ^t^ 
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CALL    LRLdXMIN.IXMAX.IYMIN.IERR) 
C   Draw    left    side 

CALL  LUD(IYMIN,IYMAX,IXMIN,IERR) 

C  Draw  right  side 

! 

CALL  LUDdYMIN,  lYMAX.  IXMAX,  lERR) 

RETURN 

END 


APPENDIX  II 

ANALYTICAL  EXPRESSIONS  FOR  KINETIC  MODELING  IN  CHAPTER  5 

A.   Kinetic  Differential  Rate  Equations  for  the  Scheme  Shown 
in  Figure  5.3. 


I      =1°      -° 

S«3'     S«3'     ~  -3-3 


+  1°     [^  +   6<j./(ioe. )]  [export)  -  1] 


I     ^  =  31- 


_C  „  ^    [->^P(92^)  -  ^-P  (^^H  -  To  e*^°     ^  {[exp(9^t) 


£-C3H3 


-i]/e^  -  3[exp(0  t)  -  i]/e  } 


=  31 


C  HD  ^   "a-c   H  +f^"^^^^^  -  ^expO^t)  +  exp(e3t)]  + 


Yo6*i°  {[exp(e  t)  -  i]/e    -  6[exp(e  t)  -  il/e 

£-0    H  '  2  J'     : 


S°3  ^-<=3"3 


+  6[exp(e3t)   -  ij/e^} 

+  [3  export)    +  expO    t)    -  exp(e.t)    -    3  exp(e,t)l 


-   6<t.I  +([exp{0^t)    -    l]/(10e    )    -   9[exp(9    t)    - 

^"^3^3 


ij/doe^)  +  gfexpce^t)  -  i]/(503)  -  [expo  t)  -  i]/e  |, 


u^- 
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m  Which  5  =  Vk,,  ^  =  k,)CfPc^,^/(k^  .  k,),   e,  =  k,^/io  -  k,Pc^,^  -  ^^^a,u,i, 

92  =  3k^4./10  -  kfP^^^^   -  kpP^^  63  =  3k^^/5   -  kfP  -  kpPc  „  I.  and   64   = 


k^((.  -  kfP^^p^  -  kpP(,^,j^j..     Theflttlng   parameters    for    this    model    are 


s 


l-C   H  ■"'  ^'  ^f^C2D2'and  Vc3H3l-   1°   +  is  fixed  at  it 
3  3  C3H3 

experimental  value  If  the  Initial  point  Is  at  t=0;  otherwise 
it  is  also  a  fitting  parameter. 


B.   Kinetic  Differential  Rate  Equations  for  the  Scheme  Shown 
In  Figure  5.7. 


dl 


^dt  =   81 


^-^3"3'  ^-^3"3 


dl  /dt  =   BI  /(lOF^) 

^-S«3^  ^-^3«3 


dl 


/dt   =  -    (f\'    +  k' )    I  +  +   *2^^  + 

£-HDC,H  I  ^         Jl-HDC   H  A-HDC   H 

3  J  J 


FI  )/(5F    )    +    2(^      I  /(5F) 


dl  /dt 


-f(f\;  +  k^)  I       ^  +  *2  (^ 

^         ^""23°  £-HDC   H 


FI  )/(10F^)    +    (t,    I  /(SF^) 


dl 


C-C3H2D 


^/dt    =    3B[4'2(I  +   + 


FI 


)/(2F)    + 


£-HDC    H  £-H   C   D 


hr-.. 


*,! 


3    3 


J/(5F) 
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dl 


£-HDC  D 


/dt   =   -F(F^k'    +  k' )1 
+  r  n 


£-HDC   D  £-HDC    D 


:  )/(5F    )    +    2^    (I  + 

ii-D^C^H  "^       £-HDC    h'^ 


FI  .)/(5F'')    +    ,^1  /{5F) 

i-H_C.D  '    £-C   H, 


2    3 


3    3 


dl 


2     '  ' 

/dt   =    -(F   k,   +  k    )I  +    A    /pT 

£-D   C    H  £-HDC^D 


2   3 


£-HDC   h"*" 


l-D^C^H 


FI  ^)/{5F)    +    *   I  /lo 


^-"2^° 


N-C3H3-' 


dl  /dt 


'^t*3^"  ^   ^    I  )/f2    -H    mi 

£-HDC3D-  ,-o^C3H-  ^2 ^    ^_^^^^^, 


+   FI 


^)/F    +    (j),!  +/2I/5 


dl 


2    3 
.2  .    2     ' 


3  "3 


^-S°2 


+/dt    =    [<|,^/F      -   F(F   k^    +  k    )]l  + 


£-0303 


2.i>3(FI 


£-HDC,D 


+    I 


£-D    C   H^'^^'"    '    "    ^2^'  .   ^   "  J/(10F) 

^  °2S"  £-HDC   H^  £-H   r   n+ 

-"  23 
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dl 


c-C   D 
3    3 


^/dt  =    8[,|>^I  ^/f2    +   24,3  (FI  ^ 


l-C   D 
3    3 


£-HDC.D 


I  _|_)/(5F)    +    <{,    (I  .     +   FI  .)/10, 

£-D   C    H  £-HDC    H  Z-HCD 

^    -^  3  2    3 


.2,  • 


k      =  k   P_  ..   _,    (|)^    =  F'k  / 


in   which  6  =  <!>  /(lOF^)    -  F^k'    -  k',   k'  =  k^P^  ^    ,   ^      -  -  -„  ,9 

1  f         P       f  f  C^D^'     p         p  C3H3I      *1 

(A^    +    B\^),    <^2   =   F^V(X3   +    6A^),    ,(,3    =  f\^/(X^    +    SX^),    <J)^    =   f\^/(1    +   Fg), 
and    the    ^'s    are    given    as: 


X^    =    (1    +    6/F  +    3/F    )/10 
X      =    (3   +    6/F   +    1/F^)/10 


X3   =   X^/F 

X^   =    (3   +   2/F)/5F^) 

X^    =    (2   +   3/F)/(5F). 


o  '    I 

The  fitting  parameters  for  this  model  are  ^      .»  B'  k  ,  k  ,  and  P 

£-C  H  ^  P 
1°     is  fixed  if  the  initial  point  is  at  t  =  0 ;  ^otherwise  it  is 

C  H 

also  a  fitting  parameter. 


C.   Kinetic  Differential  Rate  Equations  for  the  Scheme  Shown 
in  Figure  5.10. 


3  3       3  3      *  ^3"3 


^.^  u  +  "  ''-11*1^       [exp(e  t)  -  exp{9.t)]/( 
^5"3  ^-"^2^3 


12 


I- 


'*r- 


■:Tf!! 
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^C  H   ■*■  ~  ''-11^-21*1*2^°  +[exp(e^t)/e^3    -   exp(e2t)/923   + 


S"3 


^-S"3 


(1/623    "    1/9i3)exp(93t)]/e^ 


'c  H   ^s)    =  '-1l'sl*1  V°       „  ,{[exp(e^t)    -   l]/e,    -    [exp(e_t)    - 


S«5    ''' 


^-S"3 


iJ/ezl/Qi: 


'    H  +fo  ~  ''-ii''-2i''s2*i*2V°         ^{[exp(e  t)  -  i]/(e  e    )  - 
[expce^t)  -  ^Vie^Q^^)  -  (1/823  -  V6^3)[exp(e  t)  - 


l]/93}/9i2' 


H,+  =  ''-''''p2'^^l_^  ^  Jexp(e^t)/e;  -  expce^tve^"  .  d/e^ 


8    6 


•i/e^)exp(-k2t)]/e^2' 
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in  which  ^r^,,^cn/iKu    -  K,2   -   ^-13^'  *2  =  ^f2^C  H  /(^-21  -  '^-22  "^ 

^    ■^  4  2 

'^sr'  *3  "  ^fS^C  H  ^^^-3^   "^  ^_32  +  ^s2''  ~^1 '  "®2'  ^"'^  "®3  ^'^^  parameters 
(1),  (2),  and  (3),  in  that  order,  in  equations  (a).  9,2  ~  ^i  "  ^2'  ®13  '' 

9l  -  ,S'  ^23  =  ^2  -  ^3'  ^1  =  ^1  +  ^d'    ^""^   ^2  =  ^2  ""  ""d*   The  fitting 
parameters  are  9^.  e^,  63,  l°     ^(1  +k_3^^^/e^),  k_^^*^l°     ^/e^^' 


i-c^n. 


^s2*3'  ^^-21*2'  ^51*2'  and  k^.   It"  can  b 


^-S"3 


e  seen  that  k^^  =  -9i  "  ^^-21*2  " 


^sl'^2'  i"  which  k 


P2    ^p2Pc3H3l* 


;::ti-7^ 


»  ■ 
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